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Calculations  also  show  that  the  energy  barrier  to  intramolecular  conversion 
from  the  linear  to  the  cyclic  form  is  prohibitively  high.  The  linear  form  can 
react  with  to  yield  four  possible  isomers  of  CjHj+. 

The  section  on  theoretical  studies,  which  includes  results  of  quantum 
mechanical  calculations,  also  presents  a  new  algorithm  for  the  geometry 
optimization  searches.  The  algorithm  is  designed  to  reduce  the  ntnnber  of 
iterations  required  to  find  the  lowest-energy  geometry  for  a  stable  molecule 
and  for  transition-state  structures  in  chemical  reactions.  This  alogarithm  is 
used  in  the  energy  calculations  presented  in  this  report. 

Quantum  mechanical  calculations  to  determine  ultiaviolet  and  Iniared 
spectral  line  positions  are  presented  for  some  isomers.  These  calculations, 
were  done  to  determine  which  laser  frequencies  are  required  to  fragment  the 
various  isomers.  Laser  studies  to  determine  which  Isomers  might  be  present, 
using  the  results  of  these  calculations,  were  done.  However,  these  studies 
did  not  successfully  determine  which  isomers  of  the  various  cations  were 
present  in  the  ICR  reaction  chamber. 

ICRMS  studies  reveal  that  linear  C3B3'*‘  reacts  with  C4H2  to  yield 
higher  molecular  weight  cations,  but  C3H3^  and  C2H2  do  not  react  to 
yield  detectable  levels  of  higher  molecular  weight  cations.  However,  linear 
C3H3'*’  does  isomerize  to  cyclic  C3H3''‘  and  in  the  presence  of  C2D2 
yields  an  isotopic  mixture  of  cyclic  C3H3'*’.  Therefore,  some  isomer  of 
C5H5'*‘  is  forming,  but  apparently  falls  apart  rapidly.  The  reaction  of 
C3H3'^  and  C4H2  leads  to  the  formation  of  C5H3'*‘, 

C9H5'^,  and  all  of  which  are  detected  in  the  reaction  chamber. 

Kinetic  models  presented  in  the  report  show  the  mechanistic  possibilities 
for  the  reactions  of  linear  C3H3'''  with  C2H2,  C2D2  and  C4H2. 

The  model  ^ich  best  fits  the  data  for  the  reaction  of  linear  C3H3'*'  with 
C2D2  incorporates  a  kinetic  isotope  effect  factor.  This  factor  indicates 
there  is  some  bonding  rearrangement  in  the  traxisition-state  structure  leading 
to  the  formation  of  C5H5'*’.  The  model  idiich  best  fits  the  data  for  the 
reaction  of  linear  C3fl3^  with  C4H2  Incorporates  collision-induced 
stablizatlon  steps  for  the  C9H5''*  and  Cj^iHj'*'  cations.  These  are  the 
only  two  cations  in  the  scheme  showing  a  build-up  in  concentrations  not 
followed  by  an  attenuation. 

Reactions  of  C2H2  and  C4H2  with  C5H3+,  C5H5+  and  C7H7+ 
cations  formed  from  various  neutral  precursors  were  also  studied  as  a  part  of 
this  work.  Reactions  with  C4H2  are  generally  faster  than  reactions  with 
C2H2.  However,  the  degree  of  reactivity  of  the  various  cations  depends 
upon  what  the  neutral  precursors  are.  Apparently,  different  proportions  of 
the  more  reactive  isomers  are  formed  from  different  precursors.  Threshold 
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Studies  presented  in  this  report  show  that,  at  low  pressures,  cationic 
reactions  with  C4H2  are  more  Important  than  reactions  with  C2H2. 

However,  at  flame  pressures,  collision- induced  stablizatlons  may  make 
reactions  with  C2H2  more  Important. 
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SUMMARY 


An  understanding  of  the  mechanisms  that  cause  soot  formation  in  flames  and 
other  combustion  environments  could  help  to  eliminate  this  formation  and  lead 
to  more  efficient  fuel  consumption,  lower  atmospheric  pollution,  and  even 
eliminate  plumes  seeded  through  exhaust. 

Two  soot  formation  mechanisms  have  been  postulated.  One  invokes  the 
formation  of  small  radicals,  that  further  react  with  small  hydrocarbons  to 
form  the  larger  particles  that  make  up  soot.  A  second  mechanism  Involves  the 
reaction  of  small  positive  ions  with  neutral  hydrocarbons,  again  leading  to 
soot  through  condensation  and  condensation/elimination  reactions.  This 
mechanism  is  posttilated  due  to  the  observation  of  high  concentrations  of 
ions  in  sooting  flames  and  the  precursors  of  this  mechanism  have  been  examined 
in  the  work  reported  here. 

The  small  hydrocarbon  ions  found  in  sooting  flames  have  been  created  in  a 
Fourier  Transform  Ion  Cyclotron  Resonance  (FTICR)  mass  spectrometer,  and  the 
formation  of  larger  hydrocarbon  ions  has  been  examined  in  a  fashion  more 
systematic  than  can  be  obtained  in  flames  themselves.  It  is  found  that  there 
are  at  least  two  forms  of  C2H3^  that  might  be  present  in  sooting  environments, 
and  that  one  of  these  forms,  probably  the  cyclic  isomer,  is  not  reactive 
toward  either  acetylene  or  diacetylene.  Through  deuterium  exchange  studies, 
at  least  one  form,  the  linear  propargylium  cation,  has  been  shown  to  react 

with  acetylene,  as  well  as  diacetylene.  Studies  of  the  rates  of  these 

reactions  have  been  undertaken  to  assess  their  role  in  actual  sooting 

situations.  Reactions  of  C5H5''’,  and  CyHy'*’  have  also  been  studied  under 

the  terms  of  this  contract,  indicating  the  possible  course  sooting  reactions 
might  follow  in  the  formation  of  larger  molecules  responsible  for  soot. 

Quantum  chemical  methods  have  been  developed  for  studying  these  molecules, 
and  calculations  have  been  performed  on  the  lowest  energy  isomers  of  and 

Four  stable  structures  of  and  nine  of  have  been  examined 

in  detail,  and  predictions  made  of  their  vibrational  and  electronic  spectra 
for  possible  identification  in  sooting  flames  and  exhausts.  Their 

Identification  would  Indeed  lend  credence  to  the  ionic  mechanism  of  soot 

I 

I 
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formation.  Reactions  between  the  isomers  of  and  acetylene  have  been 
f'xamined  theoretically,  and  it  was  found  that  the  cyclic  form  does  not  react, 
whereas  the  linear  form  does.  In  addition,  it  was  found  that  only  one  of  the 
reactions  between  the  propargyl  ion  and  acetylene  leads  to  a  product 
thermodynamically  stable  relative  to  disproportionation  back  into  the  cyclic 
form  of  and  acetylene.  These  findings  are  in  good  agreement  with  those 
from  the  FTICR  experiments,  and  have  allowed  the  mechanisms  that  might  be 
operative  in  soot  formation  to  be  modeled. 
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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

The  purpose  of  this  research  was  to  provide  information  about  the 
reactivities,  structures,  and  spectra  of  gaseous  ions.  This  information 
allowed  an  assessment  of  the  proposed  ionic  mechanism  of  soot  formation,  in 
part,  testing  its  validity  when  compared  to  various  schemes  which  rely  on  free 
radical  mechanisms.  A  detailed  understanding  of  soot  nucleation  can  guide 
soot  abatement  efforts  in  many  combustion  systems,  leading  to  improved 
efficiency  in  burners  and  jet  engines  and  to  reduced  particulate  emissions 
(and  thus  lowered  visibility  and  signatures)  from  various  combustors. 

B .  BACKGROUND 

1.  Ions  in  Flames 

It  has  long  been  known  that  significant  concentrations  of  ions  exist 
in  flames.  For  many  years  their  presence  was  considered  a  scientific 
curiosity.  More  recently  the  practical  importance  of  flame  ions  has  been 
recognized  (1).  The  presence  of  high  ionization  levels  in  jet  engine  exhausts 
and  rocket  plumes  can  serve  both  to  impede  radio  communication  and  to  provide 
unwanted  exhaust  signatures.  Intentionally  seeded,  conducting  flames  have  a 
high  potential  for  practical  use  in  MHD  reactors.  Essentially  all  hydrocarbon 
flames  are  characterized  by  a  predictable  and  constant  level  of  ions;  this 
fact  forms  the  basis  for  the  nearly  universal  hydrocarbon  flame  detector  used 
in  modem  gas  chromatography.  There  has  been  a  controversy  for  many  years 
involving  the  importance  of  ions  in  flames  to  soot  formation.  Almost  equally 
compelling  evidence  has  been  given  for  the  primary  importance  of  radical- 
molecule  reactions  and  for  ion-molecule  reactions  in  the  condensation 
reactions  which  lead  to  nucleation  and  ultimately  to  soot  formation  (2). 

For  these  and  other  reasons  there  has  recently  been  rapid  development 
years  of  techniques  for  the  study  of  ions  and  electrons  in  high  temperatures 


1 


environments.  A  variety  of  techniques  have  been  employed.  These  include 
microwave  absorption  and  phase  shift  techniques  for  measuring  electron 
densities.  For  more  sensitive  point  source  measurements  in  the  relatively 
thin  flame  zones  Langmuir  probes  have  been  extensively  used.  Stark  broadening 
has  been  used  for  measurements  of  ionization  in  seeded  flames.  These 
techniques  are  sensitive  to  local  concentrations  of  charge  but  provide 
relatively  little  Information  concerning  the  nature  of  the  ions  or  their 
formation  or  destruction  mechanisms. 

Combustion  products  present  in  flames  all  have  relatively  high  ionization 
potentials.  Therefore,  the  presence  of  high  ion  concentrations  far  exceeds 
that  expected  from  thermochemical  equilibrium.  It  is  generally  accepted  that 
primary  ions  are  created  by  chemi-ionization  at  a  very  thin  reactive  zone  in 
flames,  undergo  extensive  secondary  reactions  and  ultimately  evolve  primarily 
to  NO"^  and  electrons  in  the  burnt  gas  region. 

At  the  present  time  the  most  sensitive  probing  of  ion  concentrations  in 
flames  is  done  with  flame-ion  mass  spectrometers.  In  this  technique 
microprobes  are  used  to  sample  gases  in  the  reactive  zones  of  flames  (2-10). 
Typically  the  pressure  is  reduced  to  expand  the  volume  of  the  reactive  zone  as 
an  aid  in  probing.  Both  positive  and  negative  ion  concentrations  above  10^ 
ions  cm  have  been  characterized  in  this  way  in  a  variety  of  flames.  A 
serious  potential  exists  for  sampling  artifacts  when  using  this  technique 
because  of  catalyzed  reactions  in  the  hot  sampling  probe  and  in  the  mass 
spectrometer.  These  effects  are  discussed  extensively  by  Hayhurst  and 
Kittelson  (6). 

It  is  apparent  from  References  2-10  that  positive  ion  concentrations  in 
flames  are  significant,  with  C3H3'*’  and  H30‘*'  reaching  concentrations  near  10^® 
ions*cm~^  in  two  torr  flames.  Concentrations  of  negative  ions  are  believed  to 
be  on  the  order  of  100  times  lower,  with  electrons  contributing  the  remainder 
of  the  negative  charge  required  for  overall  neutrality. 

Positive  ion  concentrations  are  relatively  insensitive  to  flame 
temperature  (11).  This  is  consistent  with  the  nonthermochemical  equilibrium 
nature  of  their  production  mechanism  and  reflects  the  probability  that 
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secondary  ions  are  generated  from  a  small  number  of  ionic  precursors  which  are 
created  by  chemi-ionization.  Studies  of  the  concentration  dependence  of 
several  positive  ions  upon  fuel-oxidizer  stoichiometry  (4)  give  results  which 
again  are  consistent  with  chemi-ionization  production  of  one  or  two  primary 
ions  followed  by  a  large  number  of  ion-molecule  react  ions - 

2.  Possible  Formation  Mechanisms 

It  is  generally  agreed  that  primary  ions  in  hydrocarbon  flames  arise 
from  chemi-ionization  rather  than  from  pyrolysis  although  evidence  for  a 
pyrolysis  origin  continues  to  surface  (5).  A  score  of  different  chemi- 
ionization  reactions  have  been  proposed  and  investigated  for  hydrocarbon 
flames  (2)  as  potential  primary  ion  formation  mechanisms.  Most  workers  are 
currently  in  agreement  that  two  radical-molecule  reactions  are  dominantly 
responsible  for  creation  of  primary  ions: 

CH(X^n)  +  0(^P)  — HCO'*’  +  e~  (1) 

CH(A^A)  +  -^>  +  e".  (2) 

Contrary  to  what  is  reported  in  the  older  literature,  Reaction  (1)  is 
currently  thought  to  be  slightly  exoergic  (7).  This  is  consistent  with  the 
ion  formation  process  being  nearly  temperature  independent.  Considerable 
controversy  abounds  concerning  the  rate  of  Reaction  (1).  Values  range  from 
10“^^  to  10~^^  cm^ •molecule”^ •sec“^ .  The  most  recent  careful  measurement 
under  flow  tube  conditions  by  Vinckier  (12)  gives  a  value  of  k^  =  2.4  x 
cm-^ ‘molecule  ‘sec  ^  at  room  temperature.  It  is  generally  concluded  that  the 
reaction  may  be  a  factor  of  ten  faster  at  flame  temperatures.  The  rather 
inefficient  channel  for  ion  formation  competes  with  the  reaction  to  produce 
neutrals: 

CH(X^n)  +  0(^P)  — —>  CO  +  H  (3) 

OH  +  C 

In  laser  experiments  Fllseth  (13),  ^  al. ,  have  followed  the  total 

disappearance  rate  of  CH(X^n)  and  report  a  value  for  k3  =  8.8  X  10“^^ 
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cm^ •molecule'^ ‘sec”^ .  It  is  apparent  that  the  branching  between  Reactions  (1) 
and  (3)  favors  formation  of  neutrals  by  10^  -  10"^ . 

Reaction  (2)  is  about  25  kcal/mole  exothermic  as  written  and  would 

therefore  be  very  endothermic  for  the  reaction  of  ground  state  CH(X^II)  to 
p-^oduce  C3H3'*’.  Although  the  Importance  of  Reaction  (2)  in  flames  has  been 
questioned  (2),  it  is  currently  the  only  chemi-ionization  reaction  seriously 
considered  to  be  responsible  for  formation  of  C3H2"*’,  which  precedes  the 
formation  of  HCO'*'  in  some  flames  and  therefore  is  considered  a  primary 

ion.  In  fact,  it  has  been  used  as  the  primary  ion  formation  mechanism  in 
recent  crude  models  of  sooting  flames  (14)  although  the  modelers  admit  to  many 
unanswered  questions  about  its  importance.  Of  course,  can  also  be 

produced  by  secondary  reactions,  i.e., 

CHs'*'  +  (4) 

or  more  complicated  schemes  (7).  C3H3'*’  appears  very  early  in  all  of  Goodings 

and  Bohme's  (7)  profiles  and  the  concentration  profiles  do  not  show  shoulders 
or  secondary  peaks. 

3.  Secondary  Ions  and  Soot  Formation 

Goodings  and  Bohme  (7)  have  mapped  out  a  consistent  flow  picture  for 
the  formation  of  secondary  ions  in  both  fuel-rich  and  lean  conditions  for 
methane-oxygen  flames.  Large  numbers  of  reactions  involving  HCO'*’  with 
neutrals  dominate,  which  create  larger  intermediate  (usually  oxygen 

containing)  ions.  The  ultimate  fate  of  most  of  these  ions,  particularly  in 
fuel-lean  conditions,  involves  formation  of  HCO'*',  e.g. 

+  OH  - HCO'*’  +  +  H  (5) 

followed  by  conversion  to  H30^  via  the  reaction 

HCo'*’  +  H^O  - H^O"^  +  CO,  (6) 
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where  Che  equilibrium  lies  to  the  right. 

Certain  secondary  ion-molecule  reactions  in  flames  may  be  important  in 
understanding  soot  formation.  Studies  by  Prado  and  Howard  (15)  and  by  others 
(9,10)  have  demonstrated  the  presence  of  large  hydrocarbon  ions  (probably 
polyaromatic)  in  sooting  flames.  Experiments  involving  the  effects  of  ionic 
additives  (16)  and  of  electric  fields  (17)  on  the  amount  of  soot  produced  in 
such  flames  also  indirectly  point  to  participation  of  ions  in  these 
processes.  Whether  these  large  ions  are  forced  directly  from  smaller  ions  by 
ion-molecule  reactions,  or  whether  they  are  formed  by  charge,  proton,  or 

hydride  transfer  after  large  hydrocarbon  neutrals  have  been  formed  remains 
unclear.  In  recent  papers  discussing  ions  in  flames  Hayhurst  and  Kittelson 
(6),  without  much  discussion,  conclude  that  the  presence  of  ions  like  C^H^^ 
erd  G^H2'*'  indicates  the  presence  of  large  uncharged  hydrocarbon  fragments  in 
the  burnt  gas  region,  Goodings,  Bohme,  and  Ng  (7)  argue  that  both  neutral  and 
ion-molecule  condensation  mechanisms  may  lead  to  the  larger  charged  particles 
important  in  soot  formation,  but  do  not  assess  the  relative  importance  of  the 
two  competing  pathways.  Bayes  and  co-workers  (18)  state  rather  unequivocally 
that  "the  rapid  ionic  polymerization  of  acetylene  ...  is  a  probable 

Initiator  of  soot  formation  in  rich  acetylene  flames."  In  his  Plenary  Lecture 
(19)  on  "Soot  Formation  in  Combustion",  at  the  17th  International  Combustion 

Symposium,  H.  G.  Wagner  considered  the  matter  still  unresolved:  "the  rapid 

growth  of  the  large  hydrocarbon  molecules  which  are  necessary  intermediates  of 
soot,  indicated  that  their  chemical  growth  must  proceed  via  fast  reactions. 
These  can  be  either  radical  reactions  or  ion  reactions,  or  a  combination." 

Recent  investigations  of  fuel-rich  and  sooting  ethylene  (20),  acetylene 
(9,10,20),  and  benzene  (10)  flames  have  provided  evidence  which  further  fuels 
the  ionic  versus  neutral  soot  formation  controversy.  Delfau,  Michaud,  and 
Barassin  (20)  found  evidence  for  two  independent  ionization  regions  in  futl- 
rich  ethylene  and  acetylene  flames,  the  first  corresponding  to  small  chemi- 
ions,  and  the  second  to  larger  ions  at  the  early  stages  of  soot  formation. 
They  argue  that  the  large  ions  can  be  formed  by  direct  thermoionization  of 
heavy  hydrocarbon  molecules  and  young  soot  particles,  which  are  formed  by 
neutral,  and  not  ionic  mechanisms.  They  further  believe  that  the 
concentration  of  high  mass  ions  is  too  large  to  be  accounted  for  by  an  lon- 
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molecule  pathway  from  the  relatively  low  concentrations  of  small  chemi-ions 
they  observe.  The  same  workers  (9)  have  observed  C2H3'*’,  C5H3'*’, 

^9^3^’  large  polycyclic  aromatic  hydrocarbon  (PCAH)  ions  in  sooting, 

premixed  acetylene  low  pressure  flames,  and  have  explained  the 
2,  3,  and  4)  concentrations  in  terms  of  equilibria  between  these  ions,  the 

dominant  C3H3'*’  ion,  and  C2H2  and  H2  neutrals  present  in  the  flame. 

Based  on  their  own  measurements  (10)  and  a  different  interpretation  of  the 
above  mentioned  work,  the  group  at  Aerochem  Research  Laboratories  has 
developed  an  ionic  model  for  soot  formation  (14a, b).  They  also  see  C2H3^  as  a 
major  ion  in  sooting  flames,  plus  relatively  high  concentrations  of  C5H3'*’, 

CyHy"^,  0^2Hg^,  ^13^9"^>  ^15^11"^’  and  ^i9^1i^  ions.  They  argue  for 
sequential  ion-molecule  reactions  Involving  the  smaller  ions  and  neutrals  such 
as  the  C2H2  fuel  as  well  as  C4H2  other  carbon  rich,  hydrogen  poor  species 

formed  in  relatively  high  abundance  in  the  flame.  Even  with  assumed  ion 
molecule  reaction  rates  of  only  2  X  10“^^  cm^ *3”^  for  all  relevant  reactions, 
their  admittedly  crude  and  incomplete  model  shows  that  adequate  concentrations 
of  large  ions  which  may  be  direct  soot  precursors,  as  well  as  reasonable 

levels  of  neutrals  formed  by  ion  neutralization,  can  be  predicted.  However, 

they  repeatedly  stress  the  inadequacy  of  the  input  data  such  as  ion-molecule 
reaction  rates,  ion  structures,  and  ion  thermochemistry  used  in  their  model. 

4.  Recent  Developments  in  Neutral  Soot  Formation  Mechanisms 

One  of  the  strongest  arguments  for  an  ionic  soot  nucleation  mechanism 
has  been  that  it  is  difficult  in  any  of  the  proposed  neutral  schemes  to 
produce  small  cyclic  and  polycyclic  species  with  a  sufficiently  high  rate  to 
model  the  quite  rapid  soot  formation  seen  in  flames  (14b).  lon/molecule 
reactions  are  known  to  proceed  quite  rapidly,  and  to  yield  cyclic  products  if 
those  structures  are  the  most  stable.  Thus  ion/molecule  reactions  seem  almost 
necessary  in  the  early  stages  of  soot  nucleation.  Recently,  however,  a 
neutral/radical  soot  formation  model  has  been  proposed  which  contains  many 
steps  and  incorporates  quite  recently  obtained  shock  tube  kinetic  data  (21). 
This  model  is  capable  of  yielding  small  ring  compounds  at  a  higher  rate  than 
previously  suggested  radical  mechanisms,  although  there  is  some  question  as  to 
the  correct  rate  constant  for  one  or  two  of  the  crucial  steps  in  formation  of 
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the  first  ring.  In  comparing  an  ionic  model  with  this  radical  scheme,  it  is 
quite  important  to  be  able  to  know  where  the  first  ionic  ring  compound  is 
formed.  Is  a  f ive-membered  ring  formed  as  one  of  the  isomers,  or  must 

the  system  progress  to  the  Cy  or  Cg  size  in  order  to  form  six-membered  or 
larger  rings?  One  recommendation  of  a  recent  workshop  on  mechanisms  of  soot 
formation  (22)  was  that  improved  data  for  hydrocarbon  ion  reactions  and 
structures  be  obtained  so  that  a  more  complete  ion/molecule  model  might  be 
developed  and  compared  with  the  latest  neutral/free  radical  model. 

C.  SCOPE /APPROACH 

This  project  involved  combined  experimental  and  theoretical  study  of 
possible  ionic  precursors  of  soot  in  hydrocarbon  combustion  systems.  The 
reactivities,  structures,  and  spectra  of  a  number  of  small  hydrocarbon  ions 
which  have  been  postulated  as  important  in  soot  formation  were  investigated. 
The  ions  studied  have  all  been  observed  in  fuel-rich  and  sooting  flames  and 
included  C3H3'^,  C5H3’^,  C5H5'^  and  CyKy'*’. 

Experimental  studies  were  carried  out  using  Fourier  Transform  Ion 
Cyclotron  Resonance  Mass  Spectrometry,  and  concentrated  on  kinetic  and 
spectroscopic  characterization  of  the  gaseous  ions.  Rate  coefficients  for 
ion/molecule  reactions  of  the  ions  listed  above  with  a  range  of  neutral 
molecules  found  in  flames  were  determined.  Measurements  were  generally  made 
at  quite  low  pressures  and  near  room  temperature,  but  temperatures  up  to  ca. 
550  K  and  pressures  up  to  ca.  5  X  10”^  torr  were  also  investigated. 
Electronic  spectra  for  ions  were  probed  by  laser  photodissociation 

techniques,  using  a  Nd:YAG  pumped  dye  laser. 

Theoretical  studies  were  carried  out  by  Dr.  Zerner  and  his  research  group 
using  facilities  of  the  University  of  Florida  Quantum  Theory  Project.  These 
Included  a  VAX  11/780  computer,  a  SUN  System  composed  of  3/50  work  stations 
and  3/280  File  and  Compute  servers,  and  a  variety  of  resident  programs  for 
calculating  energetics  and  spectra  of  chemical  species.  Including  gaseous 
ions.  Supercomputer  facilities  at  Florida  State  University  and  the  University 
of  Pittsburgh  were  also  used.  The  energies  and  properties  of  several  isomers 
of  C3H3^  and  ^585^  were  calculated  at  successively  higher  levels  of 
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was  also  followed 


sophistication,  and  reactivity  of  with  ^2^2 

theoretically.  Spectra  of  the  and  C5H5'*’  ions  in  both  the  IR  and  UV-VIS 

regions  were  predicted. 
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SECTION  II 


THEORETICAL  STUDIES 


A.  INTRODUCTION 


Theoretical  studies  were  performed  on  the  possible  precursors  of  soot  and 
the  elementary  reactions  between  small  carbocations  and  acetylene.  Quantum 
mechanical  methods  of  the  semiempirical  and  ab  initio  types  were  used,  and 

these  are  described  in  some  detail  in  the  sections  that  follow. 

In  the  early  stages  of  this  study  it  was  realized  that  truly  effective 
ways  of  calculating  molecular  geometry  from  an  initial  configuration  of  atomic 
positions  and  the  number  of  electrons  were  not  available.  Methods  were 

designed  to  accomplish  this  under  the  terms  of  this  contract,  and  these 
methods  are  now  in  widespread  use  by  many  groups  as  a  result  of  our 
findings.  The  first  two  subsections  of  this  section  deal  with  the  methods 
that  we  have  uncovered.  They  are  basically  reprints  of  papers  published  in 
part,  or  in  toto.  A  review  paper  has  also  been  written  on  this  subject  by 
John  Head  and  Michael  Zerner  and  will  soon  appear  in  the  ^Vivances  in  Quantum 
Chemistry  (Academic  Press),  a  book  which  publishes  only  invited  contributions. 

Calculations  for  this  contract  were  performed  initially  on  the  Quantum 

Theory  Project's  VAX  11/780,  augmented  with  the  3-RA81  disk  subsystem  bought 
with  funds  for  this  work.  As  we  discovered  that  only  the  ab-initio  big  basis 
set  calculations  were  really  reliable  in  predicting  the  relative  energies  of 
the  species  in  which  we  were  Interested,  some  of  these  calculations  were 
performed  on  the  Pittsburgh  Supercomputing  Center  facilities,  and  some  through 
a  separate  grant  for  this  purpose  on  the  Florida  State  University  Cyber  205. 
In  addition,  supplementary  funds  obtained  during  this  contract  allowed  the 
purchase  of  a  SUN  3/280  computing  system,  which  was  totally  dedicated  to  these 
calculations  during  the  final  period  of  this  contract.  We  had  grossly 
underestimated  the  amount  of  computing  resources  that  would  be  required  to 
lead  to  definitive  results  in  the  relative  stabilities  of  these  systems,  and 
to  characterizing  these  systems  with  respect  to  their  IR  and  UV-visible 

spectroscopy.  By  devoting  supercomputer  resources  to  this  work,  in  addition 
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to  the  Quantum  Theory  Project's  considerable  resources,  however,  we  have 
obtained  state-of-the-art  results,  and  have  even  shown  (Section  IIP)  that  many 
of  our  colleagues  reporting  on  these  structures  in  the  past  have  been  wrong. 

As  will  be  presented  in  section  HE,  we  have  begun  an  examination  of  the 
actual  reaction  of  acetylene  with  C2H2'*’-  Our  results  are  in  agreement  with  the 
FT-ICR  experimental  results  presented  In  Section  III.  This  work  will  carry  on, 
even  in  the  absence  of  further  funding  on  this  project,  for  we  believe  that  such 
studies  will  lead  to  a  clearer  understanding  of  the  possible  role  that  these 
small  carbocations  have  in  the  formation  of  soot  particles.  Other  studies  that 
should  be  undertaken  Involve  the  reactions  of  with  diacetylene,  and  of 
other  small  ions  with  both  acetylene  and  diacetylene. 

B.  A  BROYDEN-FLETCHER-GOLDFARB-SHANNO  OPTIMIZATION  PROCEDURE  FOR  SEARCHING 

QUANTUM  MECHANICAL  SURFACES  FOR  MOLECULAR  GEOMETRIES 

1.  Introduction 

The  use  of  molecular  quantum  mechanics  to  obtain  molecular  geometry  has 
become  reasonably  widespread.  The  structures  of  molecules  as  yet  unmade  and  of 
excited  state  of  molecules  present  few  addtional  problems.  Transition  states 
that  can  only  be  inferred  from  experiments  are  somewhat  more  difficult  to  obtain 
but  utilize  the  same  general  techniques. 

The  most  common  techniques  used  to  search  a  molecular  energy  hypersurface 
obtained  either  through  semi-empirical  or  ab  initio  quantum  mechanics  require 
analytic  first  derivatives  of  the  energy  with  respect  to  coordinate  displace¬ 
ment,  and  then  proceed  to  build  up  curvature  information  on  the  surface  as  they 
move  along  ir.  The  most  common  of  these  methods  used  are  of  the  Murtagh  Sarrent 
(MS)  (23)  type  first  introduced  into  quantum  chemistry  by  Mclver  and  Komornicki 
(24).  For  small  molecules,  this  method  has  proven  quite  successful.  For  larger 
molecules,  gradients  are  generally  reduced  to  below  10"^  Hartree/Bohr  with 
reasonable  efficiency,  but  a  further  monotonlc  decrease  of  the  gradients  to  an 
acceptable  value  is  quite  slow.  It  is  this  reduction  difficulty  that  has 
prompted  us  to  reexamine  the  various  update  procedures  that  are  available  and  to 
test  their  performance,  especially  with  respect  to  the  popular  MS  method.  In 
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doing  this  ve  have  found  one  method,  a  variant  of  the  Broyden-Fletcher-Godlfarb- 
Shanno  (BFGS)  algorithm  (25-29)  is  systematically  superior,  especially  so  for 
large  molecules,  and  does  proceed  to  reduce  the  gradients  in  a  more  or  less 
monotonic  fashion.  Ue  describe  this  method  here  and  compare  it  with  our  version 
of  the  MS  method.  A  more  complete  report  on  these  methods  as  well  as  others  ve 
have  examined  is  in  progress  (30). 


2.  Quasi-Newton  Method 

Ic  k. 

The  energy  E(x  )  and  the  gradient  vector  g(x  )  are  expanded  in  a  Taylor 
series  about  x*^  yielding 


k+1.  n/  k  k  1  k  „k  k 

E(x  )  =  E(x  )*g  s  H  s 


k+1  /  k.  t,k  k 

g  3  g(x  )  +  H  s 


(7a) 

(7b) 


with 


s 


k 


k+1 

X 


(7c) 


k.  k  * 

and  where  H  =  H(x  )  is  the  matrix  of  second  derivatives,  of  the  "Hessian" 

matrix.  Extreme  points,  such  as  minima  saddle  points  or  maxima,  are 

k+1 

characterized  by  g  =  0.  Assuming  that  g  ~  2  large  enough  k  and 

truncating  Eq.  (7b)  after  the  second  term,  yields  the  familiar  Newton  equation 

k 

for  the  search  direction,  s 


k 

s 


■  1  k 
g 


(8) 


Jr 

For  an  exact  H,  s  is  said  to  be  a  true  "Newton"  step  toward  the  extremum.  A 
sequence  of  such  steps  determines  a  Newton  procedure. 


Since  the  calculation  of  H  eithe’-  analytically  or  numerically  is  a  process 
many  times  slower  than  that  for  obtaining  the  energy  and  gradients.  Many 
methods  can  be  used  to  approximate  H  the  most  common  of  these  are  the  quasi- 
Newton  methods.  Since  most  optimization  algorithms  repeatedly  solve  Eq.  (8)  for 

Ir 

s  ,  the  inverse  Hessian  is  usually  sought  directly .  An  analogous  development 
for  the  Hessian  itself  is  easily  developed. 
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The  quasi-Nevton  methods  relate  the  differences  of  gradients  obtained  at  two 
different  points  on  the  potential  energy  hypersurface 


k  k  +  1  k 
Y  =  g  -  g 


(9) 


to  the  differences  in  the  coordinates 


jk  k  k  k+1  k 
0  =  a  s  =  X  -X 


(10) 


which,  for  a  quadratic  function,  should  obey  the  relation 

Y  =  Ho 


(11) 


Quasi-Newton  methods  are  constrained  to  obey  the  "Quasi-Newton"  condition 
obtained  from  rearranging  Eq.  (11). 


„k+l  k  ,k 
G  Y  =  6 


(12a) 


(12b) 


The  Murtagh-Sargent  method  is  a  Rank  1  method  defined  by 

„k+l  „k  „k  ^k  I  1 

G  “  G  c  =  G  +  s 


(13) 


where  we  have  used  the  common  Dirac  notation  and  a  is  a  simple  scale  factor. 

Use  of  the  Eq.  (12a)  yields 

g‘^|y’^>  +  hIuXuIyS  =  IsS  (14) 

k  -1 

Choosing  the  constant  a  in  Eq.  (14)  as  <u | y  >  and  substituting  into  Eq.  (13) 
yields  the  update  equation  (25-28) 


„k+l  „k 
G  =  G  -t- 


(|5'^>  -  g'^I  y’^>)(<5'^1  -  <y’^|g'^) 

k .  j.k  r>k  k. 

<Y  I  6  -  G  Y  > 


(15) 


that  forms  the  basis  of  all  Rank  1  methods. 


The  BEGS  method  is  a  Rank  2  method 


G*^^^  =  G*^  +  aluXul  +  blvXvl 


(16) 


Various  choices  of  | u>  and  | v>  and  of  a  and  b  distinguish  these  methods.  The 
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BFGS  choice  (29)  leads  to 


pk+1  k 

-BFGS  ” 


(<6*^1  <y'^|G*^|  y‘^>)  I  sS<&'^| 

ixrk,  k  1 2 

I  <6  I Y  >  I 


5*^><y*^|g'^  +  G*^|y'^><5*^| 


(17) 


If  cartesian  coordinates  are  used,  G^  is  usually  set  to  the  unit  matrix,^. 
The  first  step,  s^,  is  then  in  the  direction  of  steepest  descent.  Subsequent 
cycles  update  G  toward  greater  accuracy  as  they  move  along  the  surface. 


The  Rank  1  formula  has  two  principal  disadvantages.  The  positive 
definiteness  of  G  may  be  lost  on  updating.  Secondly,  the  denominator  of  Eq. 
(15)  may  approach  zero  sufficiently  rapidly  as  to  lead  to  numerical  inaccuracy. 
Murtagh  and  Sargent  have  devised  simple  methods  to  detect  these  problems.  There 
seems  to  be  no  really  successful  solution  to  either  of  these  problems  when  they 
do  occur,  however,  and  the  inverse  Hessian  matrix  is  reset  to  the  unit  matrix. 


In  contrast,  the  BFGS  formula  ensures  that  the  Hessian  remains  positive 

definite  on  all  updates.  A  positive  definite  Hessian  is  desirable  because  it 

Ic 

always  generates  a  search  direction  s  ,  given  by  Eq.  (18)  below,  which  is  in  the 
direction  of  energy  decrease.  The  BFGS  formula  can  even  be  advantageous  over 
the  exact  Hessian  H  in  regions  where  the  exact  Hessian  has  negative  eigenvalues. 
A  search  direction  will  have  uphill  steps  along  any  normal  mode  of  the  molecule 
which  has  a  negative  eigenvalue. 


3.  Line  Searches 


With  G  updated  from  either  the  MS  or  BFGS  formula,  a  search  direction. 


s  ,  is  determined  by 


k  ^k  k 
s  =  -2  S 


(18) 


A  line  search  is  then  performed  to  determine  how  far  along  this  direction  a 
step,  6*^,  should  proceed. 
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(19a) 


k+1  k  k  k 
X  =  X  +  a  s 


jk  k  k 
6  s  a  s 


(19b) 


k+1  k 

A  weak  line  search  requires  only  that  E  <  E  if  searching  for  a  minimum.  An 

Ic  k,  Ic  k. 

exact  line  search  finds  that  value  of  a  for  which  E(x  +  a  s  )  is  a  minimum. 

there  is  thus  a  competition  between  the  number  of  energy  evaluations  that  are 

k  k  k  k 

required  to  minimize  E(x  +  a  s  )  along  a  given  s  ,  or  accepting  nearly  any  move 

along  s  that  reduces  the  energy  or  the  gradient  norm  and  then  finding  a  new 
k+1 

direction  s  at  the  new  geometry. 

Ue  have  adopted  the  following  line  search  strategy  (30).  In  the  k-th  line 

k  k  k 

search  cycle  we  examine  the  interval  (a^,  of  acceptable  a  values. 

k  k  ^  ^ 

Initially  =  0  and  0-2  is  determined  by 

I  <A  =  0.4a.u.,  ^20) 

a  condition  preventing  arbitrarily  large  coordinate  moves.  The  value  of  0.4 
a.u.  is  itself,  of  course,  arbitrary,  and  larger  values  might  be  acceptable  for 

larger  molecules.  The  initial  a  value  is  set  to  1  except  for  the  first  cycle 

1  1 
where  a  is  set  to  0.4.  In  this  first  cycle  G  =1  and  the  move  is  along  the 

direction  of  steepest  descent.  The  value  =  0.04  is  empirical  and  nearly 

always  satisfies  the  line  search  conditions  even  for  poor  starting  geometries. 


k  k  k  k 

If  the  initial  a  >  a2  then  a  is  reset  to  0.2- 

evaluated  and 


„k+l  c,,  k  k  k.  . 

E  =  E(x  +  a  s  )  IS 


„k+l  .  k  .  k,  k. 
E  >-pa  <slg> 


k  k 

checked.  <s  |g  >  is  the  first  term  in  the  Taylor  series  expansion  for  the 
energy  and  is  negative  for  a  decrease  in  energy.  From  Eq.  (7)  and  Eq.  (8)  0  <  p 
<  1/2,  with  p  =  1/2  corresponding  to  an  exact  line  search  and  any  positive  value 
simply  leading  to  a  decrease  in  energy.  After  considerable  numerical  testing, 

we  choose  p  =  0.01  corresponding  to  a  "weak"  search.  If  this  test  fails, 

k  k  k  k  k  k  • 

however,  a  >  a2  (true).  a2  ^  a  (a  is  "replaced"  by)  and  a  is  reduced  using 

the  information  already  available  according  to 
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k  ,  k  1/  k  k. 

a  -*  OL  +  j{OL 


,  K  k,  1  k. 
(a  -  ap  E  (a^) 


^1/  kx  /  k  k  k-+  k 
E  (a^)  =  g(x  +  a^s  )  s  = 


/'  \ 

3a 

\  y 


(22a) 


(22b) 


ot=a, 


k+1  k  k  k 

The  energy  E  ->  E(x  +  as)  is  now  computed  and  the  right  test  of  Eq. 
(21)  repeated  until  the  test  is  obeyed. 


k.  1  Ic  k  Ic 

The  gradient  g  =  g(x  as)  is  now  evaluated  and  a  left  extreme  test 


applied 


<g*^^^|s*^>|  <  -cr  <g*^|s*^>,  0  <  a  <  1 


(23) 


k+1  k 

For  an  exact  line  search  <g  js  >  =  0  (the  gradient  is  perpendicular  to  the 

search  direction)  corresponding  to  a  =  0.  The  value  0=1  only  requires  a 

k+1  k 

reduction  of  the  scalar  product  <g  |s  >.  For  the  test  of  Eq.  (23)  we  choose  o 
=  0.9,  again  corresponding  to  a  weak  line  search.  If  the  left  test  fails, 

however,  a  is  too  small,  and  a  new  value  is  estimated  from 


.  ,  (a*^  -  a*[)  E^(a*^) 

a‘^(new)  =  a  +  — r - j — j— 

E^a^)  -  E^fa*") 


(24) 


k  k  k  k 

The  energy  E  -♦  E(x  +  a  (new)  s  )  is  evaluated  and  both  right  and  left  line 

k  k  k  k 

test  repeated  with  a^  -»  a  and  a  -»  a  (new).  When  both  tests  have  been 
successfully  passed,  the  partial  line  search  is  finished  and  a  new  search 
direction  is  sought. 


In  the  MS  procedure  the  test  for  a  positive  definite  Hessian  may  fail.  If 
so  G  is  reset  to  the  unit  matrix  and  the  line  search  is  conducted  as  if  this 
were  once  again  the  first  cycle. 

4.  Optimization  Algorithms 

An  initial  geometry  x.,  is  chosen.  E^  and  g^  are  calculated.  Then 
(i)  A  search  direction  s  is  chosen  through  Eq.  (18).  Then 
(ii)  A  line  search  is  performed  to  determine  how  far  along  this  direction 
to  move  as  described  in  previous  section.  Then 
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(iii)  A  test  for  convergence  is  made.  If  convergence  is  satisfactory,  the 
procedure  is  complete.  If  not, 

(iv)  Estimate  a  new  Hessian  or  its  inverse  and  repeat  the  procedure. 

The  success  of  a  procedure  is  measured  in  the  number  of  such  cycles  that  are 
required  to  find  the  extreme  point  sought.  The  fever  energy  evaluations,  the 
better  the  method. 

5.  Results  and  Discussion 

Table  1  summarizes  results  of  geometry  optimizations  comparing  the  MS 
and  BFGS  methods.  The  notation  x(y,2)  that  we  use  indicates  that  x  line 
searches,  y  energy  evaluations  and  z  gradient  evaluations  were  required, 
starting  with  geometries  ranging  from  quite  good  to  very  poor.  An  indication  of 
the  quality  of  the  starting  geometry  is  given  by  the  maximum  value  of  the 
gradient.  Values  greater  than  0.1  Hartree/Bohr  (1  Hartree  =  627.5  kcal/mol) 
usually  indicate  one  or  more  bond  lengths  are  displace  by  more  than  0.1  A.  Two 
of  the  larger  molecules  of  this  study  were  not  successfully  optimized  by  the  MS 
procedure.  After  fifty  cycles,  for  example,  the  largest  gradient  of 
( (CH2)2P(0)0 I ”  is  still  larger  than  2  x  10”^  a. u. /Bohr  and  that  for  CH2CHCHO 
greater  than  7  x  10”^  a. u. /Bohr.  The  former  case  is  difficult  as  the  two  methyl 
groups  are  rotating  relative  to  one  another  in  an  eclipsed  conformation  along 
two  very  soft  modes. 

Table  2  shows  how  an  optimization  proceeds  for  the  Rank  1  update  using  the 
line  search  described.  Line  4  of  the  Hessian  eigenvalues  indicates  that  a 
negative  eigenvalue  of  the  Hessian  has  been  detected.  This  negative  eigenvalue 

4 

causes  the  step  s  to  go  uphill  on  the  energy  surface  even  though  <g|s>  is 

negative.  There  are,  as  a  consequence,  repeated  failures  of  the  right  extreme 

line  search  test  Eq.  (21)  and  after  the  eighth  energy  evaluation  the  line  search 

terminates  without  a  proper  a  value.  At  this  stage  another  Hessian  update  is 

4  8 

taken  using  the  gradients  evaluated  at  the  coordinates  x  and  x  and  fortunately 
the  new  Hessian  generated  is  positive  definite.  In  the  remaining  cycles  the 
Hessian  remains  positive  definite  and  the  NH.,  geometry  optimizes  smoothly.  We 

^  _4 

note  that  even  after  all  gradients  have  been  reduced  below  10  a. u. /Bohr,  the 
approximate  Hessian  eigenvalues  are  still  not  accurate. 
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TABLE  1.  SOME  EXAMPLES  COMPARING  THE  MS  AND  BEGS  UPDATE  SCHEMES. 


NH^'^ 

max  |g^ 

0.112 

NH3'^ 

0.198 

H2O2 

0.024 

*"2^6 

0.193 

H2CO 

0.264 

H^CCHCHO 

0.105 

[(CH3)2P(0)0]~ 

0.119 

Pyrimidine 

C4^2«4 

0.127 

[2,2'  ]  para 
cyclophane 

Cl6»16 

0.050 

MS 

BFGS 

7(9,8) 

7(8,8) 

7(11,8) 

7(8,8) 

5(5,5) 

6(6,6) 

17(19,18) 

19(20,20) 

9(11,10) 

11(12,12) 

NO® 

24(27,26) 

NO® 

43(44,44) 

11(14,13) 

10(11,11) 

16(20,18)^ 

10(11,11) 

The  notation  x(y,z)  refers  to  x  line  searches,  y  energy  evaluations  and  z 
gradient  evaluations. 

^The  largest  initial  gradient  in  a. u. /Bohr  giving  a  measure.of  thg  "poorness"  of 
the  initial  geometry.  Convergence  is  assumed  when  max  |g.  |<10~  a. u. /Bohr, 
except  for  [ (CHt),P(0)0]~,  pyrimidine  and  [2,2')  paracyclophane  where  max 
g^  I <10~  a.u./Bonr. 

Q 

Displaced  0.05  a.u.  along  each  normal  mode. 

*^Displaced  0.08  a.u.  along  each  normal  mode. 

0 

Not  optimized  after  50  energy  evaluations, 
f  k  -3 

For  paracyclophane  max|g  | <10  a.u. /Bohr  at  14(16,15)  but  the  subsequent  cycle 
has  an  increase  in  gradient  above  this  threshold  even  though  the  energy 
monotonically  decreases. 
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TABLE  2.  RANK  1  (MS)  UPDATE  OPTIMIZATIONOF  NH^-HESSIAN  EIGENVALUES  (A.U.). 


1 

I.OOOIHOO 

1.000I>H)0 

l.OOOD+OO 

l.OOOD+OO 

1.0000400 

1.0000400 

2 

2.72A1>+0U 

1.00UI>+0U 

I.OUUIHUU 

1.0UU04UU 

1.0Ul)l>400 

1.UUU04UU 

3 

3.227IHOO 

2.71AIHOO 

I.OOOIHOO 

1.0000400 

1.0000400 

1.0000400 

A  -2.35AD-02 

2.717IHOO 

1.8020^00 

1.0000400 

1.0000400 

1.0000400 

6 

2.719D+00 

2.331IHOO 

1.458IH00 

1 . 0000400 

1.0000400 

1.689D-01 

9 

2.587IHOO 

2.303IH00 

1.A56IH00 

1.0000400 

9.9940-01 

1. 5810-01 

10 

2.370D+00 

2.267IHOO 

1.4  380400 

9.9990-01 

9.8620-01 

1.5520-01 

11 

2.359I>K)0 

2.2621HOO 

1.4350+00 

9.9970-01 

5.3020-01 

1.506D-01 

Exact 

2 

.27AAIHOO 

2.27AAIHOO 

1.45190400 

1.59980-01 

1. 54300-01 

1.54300-01 

Eval 

E 

.  No. 
G 

Energy 

Ul 

Mnx  1  g J 

a 

1 

1 

-12.499372 

0.358164 

0.198326 

O.OOOOOD400 

2 

2 

-12.521284 

0.050917 

0.031968 

4.00000D-01 

3 

3 

-12.522080 

0.028010 

0.017822 

1.00000D400 

4 

4 

-12.322281 

0.019552 

0.012914 

1.00000I>400 

5 

-12.476193 

8.67976D-01 

6 

-12.520752 

8.018330-02 

7 

-12.522073 

1. 551100-02  - 

8 

5 

-12.522235 

0.020010 

0.012987 

3.67105D-03 

9 

6 

-12.523207 

0.004431 

0.002909 

1.000000400 

10 

7 

-12.523216 

0.000612 

0.000343 

1.000000400 

11 

8 

-12.523216 

0.000086 

0.000049 

1.000000400 
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The  MS  procedure  would  detect  the  negative  eigenvalue  at  cycle  4  and  reset 

4 

the  Hessian  matrix  to  the  unit  matrix  and  a  =  0.4.  In  spi  fa  nf  rhe  nf 

previous  Hessian  information  by  the  MS  procedure,  one  less  cycle  is  required  to 

complete  the  optimization.  In  [ (CH^l^PfOlO J  the  Hessian  is  reset  17  times  in 

50  cycles,  and  there  seems  to  be  no  more  successful  procedure  for  controlling 

the  Rank  1  update  than  that  suggested  by  Murtagh  and  Sargent.  A  somewhat  more 

successful  procedure  of  the  MS  type,  however,  does  result  by  setting  the  right 

k, 

hand  extreme,  =  1;  that  is,  never  permitting  an  a  value  greater  than  1. 

Table  3  examines  this  same  situation  for  the  BFGS  method.  As  expected  all 
the  eigenvalues  of  the  approximate  Hessian  remain  positive  on  update,  and  the 
optimization  proceeds  rather  smoothly.  At  the  second  cycle  two  eigenvalues  of 
the  Hessian  are  updated,  a  characteristic  of  the  Rank  2  nature  of  this  method, 
to  be  compared  to  the  one  value  updated  in  the  Rank  1  schemes.  In  spite  of 
this,  at  convergence  the  Hessian  eigenvalues  obtained  from  the  BFGS  scheme  are 
no  better  than  those  obtained  from  the  MS  procedure. 

For  larger  molecules  both  the  MS  and  the  BFGS  method  seem  equally  able  to 

_2 

reduce  the  gradients  to  below  about  10  a. u. /Bohr,  usually  corresponding  to 

bond  lengths  to  within  ±0.02  A  and  bond  angles  ±2°-3°  of  their  final  values. 

Dihedral  angles,  however,  usually  to  do  with  soft  torsional  or  umbrella-like 

modes,  are  not  accurate  with  this  level  of  convergence.  The  BFGS  method  is 

superior,  as  indicated  in  Table  4  for  H2CCHCHO.  MSI  is  the  procedure  described 

in  this  work.  MS2  is  a  Murtagh  Sargent  procedure  in  which  a2  =  1  briefly 

discussed  above.  This  table  is  typical  in  showing  the  erratic  behavior  of  the 

gradient  norm  (g|  with  the  cycling  for  the  MS  procedure.  To  cycle  10  the 

methods  have  reasonably  similar  behavior.  However,  after  that,  the  BFGS  method 

-4 

is  clearly  superior.  At  cycle  20  the  BFGS  method  is  10  a.u.  above  the 

-4 

converged  energy,  whereas  both  MS  algorithms  are  greater  than  10  a.u.  from 
energy  convergence  after  50  energy  evaluations. 

In  summary,  for  small  molecules,  the  MS  and  BFGS  methods,  as  we  have 
described  them,  have  similar  behavior.  For  all  the  larger  molecules  examined, 
those  reported  here  and  very  many  others,  the  BFGS  method  is  superior.  Both  the 
MS  and  BFGS  methods  seem  equally  able  to  reduce  the  cartesian  gradients  to  below 
0.01  a.u. /Bohr,  but  the  MS  method  often  has  erratic  behavior  for  larger  systems 
after  this  point.  The  BFGS  method  nearly  always  reduces  the  gradient  norm  in  a 
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TABLE  3.  BEGS  OPTIMIZATION  OF  NH^-HESSIAN  EIGENVALUES  (A.U.). 


1 

l.OOOIHOO 

l.OOOIHOO 

l.OOOIHOO 

l.OOOIHOO 

l.OOOD+OO 

l.OOOD+OO 

2 

2.723IHOO 

l.OOOIHOO 

l.OOOIHOO 

l.OOOD+OO 

l.OOOD+OO 

9.784D-01 

3 

3.719IHOO 

1.755D+00 

l.OOOIHOO 

l.OOOD+OO 

l.OOOD+OO 

6.725D-01 

4 

2.836IHOO 

2.378IHOO 

1.032IHO0 

l.OOOD+OO 

l.OOOD+OO 

2.641D-01 

5 

2.785D+00 

2.377IHOO 

1.527IHOO 

l.OOOD+OO 

l.OOOD+OO 

1.618D-01 

6 

2.577IHOO 

2.364IH00 

1.5141HOO 

1 . OOOD+OO 

9.994D-01 

1.541D-01 

7 

2.423IHOO 

2.265IHOO 

1.476IHOO 

l.OOOD+OO 

9.952D-01 

1.627D-01 

9 

2.434I>K)0 

2.347IHOO 

1.472IHOO 

l.OOOD+OO 

9.254D-01 

1.586D-01 

Exact 

2. 

.  2744IHOO 

2. 2744IHOO 

1.4519IHOO 

1.5998D-01 

1.5430D-01 

1.5430D-01 

Eval 

E 

.  No. 
G 

Energy 

Ul 

Max  1  g J 

a 

1 

1 

-12.499372 

0.358164 

0.198326 

O.OOOOOD+00 

2 

2 

-12.521284 

0.050917 

0.031968 

4.00000D-01 

3 

3 

-12.52207  2 

0.028891 

0.018206 

l.OOOOOD+00 

4 

4 

-12.522491 

0.019688 

0.015225 

l.OOOOOD+00 

5 

5 

-12.523027 

0.016294 

0.010447 

l.OOOOOD+OO 

6 

6 

-12.523199 

0.007948 

0.005281 

l.OOOOOD+00 

7 

7 

-12.523215 

0.000941 

0.000621 

l.OOOOOD+OO 

8 

8 

-12.523216 

0.000101 

0.000059 

l.OOOOOD+OO 
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TABLE  4.  A 

COMPARISON  OF 

MS 

AND  BFGS 

UPDATE  METHODS  FOR  H^CCHCHO^, 

Cycle 

AE 

MSI 

Igl 

AE 

MS2 

Isl 

BFGS^ 

AE 

Igl 

1 

0.0123 

0.195 

0.0123 

0.195 

0.0123 

0.195 

5 

2.6D-3 

0.058 

2.8D-3 

0.029 

5.0D-3 

0.028 

10 

1.3D-3 

0.012 

1.3D-3 

9.9D-3 

2.8D-3 

0.013 

20 

5.3D-4 

0.024 

4.6D-4 

6.6D-3 

l.OD-4 

4.0D-3 

30 

2. 7D-4 

1.7D-3 

3.7D-4 

0.028 

l.OD-7 

2.6D-5 

40 

1.8D-4 

1.2D-3 

1.9D-4 

3.4D-3 

50 

1.3D-4 

0.011 

1.4D-4 

2.4D-3 

^Ae  in 

a. 

u.  from  the  converged 

energy  value 

,  |G|  =  <g|g> 

1/2 

in  a. u. /Bohr, 

'^In  32 

cycles  the 

BFGS 

method  has 

reduced  |i 

g|  to  2  X  10" 

6 

a. 

u. /Bohr. 

monotonic 

fashion, 

whereas  this 

is 

seldom 

the  case  for 

the 

MS  method. 

Neither 

update 

method  produces 

a  reliable 

Hessian 

matrix  when 

the 

maximum  gradient  has 

been  reduced  to  10 

-4 

a 

,.  u,  /Bohr 

t 

a  value 

at  which  the  geometries  are 

usually 

reliable. 

If  a  more  detailed 

knowledge  of  the  second 

derivatives  is  required, 

these 

derivatives 

must  be 

evaluated 

analytically 

or 

through  numerical 

differentiation  of  the  first  derivatives. 

C.  A  SURVEY  OF  OPTIMIZATION  PROCEDURES  FOR  STABLE  STRUCTURES  AND  TRANSITION 
STATES 

1.  Introduction 

Theoretical  determination  of  molecular  structure  by  both  empirical  and 
semiempirical  methods  is  becoming  a  powerful  tool  for  the  chemist.  Reliable 
structures  for  stable  molecules  can  be  calculated  as  well  as  geometries  of 
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experimentally  inaccessible  species,  such  as  transition  states.  Derivatives  of 

the  energy  with  respect  to  nuclear  coordinates  play  an  essential  tale  in 

efficient  optimization  algorithms.  This  section  surveys  some  of  these 

optimization  algorithms  which  are  included  in  the  QUIPU  quantum  chemistry 

2 

package  at  the  University  of  Florida^".  However,  these  optimization  algorithms 
are  also  appropriate  for  any  potential  energy  hypersurface  regardless  of  the 
Hamiltonian  used,  providing  the  correct  derivative  expressions  are  available. 
There  is  now  a  large  literature  on  energy  derivatives  and  for  recent  reviews  see 
references  (31-34). 

In  all  optimization  methods  examined  here,  the  first  derivative  of  the 
energy  is  needed.  Second  derivatives  are  also  required  for  some  of  the  schemes 
we  discuss.  To  calculate  the  exact  second  derivative,  however,  is  costly  and 
time-consuming.  If  one  is  interested  in  vibrational  frequencies,  it  is 
important  to  calculate  the  second  derivative  matrix  exactly.  However,  for 
geometry  optimizations,  an  inexact  second  derivative  obtained  by  systematic 
approximations  to  the  coupled  perturbed  Hartree-Fock  equations  is  reliable 
enough  for  successful  application  of  most  algorithms  (30). 

The  next  section  discusses  three  algorithms  for  determining  stable 
geometries.  The  third  section  presents  three  methods  for  transition  state 
optimizations.  The  final  section  gives  some  examples  and  comparisons  of 
application  of  the  methods. 

2.  Minimization  Techniques 

Determining  the  lowest  point  on  a  potential  energy  hypersurface  is 
straightforward.  One  simply  follows  the  decrease  in  energy  until  a  local 
minimum  is  found  with  the  gradient  vector  g  zero.  Verifying  that  the  surface  is 
a  local  minimum  is  also  straightforward.  The  analytical  second  derivative  matrix 
can  be  computed  at  the  appropriate  point  and  tested  for  positive  definiteness. 
Or,  more  simply,  the  geometry  of  a  stationary  point  can  be  distorted  to  give 
another  initial  coordinate  set  for  reoptimizing  the  geometry.  This 
reoptimization  should  return  to  the  same  local  minimum.  A  distortion  which 
lowers  the  symmetry  of  the  molecule  will  be  the  most  effective  test  for  a  true 
local  minimum.  For  full  optimizations,  where  all  the  molecular  coordinates  are 
allowed  to  vary,  there  is  no  real  advantage  in  using  internal  rather  than 
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cartesian  coordinates:  if  the  procedure  is  carried  out  correctly  the  same 

number  of  opfimization  steps  is  required.  For  convenience,  however,  mass 
weighted  coordinates  are  used  and  the  lighter  atoms  move  more  than  do  heavier 
ones.  Thus.  the  optimization  follows  a  -ath  more  in  accord  with  physical 
intuition. 

Three  minimization  algorithms  are  described  below.  Each  algorithm  consists 
of  two  essential  parts:  (1)  selecting  a  search  direction,  and  (2)  determining 
by  a  line  search  a  step  length  to  give  an  acceptable  reduction  in  the  energy. 
An  exact  line  search  finds  the  minimum  of  energy  along  the  chosen  direction; 
partial  line  searches  are  less  demanding,  and  generally  accept  any  reduction  of 
the  energy. 

A  compromise  can  be  reached  between  exact  line  searches  and  accepting  any 
reasonable  reduction  in  energy,  and  then  proceeding  to  estimate  a  new  search 
direction.  In  general,  we  have  found  that  exact  line  searches  are  seldom  worth 
the  extra  energy  calculation  to  obtain  them  (35). 

a.  Newton  and  Quasi-Newton  Methods. 

The  Newton  methods  can  be  derived  from  the  Taylor's  series 

expression 

g(x  +  5)  =  g(x)  4.  g(x)5  (25) 

where  the  minimum  g(x4-5)  =  0  so  that  a  search  direction  s  is  obtained  as 

S  =  as  =  -a(G~^)  g(x)  (26) 

where  G  is  either  the  exact  second  derivative  (Hessian)  matrix  for  the  Newton 
method  or  some  matrix  approximating  the  second  derivatives  for  the  Quasi-Newton 
method.  The  parameter  a  is  determined  by  a  line  search  algorithm  so  that  energy 
exhibits  a  reasonable  reduction  in  value.  Both  the  Newton  and  Quasi-Newton 
methods  use  the  same  line  search  algorithm,  the  details  of  which  are  given 
elsewhere  (35).  When  the  Hessian  is  positive  definite  the  Newton  method  often 
shows  quadratic  convergence  to  a  local  minimum.  Unfortunately  the  Newton  method 
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requires  explicit  calculation  of  the  second  derivative  matrix  which  require 
computational  time  comparable  to  N  Hartree  Fock  Self-Consistent  Field  (SCF) 
calculations,  where  N  is  the  size  of  the  basis  set  used.  Quasi-Newton  methods, 
on  the  other  hand,  only  depend  on  gradient  information  to  reach  the  minimum. 
They  use  the  gradients  calculated  at  the  different  search  geometries  to  build  an 
approximate  inverse  second  derivative  matrix  H  ~  (G  ^).  Using  the  "Quasi- 
Newton"  conditions. 


Y  =  g(x  *  5)  -  g(x) 


(27a) 


Hr  =  S 


(27b) 


as  a  constraint  to  obtain 


H(x  5)  =  H(x)  +  U 


(27c) 


The  search  direction  s  is  then  the  matrix  product 


s  =  -H(x  *  5)g 


Thus  the  approximate  Hessian  is  updated  at  each  cycle.  An  advantage  to  this  is 
that  a  great  deal  of  time  is  not  spent  evaluating  an  accu-'  e  second  derivative 
matrix,  which  is  useful  only  locally  on  the  surface.  Such  detailed  information 
is  only  useful  during  the  course  of  an  optimization  about  the  local  quadratic 
region.  Ve  have  examined  many  such  update  formula  (35)  and  find  the  most 
successful  given  by  Broyden-Fletcher-Goldf arb-Shanno  (BFGS)  (25-28)  as 


=  + 


<Y*^  Ih*^  1  Y*^>: 


!sS<y‘^|H*^  h‘^IyS<6*^| 

Ve  note  from  Eq.  (28)  a  minimum  is  obtained,  s=0  and  g=0,  providing  H(x  +  5)  is 
reasonable.  The  inverse  Hessian,  H,  need  not  be  accurate.  An  important  feature 
of  the  BFGS  formula  is  that  it  retains  positive  definiteness  on  update  if  it  is 
used  with  the  appropriate  line  search  algorithm.  This  positive  definiteness 
feature  implies  that 
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(30) 


<s|g>  =  -<glH|g>  <0 

which  is  the  condition  for  descent  in  energy.  That  is,  the  BFGS  formula  has  the 
property  of  always  providing  a  search  direction  towards  the  local  minimum.  This 
will  happen  even  when  the  true  Hessian  is  indefinite. 


b.  Augmented  Hessian  or  Restricted  Step  Methods. 

This  algorithm  is  useful  for  obtaining  minimum  even  when  the  exact 
second  derivative  is  not  available.  The  search  direction  s  is  found  by  a 
Newton-like  expression 

s  =  -(G  -  XI)‘^g  (31) 

where  X<X(low),  the  lowest  eigenvalue  of  G,  ensures  that  the  matrix  (G  -  XI)  is 
positive  definite.  This  choice  of  X  forces  the  descent  condition  <g|s>  <0  to  be 
satisfied.  The  step  length  |s|  is  a  function  of  X.  A  variety  of  algorithms, 
usually  based  on  the  extent  of  the  quadratic  nature  of  the  potential  energy 
hypersurface  at  the  point  of  interest,  exist  for  selecting  the  maximum  allowed 


In  our  implementation  the  step  length  is  found  by  diagonalizing  the 
"augmented  Hessian" 


fG 


yoLg 


ag" 

=  X 

0. 

(32a) 


2  2 

for  the  lowest  eigenvalue  X  and  lowest  normalized  eigenvector  v  +3  =1.  The 
form  of  the  augmented  Hessian  automatically  requires  X<X(low).  By  expanding  Eq . 
(32a)  two  new  equations  are  obtained 


Gv  +  a6g  =  Xv 
a<g|v>  =  X3 

Rearranging  (32b)  gives  the  step 


s  =  .  -(G  -  XI)-lg 


(32b) 

(32c) 


(33) 


where  a  is  a  parameter  which  can  be  varied  to  give  the  required  step  lengths,  3 
is  defined  by  Eq .  (32c)  and  I  is  the  iinir  matrix.  In  the  d i agonal  1 7a r i on  one  of 
the  following  three  situations  are  obtained. 

1.  A  <  0  and  g  =  0.  In  this  case  v  =  v(low)  is  the  eigenvector  of  G 

corresponding  to  the  lowest  eigenvalue.  <v(low) jg>  is  zero  and  X  = 
X(low).  A  displacement  along  v(low)  'will  reduce  the  energy. 

Alternatively,  if  a  displacement  along  v(low)  is  not  desired,  then 
v(low)  is  shifted  to  a  higher  eigenvalue  by 

G'  =  G  +  X( shi f t ) I v( low)><v( low) I  (34) 

and  the  Augmented  Hessian  reformed. 

2.  X  <  0  and  g  =  1.  G  is  now  positive  definite.  The  step  will  vary  from 

Newton  like  for  small  X  to  steepest  descent  for  large  X.  The  step 

length  will  usually  be  inside  the  restricted  region. 

3.  Intermediate  between  the  two  conditions  is  X  <  0  and  g  ^  0.  G  has  a 

negative  eigenvalue  but  s  given  by  Eq.  (33)  will  give  a  reduction  in 
energy  along  the  negative  eigenvector  mode. 

A  great  deal  of  thought  has  been  given  to  choosing  an  effective  X  for  the 
restricted  step  algorithm  of  Eq.  (33)  (36-38),  but  we  find  the  simple  matrix 
diagonalizat ion  suggested  by  Eq.  (32a)  as  good  as  any  for  finding  minima. 

c.  Norm  of  Gradient  Squared  Methods 


Here  the  sum  of  the  squares  of  the  gradient  is  minimized 


=  <g|g>  =  E  si 


This  method  can  be  applied  to  find  any  stationary  point  and  will  not  necessarily 
find  a  local  minimum  with  respect  to  the  energy;  rather  the  nearest  stationary 
point  is  usually  uncovered.  Proceeding  as  before 


n  -  -  9  V  ^^i 

j  '  dx.  ~  E  g^  (jx. 
3  ^  J 


(36a) 


D  =  2Gg 


(36b) 


where  D  requires  knowledge  of  the  second  derivative  matrix  G.  Second 
derivatives  of  a  require  third  derivatives  of  the  energy.  To  circumvent  this 
requires  a  Newton-like  selection  of  a  search  direction. 


s  =  -G 
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where  G  is  explicitly  calculated  or  obtained  by  the  BFGS  formula,  Eq.  (29).  In 
the  line  search  algorithm  we  require 


CT,  ,  <  CT, 

k+1  k 

3.  Transition  State  Search  Algorithms 


(38) 


The  determination  of  transition  space  structures  is  necessarily  a  more 
difficult  problem  relative  to  local  minimum  optimization.  Here  one  attempts  to 
find  a  maximum  in  the  energy  along  the  reaction  direction  and  minima  for  all 
remaining  directions.  For  a  true  transition  state  structure  the  second 
derivative  matrix  should  have  one  and  only  one  negative  eigenvalue.  The  three 
algorithms  are  then  documented. 

a.  Constrained  Internal  Coordinates 

For  unconstrained  optimization  most  algorithms  can  be  forced  to  be 
invariant  to  coordinate  selection  and,  as  mentioned,  there  is  no  advantage  of 
using  internal  over  cartesian  coordinates.  However,  internal  coordinates  become 
a  valuable  tool  in  constrained  optimizations.  Internal  coordinates  allow 
natural  valence  bond  parameters,  such  as  bond  lengths  or  bond  angles,  to  be 
constrained  as  the  remaining  structural  parameters  are  optimized. 

When  using  internal  coordinates  for  finding  a  transition  state  structure  one 
performs  a  series  of  minimizations  on  a  molecule  with  some  coordinates 
constrained.  A  simple  example  is  the  ring  opening  reaction  of  cyclobutene  to 
cis-butudiene 


One  performs  several  minimization  with  r^^^  held  constant  at  different  distances. 
The  transition  state  corresponds  to  the  maximum  of  the  energy  with  respect  to 
r^2  while  the  energy  is  minimized  with  respect  to  all  the  other  coordinates.  An 


27 


advantage  of  this  approach  is  that  no  second  derivatives  are  needed  to  reach  the 
saddle  point.  The  disadvantage  is  that  one  needs  to  be  able  to  identify  a 
coordinate  important  in  the  reaction, 

b.  Augumented  Hessian 

Several  algorithms  for  walking  uphill,  based  on  the  Augmented 
Hessian  have  been  suggested  (36-38).  Ve  have  implemented  the  procedure  of 
Nguyen  and  Case  (39).  The  search  direction  is  obtained  by  a  similar  expression 
used  for  the  minimum 

s  =  -(G  -  XI)'^g  (39) 
It  is  convenient  to  write  the  step  direction  as 


s  =  -E 


v.Xv.  |g> 
1  1 

X.  -  X 
1 


where  X^  and  |v^>  are  the  eigenvalues  and  eigenvectors  of  G,  respectively.  If 
we  now  choose  X^  <  X  <  X2  the  step  will  go  uphill  along  |vj>  and  downhill  along 
the  remaining  eigenvector  (v^>  modes. 

In  general  the  lowest  eigenvector  |v^>  is  not  the  direction  we  seek  for  an 
uphill  walk.  If,  instead,  |v^>  overlaps  most  strongly  with  the  preferred  uphill 
search  direction,  we  then  scale  X  so  that 

X 

X^  =  X^/A  (41) 
The  search  direction  is  now  given  by 


nVx>  <nv^|g>  ^  |v.><vjg> 

/  X.  -  X 

X  -  X  1 

X  1 


with  the  scale  factor  n  =  2'l X^/X^ .  Ue  use 


X  =  (X2  +  1/2  X^)/2 


to  determine  s.  The  step  is  further  scaled  to  get  a  step  length  appropriate  to 
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the  extent  of  quadratic  region.  Once  the  Hessian  G  develops  a  negative 
eigenvalue  we  use  the  square  of  the  norm  of  the  gradient  to  find  the  stationary 
point . 

c.  Norm  of  the  Gradient  Squared  Methods 

Once  we  are  at  a  point  on  the  potential  energy  hypersurface  where  G 
has  a  negative  eigenvalue  we  minimize 

a  =  <g  |g>  ,  (44) 

as  discussed  above.  This  approach  for  transition  states  was  originally 
suggested  by  Mclver  and  Korminichi  (24). 

4.  Some  Examples 

We  summarize  some  examples  of  minimizations  in  Table  5.  Methods  which 
require  the  calculation  of  the  exact  Hessian  converge  with  fewer  energy 
evaluations,  but  since  each  second  derivative  evaluation  takes  a  factor  of  N 
times  as  long  as  the  SCF  calculation  itself,  (a  factor  of  15  for  cis  1,2 
-bromof lurocyclopropane)  the  approximate  Hessian  methods  are  faster.  Especially 
effective,  even  for  these  small  systems,  is  the  "VI  method"  (30)  that 
approximates  the  coupled  perturbed  Hartree-Fock  solution.  The  Augmented  Hessian 
usually  requires  a  greater  number  of  energy  evaluations  than  does  the  Newton 
method.  The  reason  for  this  is  the  restricted  step  nature  of  the  augmented 
Hessian  (Eq.  (40)).  when  searching  for  a  minimum  X(low)  <  0,  all  other  >  0, 
and  the  step  is  damped. 

As  a  first  example  of  a  transition  state  calculation.  Table  6  presents  a 
search  for  the  transition  state  in  the  reaction. 

H 

A. 

H — C  —  N  C - N  C  —  N — H 

I  II  III 

The  initial  structure  is  the  optimized  linear  structure  of  HCN.  The  direction 
of  search  is  given  as  the  bending  mode.  Nine  searches  of  the  Augmented  Hessian 
type  are  then  required  to  find  a  single  negative  eigenvalue  of  the  Hessian. 
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TABLE  5.  A  SUMMARY  OF  SOME  OPTIMIZATIONS. 

For  each  type  of  optimization  the  number  of  cycles  is  given.  The  number  in 
parenthesis  is  time  (mintsec)  on  a  Vax  11/780,  and  because  of  system,  are  in 
possible  error  of  20%.  The  test  calculations  are  INDO  SCF . 


System 

NH^ 

CH^OH 

cis  l,2-Br,F 
cyclopropane 

Initial^ 

Grad  Norm 

0.169 

0.  160 

0. 143 

BFGS,  Eq .  5 

8 

25  (4:45) 

29  (10:24) 

Minimize  g^g 

Eq.  12 

7 

25  (7:12) 

Newton  Method*^ 

Approx.  Hessian, 

Eq.  2 

6 

4  (1:42) 

9  (  5:46) 

Aug.  Approx.*^ 
Hessian,  Eq.  8 

5 

9  (2:28) 

17  (10:30) 

Newton  Method, 
exact  Hessian, 

Eq.  2 

3 

7  (4:30) 

8  (37:32) 

Aug.  Hessian, 
exact  Hessian 

Eq.  8 

3 

9  (5:28) 

17  (84:42) 

^Initially  all  symmetry  is  broken.  For  convergence  all  gradients  have  been 
reduced  to  below  5*10-4  a.u./bohr. 

^The  approximate  Hessian  is  formed  from  the  V(l)  method. 
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TABLE  6.  HCN  TO  CNH  TRANSITION  STATE  SFAROH. 


Cycle 

Energy  (a.u.) 

Angle 

Grad  Norm 

Max  Grad 

Alpha 

0 

-17.399659 

180. 

0.000009 

0.000007 

0.000000 

1 

-  '  r\n  rs 

—  X/  •  J  y  u  ^ 

1  /  . 

'j.Ui2il4 

0.008997 

O.OOOOoO 

2 

-17.397598 

170. 

0.013596 

0.010633 

0.000000 

3 

-17.396393 

167. 

0.016541 

0.013954 

0.000000 

4 

-17/393372 

162. 

0.022758 

0.018559 

0.000000 

5 

-17.387683 

155. 

0.031371 

0.026375 

0.000000 

6 

-17.377439 

146. 

0.041703 

0.032579 

0.000000 

7 

-17.360093 

134. 

0.054164 

0.044396 

0.000000 

8 

-17.340161 

123. 

0.060774 

0.043167 

0.000000 

9 

-17.318690 

112. 

0.061493 

0.045832 

0.000000 

0 

-17.318690 

112. 

0.061493 

0.045832 

0.000000 

1 

-17.305166 

105 

0.057417 

0.044127 

0.043937 

2 

-17.293763 

98. 

0.048153 

0.037835 

0.136931 

3 

-17.285888 

91. 

0.032596 

0.025713 

0.302938 

4 

-17.282615 

85. 

0.010536 

0.007121 

0.695352 

5 

-17.282563 

83. 

0.000301 

0.000228 

1.000000 

TABLE  7.  CIS 

1,2  BROMOFL'JOROCYCLOPROPANE  TRANSITION  STATE 

SEARCH , 

Cycle 

Energy  (a.u. 

)  C-Br,  C-C 

Grad  Norm 

Max  Grad 

Alpha 

0 

-60.921237 

1.92,1.49 

0.000000 

0.000483 

0.000000 

2 

-60.919410 

1.94,1.49 

0.012906 

0.007870 

0.000000 

4 

-60.910574 

1.98,1.49 

0.030662 

0.019198 

0.000000 

6 

-60.895073 

2.02,1.49 

0.046260 

0.028827 

0.000000 

8 

-60.879683 

2.06,1.49 

0.055757 

0.034111 

0.000000 

10 

-60.861925 

2.07,1.49 

0.063907 

0.037261 

0.000000 

12 

-60.841640 

2.07,1.50 

0.072336 

0.039998 

0.000000 

14 

-60.817903 

2.08,1.51 

0.080903 

0.042538 

0.000000 

16 

-60.791225 

2.09,1.51 

0.088479 

0.044238 

0.000000 

18 

-60.768511 

2.08,1.52 

0.092289 

0.045723 

0.000000 

20 

-60.747962 

2.09,1.53 

0.096211 

0.047915 

0.000000 

21 

-60.735117 

2.09,1.54 

0.096953 

0.048456 

0.000000 

0 

-60.735117 

2.09,1.54 

0.096953 

0.048456 

0.000000 

2 

-60.715634 

2.11,1.55 

0.091909 

0.045909 

0.044725 

4 

-60.724647 

2.11,1.55 

0.076953 

0.038289 

0.088425 

6 

-60.738067 

2.07,1.55 

0.062593 

0.030933 

0.105320 

8 

-60.749658 

2.04,1.57 

0.048845 

0.023988 

0.127181 

10 

-60.758475 

2.02,1.57 

0.036117 

0.017943 

0.154638 

12 

-60.764553 

2.00,1.58 

0.024674 

0.012363 

0.194934 

14 

-60.768260 

1.99,1.59 

0.014560 

0.007292 

0.269270 

16 

-60.770030 

1.97,1.59 

0.005689 

0.002728 

0.469534 

17 

-60.770315 

1.97,1.59 

0.001676 

0.000817 

0.796780 

From 

this  point  five 

searches  are  made 

to  minimize 

(T,  Eq.  (44). 

This  might  be 

compared  with  a  similar  search  of  Baker  on  the 

same  system  that 

required  10 

cycles  to  obtain  the 

transition  state, 

but  starting 

from  a  =  90° 

[AO] 

|. 
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Table  7  summarizes  the  transition  state  search  of  cis 
1 , 2-bromof lurocyclopropane  opening  with  loss  of  Br^, 


A  picture  of  the  transition  state  appears  in  Figure  1.  Starting  with  the 
optimized  structure  of  IV,  cis  1,2-B,F  cyclopropane,  21  Augmented  Hessian 
iterations  are  required  to  find  the  negative  eigenvalue  of  the  Hessian,  then  17 
cycles  are  required  to  find  the  stationary  point.  This  search  is  of  the  SCF 
type:  thus  the  calculated  transition  state  might  not  be  that  actually  found  as 
separation  into  ionic  products  would  require  a  quantum  chemical  treatment 
including  electron  correlation.  In  addition,  separation  into  ionic  products 
might  also  require  consideration  of  solvent  effects  (41). 

D.  STRUCTURE  AND  PROPERTIES  OF  CATIONS 


1.  Introduction 

Two  possible  mechanisms  have  been  proposed  for  the  formation  of  soot  in 
fuel  rich  flames.  One  mechanism  invokes  radical  reactions,  the  other 
ion/molecule  reactions  involving  small  ions.  The  latter  mechanism  has  been 
advanced  by  the  observation  that  ^3^3^  is  most  intense  ion  observed  in 

sooting  flames  (7,10)  and  motivates  this  theoretical  study  along  with 


32 


experimental  studies  of  ion/molecule  reactions  involving  small  hydrocarbon  ions 
using  Fourier  Transform  Ion  Cyclotron  Resonance  techniques  (Fq-  4?).  A  good 


characterization  of  possible  isomers  would  be  helpful  identifying  their 

presence  in  sooting  flames.  The  presence  of  such  species  are  experimentally 


ruscoplc  techniques  and  the  fech'^iques 


I-  - - -  :  -j 


being  used  cannot  distinguish  isomeric  structures  unambiguously. 


Various  mass  spectrometric  studies  suggest  that  possesses  reactive  and 

nonreactive  componenco.  The  most  stable  form  of  the  cyclopropenylium 

cation  I,  Figure  2,  is  relatively  unreactive  toward  acetylene  and  polyacetylenes 
as  well  as  toward  other  small  hydrocarbons  (42,^2)  The  propargylium  cation  II, 
is  considerably  more  reactive  toward  C2H2,  and  leads  to  the  formation  of 
and  ions,  both  with  reactive  and  unreactive  components.  The  propargylium 

cation  II  is  believed  to  lie  some  25  kcal/mol  (44)  higher  in  energy  than  the 
cyclopropenylium  ration  I,  but  this  species,  as  well  as  others  potentially 
important,  might  be  formed  under  the  reaction  conditions  observed  in  sooting 
flames.  Ve  thus  embark  in  this  study  to  uncover  possibly  reactive  species  of 

Carbocacions  such  as  ^3^3*  have  been  of  interest  for  some  time. 
Cyclopropenylium  ions  represent  the  first  members  of  the  2  +  4n  aromatic  series 
and,  as  such,  have  been  the  object  of  much  theoretical  discussion.  A  rather 
complete  study  of  the  relative  stability  of  many  structures  has  been 

reported  by  Radom,  Hariharan,  Pople  and  .5rhleyet  (45).  Although  the  earlier 
calculations  are  of  somewhat  lesser  accuracy  then  those  we  report  here,  in 
general  we  agree  with  their  ordering  of  stabilities  and  their  structures,  with 
one  exception  as  discussed  later. 


2.  Methods. 


The  initial  structures  of  this  study  were  obtained  through  INDO 
calculations  (46).  These  geometries  were  then  refined  using  the  ST0-3G,  then 
3-21G  and  finally  6-31G*  basis  sets  advocated  by  Pople  (Eq.  47).  The  use  of 
this  larger  basis  set  has  been  prompted  by  our  observation  that  stable  minima  at 
the  3-21G  level  can  become  "unstable"  transition  states  at  a  higher  level  of 
calculation.  Such  a  case  will  occur  here.  A  Broyden-Fletcher-Goldfarb-Shanno 
quasi-Newton  update  procedure  was  used  to  obtain  all  structures  (30,  48)  with  no 
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constrained  degrees  of  freedom.  The  ^  initio  calculations  were  performed  with 
the  Gaussian  82  program^  suitably  modified  with  this  searnh  algorithm.  All 
calculations  begin  with  structures  of  no  symmetry;  when  symmetry  was  obtained  it 
was  as  a  result  of  the  geometry  search.  The  geometries  believed  most  accurate 
are  those  obtained  from  the  6-31G*  calculation,  although  all  the  calculated 
structures  are  reported  for  comparisons.  Vibrational  frequencies  were 
calculated  using  both  the  3-21G  geometries  and  wavefunctions  and  the  6-31G* 
geometries  and  wavefunctions  (49),  and  they  are  unsealed.  When  comparing  with 
experiment  multiplication  by  a  factor  of  0.85  has  been  found  to  greatly  improve 
predictions.  At  optimized  6-31G*  structures  single  point  6-31G**  SCF 
calculations  were  performed  and  then  corrected  through  fourth  order  in  Mdller- 
Plesset  Rayleigh-Schrodinger  perturbation  theory  but  leaving  the  core 
"electrons"  uncorrelated  (50).  Our  best  calculations  on  the  relative  energy  of 
these  structures  are  thus  characterized  by  MP4-6-31G**//6-31G* . 

The  electronic  spectra  of  these  optimized  structures  were  estimated  using 
the  INDO-CI  spectroscopic  method  (51,  52). 

3.  Results 

the  relative  energies  obtain  for  five  structures  are  summarized  in  Table 
8.  Calculated  geometries  are  summarized  in  Figures  2-6  and  Tables  9  -  13.  As 
expected,  all  calculations  predict  the  cyciopropenylium  ion  (I)  lowest  in  energy 
followed  by  the  propargylium  cation  (II).  The  best  estimate  of  the  difference 
in  energy  between  those  two  structures  is  28.9  kcal/mol,  in  reasonable  agreement 
with  the  value  of  34.4  kcal/mol  calculated  by  Radom  et.  al.  (45).  Correction 
for  zero  point  vibrational  energy  reduces  this  value  to  27.7  kcal/mol,  to  be 
compared  with  an  experimental  estimate  of  25+4  kcal/mol  (44).  Structure  IV,  the 
prop-2-en-l-yl-3-ylidine  cation  is  calculated  next  higher  in  energy,  some  72.4 
kcal/mol  higher  than  the  cyciopropenylium  ion,  a  value  reduced  to  70.9  kcal/mol 
corrected  through  inclusion  of  0.85  times  the  zero-point  vibrational  energy. 

Structure  III,  the  1-propynyl  cation,  possesses  symmetry  and  a  pair  of 
degenerate  pi-like  orbitals  sharing  two  electrons.  Such  a  situation  gives  rise 
to  a  triplet  and  three  singlets.  The  lowest  energy  singlet  is  calculated  using 
complex  orbitals  and  is  estimated  to  lie  some  101  kcal/mol  higher  in  energy  than 
the  cyciopropenylium  cation  reference. 
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TABLE  8. 

RELATIVE 

STABILITIES 

OF 

ISOMERS  (A.U. )® 

I 

Cyclic  (CH)^ 

II 

H.,-C-C-CH 

III 

H^-C-C- 

S^ 

-C 

(T)^ 

IV 

H2C-CH-C 

S^  (T)^ 

V 

Cyclic 

H2CCCK 

ST0-3G 

-113.6203^ 

0.0564 

0.2900 

(0.1136) 

0.1110  (0.1680) 

0.2049 

[0.18311^ 

3-21G 

-114.3296^ 

0.0201 

0.2819 

(0.1117) 

0.0820 

0.1991 

[0.1705]^ 

6-31G* 

-115.0070*^ 

0.0559 

[0.1928]'^ 

(0.1413) 

0.1104 

0.2090® 

6-31G**^ 

-115.0133^ 

0.0564 

0.3160 

0.1111 

0.2236 

[0.1935]^ 

+MP2 

-115.3694^ 

0.0526 

0.3253 

0.1283 

0.2659 

[0.18731*^ 

+MP4 

-115.4081'^ 

0.0460 

0.2949 

(0. 16 

0.1153 

0.2261 

Radom  et. 

al.® 

0.0548 

[0.17531^ 

0.1109 

0.2039 

Exp .  ^ 

O 

fl 

o 

o 

o 

01 

a)  1  a.u.  =  627.5  kcal/mole.  S  =  Singlet,  (T)  =  Triplet 

b)  Calculated  total  energies. 

c)  Complex  calculation. 

d)  Estimated. 

e)  Reference  45. 

f)  Reference  44. 

g)  6-31G*//3-21G 

h)  6-31G**//6-31G* 
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TABLE  9.  THE  CALCULATED  GEOMETRIES  (A)  OF  THE  CYCLOPROPENYLIUM  ION,  I, 
FIGURE  2. 


COORDINATE 

INDO 

STO-3G 

3-21G 

6-31G* 

CH 

1.087 

1.095 

1.062 

1.072 

CC 

1.386 

1.377 

1.361 

1.349 

CCC 

60° 

60° 

60° 

60° 

TABLE  10. 

THE  CALCULATED  GEOMETRIES 

OF  THE  PROPARGYLIUM  ION  II, 

FIGURE  3. 

COORDINATE 

INDO 

ST0-3G 

3-21G 

6-31G* 

C1-C2 

1.407 

1.360 

1.334 

1.342 

C2-C3 

1.341 

1.214 

1.213 

1.214 

C3-H4 

1.079 

1.091 

1.061 

1.067 

C1-H5 

1.091 

1.109 

1.078 

1.078 

H5-C1-C2 

123.0° 

121.1° 

120.5° 

120.5° 

TABLE  11. 

THE  CALCULATED  GEOMETRIES  OF  THE 

Singlet 

1-PROPYNYL  CATION 

,  III,  FIGURE  4. 

Triplet‘s 

COORDINATE 

INDO® 

3-210^^ 

6-31g'^ 

INDO'^ 

3-21G 

6  31G* 

C1-C2 

1.382 

1.418 

1.434 

1.430 

1.435 

1.445 

C2-C3 

1.242 

1.357 

1.359 

1.272 

1.267 

1.257 

Cl-H 

1.110 

1.095 

1.091 

1.11 

1.082 

1.082 

H-C1-C2 

110.2° 

110.2° 

109.2 

108.9° 

112.3 

111.7 

a)  Two  determinant  ROHF  optimization. 

b)  Complex  orbital  optimization. 

c)  UHF  calculations. 

d)  A  slight  break  in  symmetry. 
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TABLE  12.  THE  CALCULATED  GEOMETRIES  OF  THE  PROP-2-EN-1-YL-3-YLIDINE  CATION,  IV, 
FIGURE  5. 


COORDINATE 

ST0-3G 

3-21G 

6-31G* 

C1-C2 

1.377 

1.372 

1.378 

C2-C3 

1.396 

1.347 

1.346 

H5-C1 

1.102 

1.074 

1.077 

H6-C1 

1.102 

1.076 

1.078 

H3-C2 

1.090 

1.081 

1.078 

H5-C1-H6 

116.9° 

117.3° 

117.0° 

H5-C1-C2 

121.2° 

120.3° 

120.9° 

C1-C2-H3 

123.6° 

122.9° 

125.9° 

C1-C2-C4 

115.8° 

118.0° 

108.0° 

TABLE  13. 

THE  CALCULATED 
FIGURE  6. 

GEOMETRIES 

OF  THE  CYCLOPRO-l-YL 

2-YLIDINE 

CATION  V, 

COORDINATE 

INDO 

ST0-3G 

3-21G 

6-31G* 

ST0-3G® 

C1-C2 

NOT 

1.478 

1.438 

NOT 

1.456 

C2-C3 

S 

1.700 

1.754 

S 

1.641 

C1-C3 

T 

1.507 

1.465 

T 

1.530 

H3-C2 

A 

1.112 

1.076 

A 

1.104 

H5-C1 

B 

1.102 

1.098 

B 

1.109 

C1-C2-C3 

L 

56.1° 

53.5° 

L 

58.8° 

C2-C3-C1 

E 

54.5° 

52.1° 

E 

54.6° 

C3-C1-C2 

69.4° 

74.4° 

66.6° 

H3-C2-C1 

123.9 

126.5° 

145.1 

a)  From  reference  45. 

The  optimized  triplet  structure  lies  at  considerably  lower  energy  than  the 
corresponding  singlets.  This  structure  is  predicted  to  lie  89  kcal/mol  above 
the  cyclopropenylium  ion  reference,  but  higher  in  energy  than  either  of  the 
singlet  structures  II  and  IV. 


Comparisons  between  states  of  different  multiplicities  such  as  we  have  made 
here  for  structure  III  are  notoriously  difficult  for  all  but  the  most  exacting 
ab  initio  theories,  as  these  calculations  tend  to  over-stabilize  triplet  state 
energies  relative  to  singlets.  Considering  the  partially  filled  degenerate 
orbitals  of  III, 

1L_  —  AL  JL  ±  J-JL  _L_L 

"^a  %  '^c  '^d  '^e  '*'f 
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the  use  of  complex  orbitals  approximates  the  —  ( ij/  -  'i',  )  component  of  the 

rj  a  b 

1  3^^ 

(the  other  component  is  —  ( The  state  is  estimated  by  a  UHF 

il  ^ 

calculation  on  ^  the  (the  other  components  are  4/r  and  —  (\J;  +  ^  ^  • 

e  t  ^  c  d 

Our  calculated  energy  differences  between  these  two  states  is  likely  too 

large.  The  INDO/S-CI  method,  however,  has  been  found  quite  reliable  in 

estimating  these  differences  at  fixed  geometries.  These  calculations, 

“  3 

performed  at  the  optimized  triplet  geometry,  suggest  that  the  state  lies 

1  1  ^ 

only  11  kcal/mol  lower  than  the  E  state,  and  that  the  A,  state  (l//2(ii/  +  )) 

2^  1  3  0 

lies  11  kcal/mol  higher  in  energy  that  the  ^E.  The  difference  between  the  ab- 
1  3 

initio  E  -  A^  energy  of  32  kcal/mol  with  that  of  11  kcal/mol  estimated  by  the 
INDO/S  calculations  corresponds  closely  to  the  observation  that  ab  initio 
calculations  of  this  type  tend  to  artificially  favor  triplets  by  about  20 
kcal/mol . 

The  cyclic  structure  V,  the  cycloprop-l-yl-2-ylidine  cation,  is  estimated  to 
lie  142.9  kcal/mol  above  the  reference  compound  I  at  the  3-21G  level.  At  the 
6-31G*  level,  structure  V  is  not  stable  and  opens  to  become  structure  IV.  This 
finding  is  in  contrast  to  previous  studies  (45).  With  the  inclusion  of 
polarization  functions  the  cyclic  form  V  opens  to  IV,  but  the  angle  that  has 
opened,  C1C2C4  of  Figure  5,  has  decreased  with  this  better  basis  from  118° 
(3-21G)  to  108°  (6-31G*),  Table  12. 

Structure  V  has  two  pi-like  electrons  (a").  A  near  lying  state  with  four 
pi-like  electrons  also  is  unstable,  leading  back  to  the  propargylium  cation,  II. 
Other  structures  have  been  examined  by  Radom  et  al.  (45)  but  those  are  all 
predicted  to  be  less  stable  than  these  we  report  here.  Although  these  latter 
calculations  are  reported  at  the  Hartree-Fock  level.  Table  8  shows  that  the  use 
of  a  slightly  better  basis  set  and  the  inclusion  of  a  fourth  order  perturbation 
theory  does  not  affect  the  differences  in  energies  by  more  than  about  107^. 
Note,  however,  that  differences  calculated  at  second  order  deviate  more  from  the 
fourth  order  results  than  do  the  Hartree-Fock  differences  themselves.  The 
inclusion  of  polarization  functions  in  the  6-31G**  calculations  preferentially 
stabilizes  the  cyclopropenyl  cation  I,  and  destabilizes  structure  V. 
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Table  9,  Figure  2.  summarize  the  calculated  geometry  of  cyclopropenylium 
cation.  As  usual,  the  better  basis  set,  yields  shorter  bond  length  then 
does  the  minimum  basis  set  ST0-3G.  Tire  INDO  results  are  in  good  accord. 


Table  Iti  and  Figure  1.  compare  rlie  calculations  for  the  propargylium  cation. 
The  INDO  and  ab  initio  results  suggest  different  bonding  schemes  for  this  ion. 

INDO 
H 

\ 

C  ^  C  =  C  -  H  '/s. 

/ 

H 


res pec  t i vely . 


Table  11  and  Figure  4,  summarize  the  calculated  structure  of  the  1-propynyl 
cation.  The  singlet  structure  is  best  described  as 


H 

\  * 

H  -  C  C  -  C: 

H 


with  the  C  -  C  single  bond  shortened  by  hyperconjugation, 

h" 

H  -  C  =  C  =  C: 

/ 

H 


The  triplet,  however,  is  best  described  as 


H 

\ 

H  -  C  -  C  s  C. 
H  ■ 


ab-ini t lo 
H  ~ 

C  -  C  5  C  -  H 
/ 

H 
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Table  12,  Figure  5,  gives  the  i^alrulateH  gpometiies  of  the  prop-2-en-l -yl -2 
yiidine  cation  IV.  There  is  a  considerable  shortening  of  the  C2  -  C4  bond  when 
employing  the  better  basis  set,  implying  a  more  important  role  to 


H 

C 

/ 

H 


than 


H 

\ 

C 

/ 

H 


H 


/ 

/ 

C 

x\ 


-  >  C: 


Table  13,  Figure  6,  presents  a  summary  of  our  findings  for  the  cyclopro-1- 

yiidine  cation,  V.  The  structure  of  the  figure  is  that  obtained  with  the  1  210 
basis  set.  As  mentioned,  this  strurtiiie  opens  to  the  prop-2-en-l-yl  3  yiidine 
ion  IV  with  the  better  basis. 


Experimental  estimates  for  the  hear  of  formation  of  the  cyclopropenylinm  ion 
range  from  256  -  25h  kcal/mol  with  an  estimated  error  of  +3  kcal/mol  (44,  53). 
Table  8  can  then  be  used  to  estimate  the  heats  of  formation  of  other  tautomers 
of  this  series. 

The  vibrational  frequencies  we  obtain  for  these  lowest  lying  conformers  of 
are  reported  in  Table  14.  They  are  unsealed  and  we  again  remark  that 
multiplication  of  these  frequencies  by  a  factor  of  about  0.85  often  gives  better 
results  when  compared  with  experiment.  Of  inteiest  is  the  observation  that 
lower  frequency  modes  are  most  affecred  by  the  lietter  basis,  and  are  general Iv 
decreased.  The  cyclopropenylinm  ion  T  has  no  low  lying  frequencies, 
distinguishing  it  from  the  other  less  st.ahle  structines  of  this  study.  Ue  note 
the  very  low  frequency  of  structure  V  at  the  3-210  level  leading  to  ring  opening 
at  the  6-310''--  level. 

The  calculated  UV-visfhle  spectra  of  ■'i'^  snecies  are  presenrr^H  in  Ta^lo  15. 
The  3-21G  optimized  geometries  were  assumed  for  INDO/S-CI  calculations  (51,  52). 
This  method  has  been  found  accurate  for  a  great  many  hydrocarbons,  and 
especially  so  for  positively  charged  species  where  diffuse  states  are  not 
important.  For  comparisons  results  of  calculations  by  Tokada  and  Ohno  for  the 
cyclopropenylinm  ion  I  (54)  and  by  F.yler,  Oddershede,  Sabin,  Diercksen,  and 
Gruner  for  propargylium  ion  II  Eq.  (55)  aie  included.  The  former  ralmlations 
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TABLE  14.  CALCULATED  VIBRATIONAL  FREQUENCIES  USING  3-2 IG  A)OPTIMIZED 

STRUCTURES  AND  UAVEFUNCTTONS  AND  B)6-31G*  OPTIMIZED  STRUCTURES 
AND  UAVEFUNCTIONS,  IN  CM~^. 


I  II 


mode 

(a) 

(b) 

I(IR)^ 

I(R)^ 

(a)  b) 

KIR; 

KR) 

a2" 

907.9 

842.5 

83 

00 

bl  351.7  280.3 

45 

00 

e' 

1041.8 

1029.2 

27 

6 

b2  396.4  364.2 

25 

03 

— 

— 

__ 

b2  842.0  741.6 

74 

01 

e" 

1136.2 

1099.6 

00 

00 

bl  1089.1  1023.9 

02 

01 

— 

— 

-- 

b2  1182.9  1131.5 

01 

00 

a2' 

1209.2 

1164.0 

00 

00 

al  1222.8  1216.6 

77 

57 

e' 

1373,6 

1418.3 

74 

04 

bl  1281.5  1265.1 

06 

00 

_ 

— 

— 

-- 

al  1613.0  1611.1 

03 

02 

al ' 

1702.3 

1795.2 

00 

53 

al  2216.5  2194.9 

816 

28 

e' 

3466,0 

3467.5 

86 

28 

al  3275.2  3331.1 

16 

101 

_ 

— 

— 

-- 

b2  3380.2  3443.3 

45 

58 

al ' 

E  (au) 
0 

3519,1 

0.0487 

3518.5  00 

0.0486 

III  Singlet 

91 

al  3535.8  3569.5 

0.0464  0.0460 

III  ' 

75 

Triplet 

47 

Mode 

(a; 

(b,d) 

I(iR) 

I(R) 

Mode  (a) 

KIR) 

KR) 

e 

519.1 

306.2 

5 

0 

e  480.9 

41 

791 

- - 

— 

- 

— 

— 

— 

al 

820,4 

885.7 

49 

0 

al  839.9 

74 

03 

e 

1266,4 

973.8 

36 

0 

e  1171.1 

10 

2096 

— 

— 

— 

- 

— 

— 

— 

al 

1560.0 

1499.2 

24 

0 

al  1542.9 

76 

103 

e 

1616,9 

1484.9 

— 

0 

e  1549.1 

106 

773 

— 

— 

49 

0 

— 

— 

— 

al 

2427.4 

1814.7 

5 

0 

al  1765.3 

17 

1414 

al 

3201.9 

3137.1 

125 

0 

al  3110.2 

171 

166 

e 

3296,4 

3223.9 

55 

2 

e  3187.5 

38 

07 

E®(au) 

o 

0.0488 

0.0440 

47 


TABLE  14.  CALCULATED  VIBRATIONAL  FREQUENCIES  USING  3-2 IG  A)OPTIMIZED 

STRUCTURES  AND  WAVEFUNCTIONS  AND  B)6-31G*  OPTIMIZED  STRUCTURES 
AND  WAVEFUNCTIONS,  IN  CM  (CONCLUDED). 


IV  V 


Mode 

a) 

b) 

I(IR) 

I(R) 

Mode 

a) 

a' 

337.0 

254.6 

09 

07 

a" 

157.2 

a" 

521.9 

463.7 

11 

02 

a' 

417.5 

a" 

1024.4 

996.1 

65 

00 

a 

630.9 

a' 

1091.3 

1013.4 

09 

02 

a' 

935.1 

a' 

1256.0 

1251.4 

13 

11 

a" 

1003.7 

a" 

1298.0 

1265.1 

10 

01 

a ' 

1051.3 

a' 

1358.4 

1383.3 

94 

29 

a' 

1194.2 

a' 

1649.3 

1648.4 

254 

15 

a' 

1311.2 

a' 

1699.3 

1685.7 

264 

04 

a' 

1333.0 

a' 

3280.8 

3333.2 

21 

52 

a' 

3087.6 

a' 

3303.3 

3372.8 

67 

63 

a" 

3114.2 

a' 

3409.6 

3444.9 

40 

40 

a' 

3342.5 

E®(au) 

0 

0.0461 

0.0458 

0.0400 

a)  3-21G  basis 

b)  6-31G*  basis 

c)  Not  a  stable  structure,  see  text 

d)  Numerical  force  constants 

e)  Zero  point,  vibrational  energy  in  Hartrees. 

f)  The  IR  and  Raman  intensities  in  km/mol  and  (Angstroms)^/amu,  respectively. 
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TABLE  15.  CALCULATED  ELECTRONIC 


-  I  - 

This  Work 

Ref . 

.  54*^ 

^A" 

^1 

56.1 

^A" 

69.0 

60.6 

Ir" 

71.0 

60.6 

71.0 

^E' 

64.4 

(0.166)^ 

^A” 

2 

71.0 

64.4 

(0.166) 

^E' 

84.7 

^A" 

62.2 

(0.003) 

84.7 

-  Ill 

- 

- 

III  - 

Singlet  (^E) 

Triplet  (^ 

1 

^2 

3.8 

^E 

14.4 

14.4 

(0.004) 

(0.004) 

18.4 

19.6 

19.6 

(0.008) 

(0.008) 

^E 

23.0 

23.0 

(0.000) 

(0.000) 

^E 

46.0 

46.0 

(0.000) 

(0.000) 

23.3 

^E 

39.5 

39.5 

(0.004) 

(0.004) 

49.0 

(0.000) 

58.9 

(1.620) 

45.3 

(0.405) 

^"2 

-3.8 

— 

^E 

3.8 

3.8 

— 

SPECTRA  OF  C^H^t  (1000 

CM 

-b. 

This  Work 

II  - 

Ref. 

.  55 

^-2 

17.9 

18.8 

41.0  (0. 

092) 

42.2 

(0.154) 

49.4 

60.1 

52.4 

66.2 

(0.000) 

59.6  (0. 

003) 

73.0 

(0.009) 

59.8  (0. 

807) 

73.8 

(0.787) 

^"2 

14.8 

— 

-  IV  - 

- 

V  - 

This  Work 

This  Work 

^A" 

31.7  (0.019) 

^A" 

24.7 

(0.000) 

'a' 

35.6  (0.001) 

^A" 

33.9 

(0.024) 

^A" 

35.6  (0.001) 

^A" 

43.8 

(0.030) 

^A' 

44.7  (0.461) 

^A" 

49.9 

(0.003) 

^A" 

56.1  (0.021) 

^A' 

57.3 

(0.114) 

^A' 

65.4  (0.000) 

'a' 

67.9 

(0.180) 

66.8  (0.097) 

^A" 

17.4 

— 

^A' 

20.3 

(a)  The  numbers  in  parentheses  are  oscillator  strengths. 

(b)  Estimated  from  Figure  6  of  reference  54. 
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are  an  ^  initio  SCF-CI  calculations,  but  employ  a  rather  small  basis  set  for 
this  purpose.  Our  results  differ  considerably  from  those  of  T.shodo  and  Dhno. 
The  results  reported  by  Eyler  et  al.  on  the  propargylium  cation  II  are  from  good 
basis  set  polarization-propagator  methodology  (56).  The  lower  two  excitation 
energies  they  obtain  are  in  remarkable  agreement  with  the  INDO/S-CI  results. 
Other  states  they  calculate  are  10,000  -  15,000  cm  ^  higher  in  energy  than  those 
we  obtain,  but  the  ordering  and  estimated  oscillator  strengths  are  again  in  good 
accord . 


Of  some  interest  in  the  formation  of  soot  is  the  interconversion  of 
propargylium  cation  to  cyclopropenylium  cation.  The  former  can  form  an 
encounter  complex  with  acetylene,  and  disproportionate  back  to  either  the  cyclic 
or  linear  The  cyclic  isomer  does  not  appear  to  react  with  acetylene, 
based  on  evidence  obtained  from  deuterium  exchange  reactions  (4).  There  is  no 
evidence  that  cyclopropenylium  cation  tautomerizes  to  propargylium  cation  in  the 
absence  of  an  encounter  complex  with  a  neutral  molecule.  In  keeping  with  this 
we  have  not  been  able  to  find  a  transition  state  geometry  linking  the  two 
structures,  and  from  constrained  optimization  calculations  estimate  a  transition 
energy  of  at  least  60  kcal/mol. 


4.  Discussion 

We  have  examined  several  stable  structures  of  ^3^3^  ^nd  report  on  the 
five  most  stable  that  ve  have  uncovered.  The  cyclic  form,  the  cyclopropenylium 
cation,  is  lowest  in  energy,  followed  by  the  linear  propargylium  ion,  [H2CCCH]’*^, 
calculated  to  lie  some  27.7  kcal/mol  higher  in  energy,  in  reasonable  agreement 
with  an  experimental  estimate  of  25t4  kcal/mol.  Other  structures  reported  are 
calculated  at  least  25  kcal/mol  higher  in  energy  than  the  propargylium  ion.  All 
of  the  structures  we  report,  however,  have  reasonably  high  vibrational 
frequencies,  suggesting  all  but  Structure  V  are  resistant  to  interconversion. 
Structure  V,  cycloprop-l-yl-2-ylidine  cation  [H2CCCHJ^,  we  predict  is  not 
stable,  either  with  two  or  four  pi-type  electrons.  Under  the  conditions  in 
which  we  are  interested,  sooting  flames,  some,  or  even  all,  structures  might  be 
observed.  In  addition,  the  lowest  energy  conformer  might  not  be  the  most 
abundant,  or  the  most  reactive  with  ^2^2  ^4^2  higher  molecular-mass 
precursors  of  soot.  Indeed,  as  mentioned  earlier,  the  most  stable  cyclic  form 


In  attempts  to  locate  various  species  present  in  flames,  or  in  models  of 
tliese  flames  simulated  in  ion  cyclotron  resonance  experiments  (42),  ve  have 
calculated  IR  frequencies  and  electronic  transition  energies.  Especially  of 
interest  are  those  species  with  low  lying  allowed  electronic  transitions  below 
about  40,000  cm  K  Such  states  might  be  excited  through  modern  laser  techniques 
and  observed  through  laser  induced  fluorescence,  or  through  excitation  of  those 
ions  using  afterglow  techniques.  Above  40,000  cm  interference  from  other 
species,  especially  aromatics,  might  be  expected. 

E.  REACTIONS  OF  ACETYLENE  WITH  ISOMERIC  FORMS  OF  THE  IONS. 


1.  Introduction 

The  reactions  of  C^H^  with  small  neutral  hydrocarbons  such  as  acetylene 
and  diacetylene  are  postulated  as  an  initial  step  in  soot  formation,  followed  by 
sequential  condensation  and  condensation-elimination  reactions,  forming 
successively  large  ions,  and  eventually  yielding  soot  particles  (7,10,14b), 
Studies  of  these  reactions  using  Fourier  Transform  Ion  Cyclotron  Resonance  (FT- 
ICR)  techniques  seem  to  indicate  that  *^3^3^  does  form  encounter  complexes  with 
acetylene,  but  under  the  conditions  of  the  experiment  no  stable  compounds  are 
formed  (42).  In  addition,  there  are  apparently  two  forms  of  present  in 
these  experiments,  only  one  of  which  is  reactive.  Reactions  with  diacetylene, 
however,  do  lead  to  the  formation  of  species  such  as  and 


With  an  eye  to  this  mechanism  of  soot  formation,  we  have  previously  examined 
the  structures  of  possible  C^H^*  isomers,  and  have  characterized  these  species 
according  to  their  relative  energies,  and  their  vibrational  and  electronic 
spectra  for  possible  identification  in  sooting  flames  (6).  The  reaction  between 
and  C2H2  should  yield  then,  could  yield  starting 

materials  back,  with  perhaps  isomerization  of  the  ion^. 


C3H3 


^2^2 


H 

C^H, 


+9 

The  most  stable  structure  of  ,  also  presented  is  a  preliminary  study  for 

possible  reactions  of  acetylene  with  the  two  lowest  energy  isomers  of  C3H3‘^, 
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cyclopropenyl ium  cation.  Figure  7,  and  propargylium  cation,  Figure  8.  Many 
reactions  might  occur  between  a  neutral  molecule  and  a  positive  inn,  and  this  is 
to  be  the  case.  Here,  however,  only  the  propargylium  ion  reacts  without 
barrier.  In  these  studies,  only  one  of  these  reactions  leads  to  a 
thermodynamically  stable  product,  relative  to  dissociation  back  into  acetylene 
and  cyclopropenylium  cation.  Each  of  the  reaction  products  uncovered  are  then 
characterized  according  to  their  vibrational  spectra  to  insure  they  represent 
stable  molecules,  and  for  possible  identifications  in  experiments  with  sooting 
flames  and  in  FT-ICR  .''xperiments  on  systems  designed  to  model  sooting  flames. 

2.  Theoretical  Methods 

Possible  reactions  using  were  examined,  the  INDO  method  (46)  searching 
the  potential  energy  surface  with  a  Broyden-Fletcher-Goldf arb-Shanno  quasi- 
Newton  update  procedure  (30,35).  The  uncovered  structures  are  then  refined 
using  the  Gaussian  86^  ab-initio  program,  first  with  the  STO  basis  sets,  then 
3-21G  split  valence  basis  sets,  and  then  with  6-31G*  split  valence  plus 
polarization  basis  sets  (47).  Vibrational  frequencies  are  obtained  using  the 
3-21G  basis  at  the  geometry  obtained  for  this  basis.  The  correlation  energy  is 
estimated  utilizing  a  6-31G**  basis  at  the  6-31G*  geometry  and  second  order 
Moller-Plesset  Rayleigh-Schrodinger  perturbation  theory  (MP2/6-31G**//6-31G*) . 

3.  Results  and  Discussion 

Previous  calculations  show  that  the  cyclic  tautomer  lies  some  27.7 
kcal/mol  lower  in  energy  than  the  propargylium  cation  ,  in  reasonable  agreement 
with  the  experimental  estimate  of  25  kcal/mol  (44).  All  other  calculated 
structures  lie  at  least  20  kcal/mol  higher  in  energy  than  the  propargyl  form. 
Ue  find  that  the  cyclic  tautomer  forms  an  encounter  complex  of  no  specific 
structure  as  shown  from  the  final  geometries  given  in  Figure  9-11.  These 
structures  are  obtained  from  initial  geometries  of  no  particular  symmetry  ot*"er 
than  the  plane  of  Figures  10  and  these  initial  geometries  are  all  rather  similar 
to  the  final  structures  obtained,  but,  of  course  with  C2H2  to  distance 
greater.  At  the  MP2/6-31G**  level  of  theory  these  complexes  are  stabilized  by 
about  6  kcal/mol.  No  attempt  has  been  made  to  correct  this  energy  for  basis  set 
superposi tion,  as  it  was  not  our  intention  to  study  these  weak  complexes.  No 
further  reaction  occurs  without  sizeable  barrier. 
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Figure  8.  The  Propargylium  Cation,  The  calculated  geometries  are  given  in 
Section  II. D. 
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Figure  9.  An  Encounter  Complex  Formed  in  the  Association  of  with  the 

Cyclic  Form  of  The  resulting  structure  reflects  the  initial 

geometry  used  for  these  calculations. 


ss 


Figure  10. 


An  Encounter  Complex  formed  in  the  Association  of  ^2^2 
Cyclic  Form  of  . 
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Figure  11.  An  Encounter  Complex  formed  in  the  Association  of  C2H2  with  the 
Cyclic  Form  of  C3H2‘^. 
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ono  mole  ol  ga.‘:,  wo  r.ui  o.';rim,at»>  tho  f  i  m>  onoi  gy  ol  t  h  i  ■;  loartioti  a.'; 

/v;  F.(r,iir)  ,  ().'»)F.('/ii))  rti.i/')H  1  ft.ovr 

wfioro  tho  Ffr.ilr)  i  ■;  at  t  lio  MP?/f)  ’,lf;'*  lovol  of  ihooty,  ami  T  1  t  lu- 

t  ompoi a  t  ui o  i u  Ko Ivin. 

With  tlu'  initi.al  gmimotiy  in  Figuio  Ma ,  ami  with  tho  arotylom'  'Might  ly 

'akowofi  .ukI  ili  lor  fly  ovoi  tho  ronii.tl  atom  of  tho  piopaigylium  ratifui,  wo 

followoil  tho  loai-  lion  iming  tho  optimization  ".rhoiiio  goiiotatod  by  tin-  Bi  oyilon 

Flotfhor  f'.oldf.nh  hh.inuo  method  (ht).  Initially  tho  arotylono  dtilt':  to  tho 
iiui.'it  po'iitivo  end  raihon  atom,  then  att.iiho';  to  t  h  i  •;  atom  to  foim  .a  thtof> 

momhoi  I'd  i  ing.  Figuio  ]  ’,h  ':how';  a  point  .tloiig  thi';  path,  and  1  \c  and  1  id  ato 

two  diffotoni  V  i  ow'i  ol  tho  piodmi  foimod.  Thin  ':poi  Ifm  i':  not  '.t.ahlo  ll■‘l.ativo 
to  fl  i  .';pi  opoi  t  i  on.i  I  i  on  into  tho  ryilit  i'.omoi  and  ,arf' t  y  1  fuio  ,  a';  mouliouod 

pr  ov  i  om:  1  y  . 


'.tt 


A  study  of  the  One  of  the  Reactions  between  Acetylene  and  the 
Propargylium  Cation,  a)  the  starting  geometry,  b)  an  intermediat 
structure,  c)  and  d)  two  views  of  the  final  product  of  this 
particular  encounter. 
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Figure  lA.  Tiie  Final  Structuit'  that  Result  1 1  om  flu-  Four  Initial  (lonriilion 
Shown  in  Figure  12.  Only  tlie  slrut  tUM'  ol  Figuie  1A<  is  stahlt' 
relative  to  the  ryr  1  opropeny  I  i  urn  ration  anrl  .ir  f  ty  I  etie  ,  sef>  text 
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TABLE  17.  CALCULATED  VIBRATIONAL  FREQUENCIES  OF  THE  REACTION  PRODUCTS 


USING  3-2IG  OPTIMIZED  STRUCTURES  AND  WAVEFUNCTIONS 


Structure 

Structure 

Fig.  14a 

Fig.  14b 

169.4 

73.2 

291.2 

293.5 

305.7 

364.8 

467.6 

519.9 

475.5 

665.3 

765.5 

699.8 

818.2 

860.1 

834.2 

1006.4 

908.3 

1010.5 

986.7 

1151.1 

1023.7 

1165.2 

1148.3 

1204.3 

1206.6 

1243.7 

1243.3 

1322.0 

1256.3 

1390.5 

1579.0 

1566.0 

1616.0 

1619.8 

1914.8 

1624.1 

2249.6 

1651.6 

3297.6 

3187.9 

3345.7 

3261.8 

3425.1 

3309.4 

3440.3 

3440.0 

3579.6 

3448.2 

Structure 

Structure 

Fig.  14c 

Fig.  14d 

150.7 

104.5 

358.9 

237.8 

425.1 

270.7 

774.2 

419.3 

793.2 

542.0 

815.9 

747.6 

954.1 

759.6 

1021.6 

888.0 

1070.6 

959.3 

1104.9 

977.6 

1190.2 

996.6 

1204.8 

1047.3 

1224.4 

1125.0 

1257.8 

1240.3 

1306.0 

1367.9 

1425.0 

1457.5 

1540.5 

1616.3 

1610.3 

1836.0 

1911.6 

2415.6 

3331.9 

3243.7 

3408.0 

3271.6 

3410.0 

3309.3 

3419.1 

3482.6 

3495.4 

3638.6 

(cm  )  . 


Structure 
Fig.  13d 

135.0 

231.3 

335.9 

517.1 
532.5 

560.7 
733.0 

869.9 

922.4 

985.9 

1043.1 

1089.2 

1115.4 

1125.3 

1159.7 

1427.1 

1534.3 

1947.1 

2106.3 
3285.0 

3368.4 

3409.3 

3500.9 

3586.2 
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F.  .STRUCTURES  AND  PROPERTIES  OF 


L .  Introduct Ion 

The  reaction  of  ^3^13^  with  small  neutral  hydrocarbons  such  as 
acetylene  and  dlacetylene  is  postulated  In  the  mechanism  for  soot  formation. 
This  is  followed  by  sequential  condensation  and  condensation-elimination 
reactions,  forming  successively  larger  ions,  and  eventually  yielding  soot 
particles. 

Previous  examinations  of  the  structures  of  possible  ^31*3^  isomers  had 

characterized  these  species  according  to  their  relative  energies,  and  their 

vibrational  .and  electronic  spectra  for  possible  identification  in  sooting 

flames^.  Ttie  reaction  between  and  C2H9  should  yield  C.^11.^^,  or 

4-  7  9 

cyclic  0.^11^  ’  • 

C^ll^  F  C2H2  y  ^5^^ 5  ^ 

C^fl3'''  F  U2 
C3H3+  +  C2M2 


(.4 


I'll  i  s  work  oxoinliicd  Liu-  musL  sl.ihli-  sL  riu- L  un-s  ot  Ikuh  iil  llu'  most 

stahh-  sLnu-mri-s  were  then  eh.i  r.ie  Ler  I  zed  leeo  rd  i  ii)>,  Lo  their  nv/v  i  s  i  ble 
siieetr.i  and  their  vibrational  sjjec-tra  Tor  possible  Idi-nt  i  1  leat  ion  in 
experiments  with  sooLinp,  1  lames  and  in  i-oupled  1  asi- r-t  T- 1 CK  experiments  on 
systems  desippied  to  motle  1  sooting;  tlame^:. 

2.  riu'oret  leal  Methods 

Possible  Isoitierlc  structures  o!  were  ex.iiriined  iisiii);  (In-  iNDO 

method  (4b,'31)  searchinp,  the  potential  enerj'y  siirlaee  wilh  a  Hroy<ien-Pleieher- 
Co  1  il  1  ar  b-Shanno  <iuas  1 -Newt  on  ii|Hlat<‘  procedure-  (dO,!')).  Pin-  uneove-red  struc¬ 
tures  wc-rt-  then  r(?t  I  lusl  usln^;  Llu?  (lausslan  Bb  ^  .ind  (laini-s?;  (')I!)  ab-Initlo  pro- 
)>rams,  lirst  with  the  STd-dd  basis  sets,  then  d-'ilC  spilt  val('nc(>  basis  si-ts, 
and  then  with  6-31(1*  split  valence  plus  polarl  z.it  Ion  basis  si-t ;;  Vibra¬ 

tional  f  renuenc  i  e-s  were  obtained  uslnp,  the  3-21(1  basis  at  the  p.eometry  ob- 
t.alned  for  this  basis.  The  correlation  energy  was  estlinatiid  utilizing  a  (>- 
31(1**  basis  at  the  6-31(1*  geometry  and  second  order  Mo  1  1  e r-P  1  essi? t  Itiy  1  e  l)>,h- 
Schrodlnger  perturbation  theory  ( MP2-6-31(l**//6-3 1  ll*  ) .  The  iiv/vlslbh-  spectra 
was  calculated  using  the  1NI)()/S  spectroscopic  mi-lhod  (31,32)  at  tlu-  optimized 
6-31(1*  geometries.  The  oscillator  strengths  of  the  transitions  w(-re  (jstlmaled 
ustu)'  the  dipole  length  oi)erator  and  only  one-center  L>-nns,  as  described 
elsewhere  (31). 

3.  Ke suits 

Twenty-one  possible  structures,  iincover*-ti  I  roin  the  INDO  calculations 
and  from  suggestions  rejjorted  In  tin-  literature,  wi-re  Initially  examined.  Pile 
nine  most  stable  structures  are  reported  In  Table  18,  and  their  optimized 
f'eometrles  calculated  at  the  6-31(;*  level  are  reported  In  Plgures  13  to  23. 
Tlu-  best  estimates  of  relative  energies  In  this  table  wi-rc-  obtained  t  rom  t  he 
MP2 -6-32(l**//6-3 1  (1*  ca  1  cu  la  t  I  ons .  It  Is  difficult  lo  comj)are  the  open  sht-ll 
structures  of  the  cyclic  species  IV,  V,  and  VI,  with  the  closed  ?!he  I  I  species, 
and  this  problem  Is  not  removed  through  the  addition  of  second  order  perturba¬ 
tion  theory  corrections.  Thes<-  sj)ecles  appear  to  be  biased  by  about  0.1)3 
llartrees,  or  .about  19  kc.al/rnol.  Phis  correction  h.is  been  Included  In  estl- 
m.ates  In  the  Tabli-  18,  .and  Is  consistent  with  re?iults  I  rom  tin-  INDO 


IN  IIAKTKKKS 


SKK 


rAlJLI'.  18.  IIIK  KNKKGIKS  OK  D1 K  KliRMNT  ISOMKRS  OF 

FICURKS  FOR  rilF,  CORRKSPONIJI  N(;  STKUCTURF.S . 


Me  tliod 

I 

1  1 

I  I  I 

IV 

V 

(Triplet) 

STO-  }(; 

-189.3968 

-189.5614 

-1 89.5724 

-189.6333 

-189.6011 

3-2  10 

-190. 81  AO 

-190.7364 

-190.7866 

-190.8433 

-189.8160 

0-310* 

1  l)e  t 

OVR 

-191.9183 

-191 .8652 

-191 .8828 

-191 .9334 

-191 .9024 
-191 .9340 

b-310**// 

6-310* 

-191.9278 

-191.8757 

-191 .8924 

-191 .9426 

-191 .9117 

+MP2 

-192.5468 

-192.5263 

-192.5124 

-192.5511 

-192.5322 

est Imated 
rel.  energy 
( kca 1 /mo  1 ) 

0.00 

12.9 

21.6 

16.4 

9.2 

Method 

VI 

VII 

VIII 

IX 

STO- 30 

-189.6006 

-189.5542 

-189.5542 

-189.5600 

3-2  IG 

-190.8160^ 

-190.8093 

-190.8984 

-190.8019 

6-310* 

1  Det. 

OVB 

-191.9021-’ 

-191.9309 

-191 .8801 

-191.8828 

-191 .8744 

6-310**// 

6-310* 

-191 .9105 

-191 .8891 

-1 91 .8919 

-191.8844 

KIP  2 

-192.5326 

-192.5134 

-192  5181 

-192.5049 

est  Imated 
rel.  energy 
( kca 1/mo  1 ) 

9.2 

21  .0 

18.0 

26.3 

This  .structure  Is  a  transition  state  leading  to  a  distortion  ol  the 
structure  given  In  Figure  20,  which  possesses  .in  energy  about  0.0002 
Hartrees  lower  tlian  the  r;tructure  shown. 


I 


Figure  16.  Structure  II 


Figure  18. 


Structure  IV. 
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H8 


Figure  19.  Structure  V. 
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Figure  21.  Structure  VII. 
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Figure  22.  Structure  VIII. 
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Figure  23.  Structure  IX. 
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spectroscopic  values  as  described  later. 


Heats  of  formation  of  ail  of  these  species  can  be  obtained  from 
consideration  of  the  following  reaction; 


C3H3''' 

HCCH 

C5H5+ 

AR 

experlment( kcal/mol) 

256±2 

54.19 

— 

— 

E(scf)  (3-21G,a.u.) 

-114.3296 

-76.3960 

-190.8140 

-0.0884 

E(scf)  (6-31G* ,a.u. ) 

-115.0070 

-76.8178 

-191.9183 

-0.0935 

E(scf+MP2)(6-31**//6-31G*) 

-115.3694 

-77.0795 

-192.5468 

-0.0979 

0.9*E(vib)(3-21G,a.u.) 

0.04382 

0.02930 

0.08131 

0.0082 

E(scf)+0.9*E(vib)(6-31**MP2) 

-115.3256 

-76.0502 

-191.4655 

-0.0897 

S  (cal/mol-K) 

54.4 

47.1 

69.7 

-31.8 

The  calculated  AE  of  the  reaction  is  -0.0897  Hartrees,  or  -56.3  kcal/mol. 
From  this  and  the  experimental  heats  of  formation  of  the  cyclopropenylium 
cation  and  acetylene,  we  estimate  aE  of  formation  of  the  vinyl- 
cyclopropenylium  cation,  structure  I,  to  be  256.7  ±  5  kcal/mol. 

Estimates  for  the  other  structures  of  this  study  can  then  be  made 
utilizing  the  energy  differences  in  Table  18.  We  remark  that  contributions 
from  the  zero  point  vibrational  energy  appear  to  make  some  difference  in  the 
energy  of  the  reaction,  reducing  it  by  5  kcal/mol.  Assuming  ideal  gas 
relations,  AH,  the  enthalpy  of  formation,  can  be  estimated  from 

AH  =  AE  -  1.98  X  10“^  T  (kcal/mol)  (^5a) 
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and  the  Gibb's  Free  energy  of  formation  from 


AG  =  AE  -  1.98  X  10"3  T  -  3.18  x  10~2  T  (45b) 

=  AE  -  0.0338  T 

The  entropy  of  each  species  is  estimated  using  the  Sackur-Tetrode  Equation 
(57),  the  calculated  vibrational  frequencies,  and  moments  of  inertia. 

The  vibrational  frequencies  calculated  for  these  low  lying  species  are 
reported  in  Table  19,  and  they  are  unsealed.  Multiplication  by  a  factor  of 
0.9  usually  gives  quite  satisfactory  agreement  with  observed  vibrational 
frequencies.  The  calculated  electronic  spectra  are  reported  in  Table  20. 

The  calculation  of  the  5  membered  cyclic  forms,  structures  IV,  V,  and  VI 
requires  further  comment.  For  the  cyclic  geometry,  as  with  the  triplet 

structure  V,  the  highest  occupied  orbitals  are  degenerate,  pi  like  and  of  Ej^g 
symmetry.  This  then  leads  to  six  possible  states, 

11-  —11  ±1  ±  L  L  L  ±  L 

A  B  C  D  E  F 

-1 

The  triplet  state,  of  A2g  type>  has  components  made  from  0(C),  0(D),  and 

(0(E)  +  0(F)),  where  0(A)  is  the  wavefunction  corresponding  to  diagram  A 
above.  The  state  has  components  (0(E)  -  0(F))  and  (0(A)  +  0(B)),  and  is 

O 

subject  to  a  first  order  Jahn-Teller  distortion.  The  ^Aj^g  state  lies 
considerably  higher  in  energy.  The  triplet  state  has  perfect  05^  symmetry  as 
Indicated  in  Figure  18.  Structures  V  and  VI  as  reported  in  Figures  19  and  20 
can  be  understood  as  arising  chiefly  from  0(A)  and  0(B),  and  these  structures 
are  estimated  from  a  two  determinant  generalized  valence  bond  calculation. 

4.  Discussion 

The  9  most  stable  C5H5'*’  structures  that  we  have  examined  are  given  in 
Table  18.  These  systems  have  been  examined  previously  (59-63),  and  we  note 
that  a  bent  cyclic  form  of  C5H5’^  is  not  stable  relative  to  the  planar  forms, 
in  contrast  to  the  suggestions  of  Dewar  and  Haddon  on  the  basis  of  HlNDO/3 
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TABLE  19.  CALCULATED  VIBRATIONAL  FREQUENCIES  USING  3-2 IG  OPTLMIZED  STRUCTURES 
AND  WAVEFUNCTIONS  (cm~^). 


I 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 

167.1 

247.4 

150.7 

505.3 

255.6 

-273.9 

167.8 

154.1 

155.9 

251.6 

— 

358.9 

— 

360.2 

389.9 

198.4 

195.4 

274.7 

412.1 

553.9 

425.1 

794.4 

469.0 

445.9 

365.4 

401.7 

380.9 

579.5 

858.6 

774.2 

940.9 

807.2 

811.5 

395.3 

453.5 

389.7 

768.1 

994.4 

793.2 

— 

895.4 

825.8 

554.1 

622.0 

574.0 

804.0 

— 

815.9 

977.4 

900.1 

931.9 

633.0 

644.1 

591.5 

994.1 

1002.4 

954.1 

— 

933.7 

833.1 

898.0 

672.2 

1089.6 

1003.5 

1021.6 

1006.1 

1007.7 

1005.3 

846.5 

1036.6 

1015.8 

1120.2 

— 

1070.6 

— 

1107.3 

1070.6 

1104.1 

1038.8 

1050.4 

1129.8 

1059.5 

1104.9 

1092.6 

1129.8 

1127.9 

1181.3 

1106.8 

1075.5 

1131.9 

1152.9 

1190.2 

— 

1158.6 

1159.0 

1212.8 

1138.9 

1106.8 

1235.7 

1155.3 

1204.8 

1290.8 

1176.2 

1178.2 

1238.1 

1211.1 

1112.5 

1245.9 

— 

1224.4 

1297.7 

1179.4 

1213.5 

1267.6 

1271.0 

1115.6 

1335.9 

1246.9 

1257.8 

— 

1267.2 

1268.6 

1430.5 

1359.2 

1173.1 

1466.4 

1259.6 

1306.0 

1549.5 

1355.0 

1356.3 

1469.9 

1434.1 

1490.6 

1479.9 

1295.4 

1425.0 

— 

1429.4 

1439.1 

1582  -  0 

1500  6 

1522.0 

1589.6 

1316.7 

1540.5 

1566.7 

1500.2 

1475.0 

1660.2 

1642.6 

1539.3 

1777.3 

1445.6 

1610.3 

1580.6 

1581.4 

1691.1 

1721.8 

1726.2 

2041.0 

1807.6 

— 

1911.6 

— 

1800.0 

1715.9 

2309.5 

2308.8 

2112.9 

3331.7 

3480.7 

3331.9 

3474.2 

3404.1 

3389.2 

3298.6 

3316.9 

3274.4 

3394.8 

3489.1 

3408.0 

— 

3404.7 

3430.7 

3332.6 

3324.4 

3279.1 

3425.5 

— 

3410.0 

3480.8 

3447.6 

3435.9 

3380.7 

3400.3 

3345.3 

3474.3 

3509.2 

3419.1 

— 

3465.7 

3465.9 

3395.9 

3425.7 

3360.6 

3512.8 

3585.9 

3495.4 

3688.8 

3474.0 

3472.2 

3431.6 

3594.0 

3360.9 
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TABLE  20.  CALCULATED  ELFXTRONIC  SPECTRA,  USING  INDO-CI  SPECTROSCOPIC  METHOD  FOR  6-31G*  OPTIMIZED  GEOMETRIES  (1000 
cm  THE  NUMBERS  IN  PARENTHESES  ARE  OSCILLATOR  STRENGTHS. 


p-  o  m 
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calculations  (62).  Similarly,  much  speculation  has  focused  on  the  stability 

of  the  pyramidal  form  of  structure  II,  postulated  to  lie  lowest  on  the 

surface  by  Stohrer  and  Hoffmann  some  time  ago  (59).  Although  this  is  most 

certainly  a  stable  structure,  as  indicated  by  the  reported  force  field  of 

Table  19,  it  is  far  from  our  most  stable  singlet  structure  I,  estimated  to  lie 
12.9  kcal/mol  lower  in  energy.  This  is  a  value  similar  to  that  reported  by 
Kohler  and  Lischka  (63),  but  an  estimate  of  theirs  using  the  CEPA  method 

predicts  the  square  based  pyramid  lower  than  the  cyclic  form  by  24  kcal/mol. 
Although  the  differential  correlation  energy  does  seem  to  favor  structure  II 
more  than  structure  IV,  this  is  only  by  about  19  kcal/mol.  Although  we  do  not 
believe  that  structure  II  is  lowest  in  energy,  it  is  a  well  defined  struci.-re, 
and  judging  from  our  calculated  force  constants,  stable  to  tautomerization. 
Experimentally  this  form  has  been  reported  to  exist  in  solution  (64,65). 
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SECTION  III 


EXPERIMENTAL  STUDIES 


A.  INTRODUCTION 


A  Nicolet  FTMS-1000  Fourier  Cransform  ion  cyclotron  resonance  (FTICR)  mass 
spectrometer  with  a  superconducting  magnet  of  fixed  magnetic  field  (3.0  tesla) 
was  used  for  the  work  reported  here.  Basic  principles  of  the  technique  and 
its  applications  in  ion/molecole  reaction  studies  have  been  reviewed  in 
several  recent  articles  (66-68) .  Reaction  pathways  were  delineated  using  the 
ejection  capabilities  (69)  of  FTICR  which  make  it  possible  to  eject  one  ion 
from  a  complicated  reaction  mixture  to  determine  its  contribution  to  the  mass 
spectrum  of  all  the  other  ions. 

For  studies  of  ion  spectroscopy,  a  Quantel  Nd:YAG-pumped  dye  laser  system 
was  used.  photodissociation  of  ions  trapped  in  the  FTICR  analyzer  cell  was 
followed  as  a  function  of  laser  wavelength.  Light  from  the  laser  was  brought 
through  sapphire  windows  Into  the  mass  spectrometer  vacuum  chamber,  and 
entered  and  exited  the  analyzer  cell  through  Cu  mesh-covered  holes  In  the  cell 
plates. 

B.  EXPERIMENTAL  METHODS 

Reaction  rate  coefficients  were  determined  by  monitoring  the  intensity  of 
the  C^Hy"^  ions  as  a  function  of  time  (typically  at  least  2  seconds)  after 
ejection  of  all  other  ions  from  the  analyzer  cell.  Neutral  gas  pressures  were 
measured  with  an  ionization  gauge.  Ionization  gauge  readings  were  then 
corrected  by  constructing  calibration  curves  of  ionization  gauge  versus 
capacitance  manometer  (MKS-Baratron)  readings  in  the  I  x  10“^  to  1  x  10”^  torr 
range.  A  correction  factor  was  required  because  the  ionization  gauge  and 
capacitance  manometer  were  located  at  different  points  on  the  vacuum  system, 
somewhat  removed  from  the  FTICR  analyzer  cell.  This  was  obtained  by 
determining  the  rate  coefficient  of  a  well-studied  reaction  (C2H^''‘  +  C2H^  -► 
products,  ^average  *  (1-0  i  0.3)  X  lO"^  cm^/s)!^  (70-74),  where  the  ionization 
gauge  pressure  readings  were  corrected  using  the  capacitance  manometer.  This 
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experimentally  determined  rate  coefficient  was  then  compared  with  the  average 
of  published  values  (70,74)  and  the  ratio  of  the  published  value  to  the 
measured,  which  was  3-3  ±  1.1^^  was  used  as  a  correction  factor.  This  factor 
was  used  in  calculting  the  absolute  rate  coefficients  given  in  this  report. 
All  calculations  of  rate  coefficients  and  95  percent  confidence  limits  were 
performed  using  a  menu-driven  FORTRAN  computer  program. 

Reactive  and  non-reactive  ions  were  produced  by  charge  transfer 

reactions  of  various  precursors  with  Xe"^,  formed  with  an  electron  beam  pulse 
of  5  ms  duration  at  an  ionizing  electron  energy  of  about  15  ev.  Propargyl 
iodide  was  used  as  a  precursor  for  C2H3'*'  in  studies  of  the  reaction  of  this 
ion  with  acetylene,  deuterated  acetylene,  and  diacetylene,  A  number  of 
different  precursors  (propargyl  iodide,  propargyl  bromide,  propargyl  chloride, 
propyne  and  allene)  were  used  to  investigate  the  reaction  mechanisms  leading 
to  CgH5^  ion  formation.  In  some  studies  C3H3'*’  ions  were  produced  directly 
from  the  above-mentioned  precursors  by  electron  impact  to  compare  the  effect 
of  ionization  technique  on  the  ratio  of  reactive  to  unreactive  isomers. 

Propargyl  iodide  was  prepared  from  propargyl  chloride  by  a  halogen  ex¬ 
change  reaction  (75).  The  details  of  purification  are  given  elsewhere  (76). 
propargyl  chloride,  propargyl  bromide,  allene,  propyne  and  acetylene  were  ob¬ 
tained  commercially  and  their  purity  was  checked  by  obtaining  wide  mass  range 
spectra.  In  the  case  of  acetylene,  some  production  of  protonated  acetone  was 
observed  at  long  delay  times  indicating  the  presence  of  a  small  amount  of 
acetone  as  an  Impurity.  Propargyl  bromide  was  distilled  before  use  in  order 
to  remove  toluene  which  was  present  as  a  stabilizing  agent.  Deuterated 
acetylene  was  prepared  from  D2O  and  CaC2.  Diacetylene  was  prepared  by  the 
method  of  dehydrochlorination  of  1 ,4-dichloro-3-butyne  in  aqueous  potassium 
hydroxide/dioxane  solution  (77). 

C3H3'*’  ions  were  produced  from  50  eV  electron  ioniziation  of  2,4-hexadiyne 
and  by  the  reaction  of  03113^  with  diacetylene.  *-3^3^  ions  used  to  produce 
^5^3^  were  formed  by  Xe"^  charge  transfer  reactions  with  propargyl  iodide  at  an 
ionizing  energy  of  15  eV.  After  30  ms  reaction  time  with  C4H2,  all  the  other 
tons  except  C3H3'*'  were  ejected  from  the  cell  to  follow  the  reactions  of  this 
ion  with  diacetylene  as  a  function  of  time.  ions  were  produced  by 


charge  transfer  reactions  of  various  precursors  (dicyclopentadiene,  cyclo- 
pentadiene,  norbornadiene,  l-penten-3-yne,  cycloheptatriene)  with  different 
charge  transfer  agents  (Xe"*",  N2'*’i  Ar"*")  formed  with  an  electron  beam  pulse  of  5 
ms  duration  at  ionizing  energies  of  13,  16.5,  and  20  eV,  respectively. 

Cyclopentadiene  was  prepared  by  cracking  dicyclopentadiene  (78)  and  was 
kept  in  dry  ice  when  not  used  to  prevent  the  dimerization  process.  All  the 
other  compounds  used  were  obtained  commercially  and  their  purity  was  checked 
by  obtaining  wide  mass  range  spectra.  All  the  samples  were  used  after 
multiple  freeze- pump- thaw  cycles. 

All  *"5^5'*'  r®^‘^tions  were  followed  at  a  cell  temperature  of  363 

K.  ^5^13"^  ^5%^  ions  used  in  rate  constant  determination  studies  were  all 

produced  by  charge  transfer  chemical  ionization  using  the  lowest  possible 
ionization  energy  in  order  to  minimize  the  internal  energy  imparted  to  the 
ions.  Also,  since  the  total  pressure  in  the  reaction  cell  was  almost  1  X  10~^ 
torr  and  the  charge  transfer  time  was  100  ms,  substantial  collisional 
relaxation  of  the  ions  took  place  before  kinetic  data  were  collected. 

ions  were  produced  from  toluene,  cycloheptatriene  and  norbornadiene 
by  dissociative  charge  transfer.  percentages  of  reactive/unreactive  isomers 
were  determined  when  the  three  neutrals  above  reacted  with  different  charge 
transfer  agents  (C^H2^>  Xe"*",  Kr"*"  and  Ar'*')  formed  with  an  electron  beam  pulse 
of  5ms  duration  at  ionizing  energies  of  11.5,  13,  15,  and  20eV, 
respectively.  To  do  this,  all  ions  except  CyHy'*'  were  ejected  from  the  cell 
and  enough  reaction  time  (typically  3s)  at  the  pressures  used  (typically  10“^- 
10“^  torr)  was  allowed  so  that  the  reactive  structures  could  be  lost  via 
lon/molecule  reactions,  thus  leaving  an  unreactive  population. 

As  reported  earlier  (79,80),  with  cycloheptatriene  unreactive  CyHy'*’  is 
formed  as  a  function  of  reaction  time  due  ion/molecule  Reactions  (46a)  and 
(46b) .  The  above  approach  thus  required  some  adjustments  in  the  experimental 
sequence  to  obtain  accurate  percentages  of  CyHy"^  populations. 

CyHy"^  CyH0  — — — >  CyHg"^  "h  CyHy  (46  a) 

CyHg  +  CyHg  — — >  CyHy  (unteactlve)  +  CyHg  (46b) 
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At  long  reaction  times  CyHg"^  was  ejected  continuously  during  the  reaction  time 
in  order  to  avoid  Reaction  46b. 

C.  REACTIONS  OF  C3H3‘''  WITH  ACETYLENE  AND  DIACETYLENE 

1.  Introduction 

As  discussed  in  Section  II.  of  the  many  possible  C3H3'^  isomers,  two 
are  most  important  in  assessing  the  role  of  C3H3'*’  in  ion/ molecule  reaction 
mechanisms.  The  first  is  the  cyclopropenylium  ^on,  recognized  as  the  most 
stable  isomer,  and  the  second  is  the  linear  propargylium  ion. 

The  importance  of  different  precursors  in  affecting  the  reactivity  of 
CjHs^  ions  was  reported  in  a  thermochetnical  study  by  Holmes  and  Lossing 
(81).  In  an  ICR  study  of  C3H3'^  reactions,  Ausloos  and  Lias  (82)  showed  that 
significant  fractions  of  the  linear  isomer  can  be  produced  by  charge  transfer 
reactions  of  small  ions  (Ar"*",  Xe"*",  CO^,  Ne'*',  etc.)  with  propargyl  chloride  and 
bromide.  Later  it  was  reported  (76)  that  even  higher  proportions  of  the 
propargylium  isomer  can  be  obtained  with  propargyl  iodide  either  by  electron 
impact  or  charge  exchange  using  Xe*^. 

A  study  by  Smyth  et  al.  (43)  demonstrated  that  the  cyclopropenylium  ion 
was  relatively  unreactive  toward  simple  hydrocarbon  fuels  whereas  the  linear, 
propargylium  ion  was  quite  reactive.  In  particular,  the  reaction  of 
propargylium  ions  with  acetylene  was  reported  to  produce  C^Hg"^  and  CgHg"^  ions 
with  an  overall  disappearance  rate  coefficient  of  1  x  10“^  cm^/s.  Since 

this  sequence  of  reactions  is  critical  in  the  initial  stages  of  the  postulated 
lon/molecule  mechanism  of  soot  formation,  it  was  studied  as  a  part  of  the  work 
carried  out  under  this  contract.  Results  substantially  different  from  those 
reported  earlier  were  found  and  the  study  was  expanded  to  investigate  possible 
production  of  C5H5^  ions  from  ionic  sources  other  than  C3H3'*’.  Reactions  of 
p’’opargylium  ions  with  C2D2  were  also  studied  in  order  to  delineate  further  a 
proposed  mechanism  for  the  CgHg'^/acetylene  reaction. 

Dlacetylene  is  another  important  flame  neutral  which  has  been  postulated 
to  react  with  C3H3'*'  ions  tn  an  ion/molerule  soot  formation  mechanism  (i4a). 
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Reaction  pathways  and  the  rate  coefficient  for  the  reaction  of  propargylium 
ions  with  diacetylene  near  room  temperature  were  thus  also  investigated.  In 
this  Section,  results  for  the  reactions  of  ions  with  acetylene, 

deuterated  acetylene,  and  diacetylene  are  reported  and  discussed  in  relation 
to  previous  work  involving  C2H3'*’  Reactions  (43,76,82)  and  the  proposed  (14a, b) 
ionic  route  to  soot  formation. 

2.  Results 

a.  0^113'^  reactions  with  acetylene. 

Despite  an  earlier  report  (43)  that  linear  C3H3'*’  is  quite  reactive 
with  acetylene,  only  very  low  intensities  of  C5H3'*"  and  C5H5'*’  produced  via  this 
reaction  could  be  found  in  this  work.  Experimental  conditions  of  the  earlier 
study  were  duplicated  as  closely  as  possible,  and  then  varied  substantially 
with  respect  to  relative  pressures  of  neutrals  (02^12/031131  ratio  =  1  to  8)  and 
overall  system  pressure  (from  5  X  10“^  to  3  x  10“^  torr) .  O3H3'''  ions  were 

formed  from  propargyl  chloride,  bromide,  and  iodide  by  both  charge  transfer 
using  Xe"^  and  electron  Impact.  In  order  to  determine  other  possible  sources 
of  observed  under  the  earlier  reaction  conditions,  binary  mixtures  of 

acetylene  and  each  one  of  the  C3H3'*’  precursors  reported  (43,76,82)  earlier 
were  used.  Intensities  of  and  were  first  measured  after  a  125  ms 

reaction  time.  Then  the  parent  ion,  C2H2^,  and  €3113'^  were  each  ejected 
separately  during  the  125  ms  reaction  period  to  assess  their  contribution  to 
CgHg'*’  and  CgHg"^  formation. 

For  each  different  neutral  precursor,  the  sources  and  amounts  of  C^Hg'*’ 
ions  produced  were  found  to  be  different.  Propyne  and  allene  were  similar  in 
producing  large  amounts  of  C^Hg"^  and  no  ion.  However,  the  C^Hg"^  +  C2H2 

reaction  was  not  responsible  for  CgHg"*"  formation.  The  main  reactions  leading 
to  CjHg''"  were 

C2H2‘'’  +  C3H4  +  CgHg'*’  +  H  (^7) 


and 
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(48) 


C3H4'*’  +  ^2^2  ^ 

in  both  cases.  On  the  other  hand,  when  propargyl  iodide,  propargyl  bromide 
and  propargyl  chloride  were  used  as  precursors,  relatively  smaller  amounts  of 
C5H5'*’  ions  were  observed,  but  ions  were  also  formed.  In  the  propargyl 

chloride  case 

C3H3CI  "h  ^2^2  ^  (49) 


and 


^2^2^  +  C3H3CI  >  C5H5‘''  +  Cl  (50) 

were  the  major  reactions  leading  to  CjHj'*’  formation.  For  propargyl  bromide 
the 


C3H3Br'''  +  C2H2  -*•  C5H5'^  +  Br  (51) 

reaction  was  the  only  source  of  C5H5'*’  ions  observed.  Any  contribution  to 
CgH5"'‘  formation  from  linear  C3H3'^  was  less  than  the  experimental  uncertain¬ 
ty.  Finally,  very  little  (almost  negligible)  amounts  of  ions  were 

observed  when  propargyl  iodide  was  used  as  a  precursor.  The  reactions 

C2H2‘^  +  C3H3I  +  €585'^  +  I  (52) 


and 


C3H3''‘  +  C2H2  C5H5''‘  (53) 

were  the  major  contributors  in  this  case.  An  upper  limit  for  the  rate  con¬ 
stant  for  Reaction  (53)  was  estimated  as  5  X  10”^^  cm^/s  by  assuming  that  the 
very  small  05115'^  signal  observed  resulted  from  this  reaction,  and  using  the 
expression  [C3H3''‘](t)  -  [C3H3‘^]o  -  [C5H5‘''](t)  -  [C3H3‘'']  (0)e“'^^*^ ,  where  n  is 
the  C2H2  number  density.  Overall  results  for  the  production  of  C3H5'’’  and 
C5H3'*’  ions  with  different  neutrals  which  have  been  reported  (43,76,82)  as  pre¬ 
cursors  of  C3H3'*'  are  summarized  in  Table  21. 
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TABLE  21. 

PRODUCTION  OF  C3H5 
AND  ACETYLENE.^ 

AND  C3H3''‘ 

IONS  IN  MIXTURES 

OF  VARIOUS 

NEUTRAuS 

Neutral 

Ionic  sources^ 
of  C3H5''’  after 
Xe"*"  charge 
transfer 

ionization  of 

a  mixture  of 

the  neutral 

and  acetylene 

Ratio  of 

prod. 

relative  to 

that  in 

allene  case 

Ionic  sources*^ 
of  C5H3''‘ 

c 

Intensity 

of 

vs.  C5H3'*’ 

Percent 

reactive 

C3H3'^ 

ions 

Allene 

02^2*^  (^0%) 

l.O 

- 

- 

<5 

C2H^'^(60%) 

propyne 

C2H2‘^(40%) 

0.75 

- 

- 

30 

C3H4'^(60%) 

propargyl 

C2H2‘''(40-50%) 

0.25 

20%) 

3.0 

15 

chloride 

C3H3C1‘'’(50-60Z) 

C3H3'^(20%) 

C2H2^(60%) 

Propargyl 

C3K3Br'''(90-100Z) 

0.08 

C2H2'^(70%) 

2.0 

85 

bromide 

[C3H3''’  + 

C3H3Br'^](30%) 

propargyl 

CzHz'^C^OX) 

<0.02 

C2H2‘^(50%) 

1.7 

90 

iodide 

C3H3'^(60%) 

C3H3'^(50%) 

®  All  ions  were  produced  by  chemical  ionization  charge  transfer  from  Xe"^. 

^  Percentages  show  the  relative  contributions  to  and  production  as 

determined  by  double  resonance  experiments  and  have  an  estimated  uncertainty 
of  *10Z. 

^  Neutral  reactants  all  had  the  same  pressure  (7  X  10”^  torr)  as  measured  by 
the  ionization  gauge.  Xenon  and  acetylene  pressures  were  5.6  X  10“^  and  1.3 
X  10~^  torr,  respectively. 


87 


Becduse  propargyl  iodide  produced  the  highest  reactive/unreactlve  ratio  of 
C3H3'*’  ions  both  in  earlier  (76,81)  and  the  present  work  (see  Table  22),  it  was 


TABLE  22.  PERCENTAGES^  OF  REACTIVE  C3H3'''  FOUND  FROM  VARIOUS  PRECURSORS  BY 
VARIOUS  IONIZATION  TECHNIQUES  (MONITORED  BY  OBSERVING  REACTION  WITH 
THE  PRECURSOR  NEUTRAL). 


Ionizing 

Precursor 

technique 

Propargyl 

Propargyl 

Propargyl 

iodide 

bromide 

chloride 

Electron  impact 

90% 

40% 

10% 

(15  eV) 

Chemical  ionization 

charge  transfer 

with  Xe"*" 

90% 

85% 

15% 

^  Estimated  error  is  i5%. 


used  as  a  precursor  for  C3H3'*'  Ions  in  these  reaction  kinetics  studies.  Since 
the  precursor  neutral  molecule  was  always  present  in  the  FTICR  analyzer  cell, 
it  was  a  competitor  with  the  reactant  neutral  of  interest  in  ion/  molecule 
reactions  involving  C3H3'*'.  To  determine  the  rate  coefficient  of  reactions  of 
^3^3^  with  the  neutral  reactant,  it  was  necessary  to  monitor  the  reactions  of 
this  ion  with  C3H3I  and  subtract  the  rate  coefficient  for  this  reaction  from 
the  total  rate  coefficient  observed  in  the  presence  of  both  the  precursor 
neutral  and  the  reactant  of  Interest.  Reactions  of  C3H3'*’  with  propargyl 
iodide  were  monitored  as  a  function  of  time  following  charge  transfer  chemical 
ionization  of  C3H3I  by  Xe"*"  and  ejection  of  all  ions  but  C3H3^  from  the 
analyzer  cell.  Results  obtained  were  identical  to  C3H3''"  reaction  channels 
with  propargyl  iodide  which  have  been  reported  elsewhere  (76). 

b.  Isomerization  of  linear  C3H3''’. 

In  addition  to  the  absence  of  any  significant  C3H3^  and 
formed  by  reaction  of  linear  C3H3^  with  C2H2,  it  was  also  observed  Chat  02^2 
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led  to  the  isomerization  of  linear  ions  to  their  cyclic  form,  thus 

rendering  them  unreactive  toward  their  parent  neutral  as  well  as 

toward  C2H2-  This  isomerization  was  followed  as  a  function  of  G2H2  pressure 
and  a  direct  pressure  dependence  was  found,  as  shown  in  Figure  24. 

c.  ^3^3"^  reactions  with  C2D'> . 

To  achieve  a  better  understanding  of  the  isomerization  of  linear 
C2D2  instead  of  C2H2  was  used  as  the  neutral  reactant.  The  following 
isotope  exchange  reactions  were  observed. 


C3H3‘'‘  +  C2D2  +  C3H2D‘^  +  C2DH 

(54) 

C3H3  +  C2D2  C3HD2  "h  C2H2 

(55) 

C3H2D''‘  +  C2D2  ->•  C3HD2‘'‘  +  C2DH 

(56) 

031120'^  +  C2D2  0303^  +  C2H2 

(57) 

C3HD2‘*'  +  C2D2  -*■  0303"^  +  C2DH 

(58) 

Using  FTICR  ejection  capabilities,  it  was  found  that  reactions  (55)  and  (56) 
contribute  equally  to  the  production  of  C2HD2'*’  while  reaction  (58)  produces 
more  of  0303"*“  (80  percent)  than  reaction  (57)  (20  percent). 

Ion  intensity  versus  time  curves  for  the  C3H3'''/C2D2  reaction  are  shown  in 
Figure  25.  The  overall  rate  coefficient  for  the  disappearance  of  C3H3'*'  was 
calculated  by  subtracting  the  observed  rate  coefficient  for  the  reaction  with 
propargyl  iodide  from  the  total  observed  rate  coefficient  in  the  presence  of 
C2D2*  This  observed  rate  coefficient  was  then  converted  to  the  true  rate  co¬ 
efficient  using  the  corrected  pressure  of  C2D2*  A  value  of  (4.5  i  1.9)  x 
10“^®  cm^/s  was  found  at  a  cell  temperature  of  373  K  for  the  disappearance  of 
^3%^  (reactions  (54)  and  (55)).  In  Figure  26,  ion  intensity  versus  time 
curves  of  C3H3''’  are  compared  for  reactions  with  and  without  C2D2* 

d.  ^3113'*’  reactions  with  dlacetylene 

After  ejection  of  all  ions  except  C3H3'*"  following  charge  transfer 
chemical  ionization  by  Xe'*’  of  a  mixture  of  dlacetylene  and  propargyl  iodide, 
the  ion/molecule  reactions  as  a  function  of  time  were  monitored.  Consecutive 
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Normalized  peak  area  of  C3H3 


Figure  24.  Isomerization  of  Linear  03113'^  Ions  at  Different  Pressures  of 
C2H2.  03113^  ions  were  produced  by  charge  transfer  reactions  with 

Xe"*".  Pr.zi.T  ■  1.1  X  10~^  torr;  was  adjusted  to  maintain  a 


PC3H3r 

constant  total  pressure  of  2.6  X  10 


-6 


torr  as  measured  on  the 


ionization  gauge, 
corrected. ) 


(All  pressures  are  capacitance-manometer 
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Figure  25.  Isotope  Exchange  Reactions  of  €3113  with  C2D2.  Disappearance  of 
C3H3”*'  ion  Includes  reactions  with  propargyl  iodide.  Note  that  the 
sum  of  all  isotopic  forms  of  031^3'^  remaining  at  the  end  of  the 
reaction  with  C2D2  approximately  equals  the  total  unreactive  C3H3'*’ 
when  C2H2  is  used  as  a  neutral  reactant  at  the  same  pressure  (see 

Figure  24).  PC3H3I  *  ^  torr;  PC2D2  *  ^  torr; 

Py-  *  6.2  X  10~^  torr. 


Peak  area/arb.  units 


Figure 


26.  C3H3 


Ion  Decay  Curves  for  Reaction  with  C3H3I  and  C2D2< 
(Pressures  are  the  same  as  given  for  Figure  25. ) 


92 


C2  and  addition  reactions  were  observed: 


C3H3“^  + 

C4H2 

-► 

C7H5+ 

(59) 

C3H3 

C4H2 

> 

C5H3''’  +  ^2^2 

(60) 

C5H3 

C4H2 

C9H5+ 

(61) 

C3H3  + 

C^H2 

> 

C^H3  ^2^2 

(62) 

C4H2 

C11H5+ 

(63) 

C4H2 

-► 

^11^5  ^2^2 

(64) 

09^5“*’  + 

C4H2 

-► 

C13H7+ 

(65) 

Some  of  these  product  ions  were  seen  to  react  further  with  propargyl 
iodide  by  displacement  of  atomic  iodine: 


^  (66) 

C7H2‘''  +  C3H3I  +  +  I  (67) 

+  C3H3I  ^ 

CgHg"^  +  C3H3I  >  +  I  (69) 

CgHj'*'  +  C^H3l  -►  0^2%''’  ^ 


Ion  intensity  versus  time  curves  for  the  C3H3'^/C^H2  ^^eaction  are  shown  in 
Figure  27.  The  rate  coefficient  for  the  disappearance  of  C3H3''’  (Reactions 
(59)  and  (60),  Figure  28)  was  calculated  as  described  earlier,  and  a  value  of 
k  -  (1.4  i  0.7)  X  10“^  cm^/s  was  found. 

propargyl  iodide,  bromide,  and  chloride  were  all  used  as  precursors  of 
C3H3'*'  in  studying  its  reactions  with  dlacetylene.  For  each  precursor,  both 
electron  Impact  and  charge  transfer  chemical  ionization  techniques  were 
used.  The  percentages  of  reactive  isomer  in  the  reaction  with  diacetylene  are 
shown  in  Table  23*  When  these  percentages  of  reactive  isomer  were  compared  to 
those  in  the  absence  of  C4H2  Table  22),  it  was  clear  that  some  isomer¬ 
ization  of  the  reactive  linear  0383"^  ion,  as  well  as  Reactions  (59)  and  (60) , 
had  taken  place  (see  also  Figure  5).  This  isomerization  was  followed  as  a 
function  of  0^82  pi^^ssure  and  a  direct  pressure  dependence  was  found,  as  can 
be  seen  in  Table  24. 
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Normalized  peak  area  of  C3H 


+  :  C3H3  with  C3H3I 
□  C3H3  with  C3H3  I  and  C^H 


Figure  27.  Reactions  of  C3H3''’  with  C4H2.  Disappearance  of  C3H3'^  and  product 
ions  include  reactions  with  propargyl  iodide.  Pc^H^I  ”  ^ 

_ _  _  _  o  V  1  n—7  _  t  •)  V  1  fnrr.  ^’All 


torr;  p^ 


4.8  X  10"'  torr;  pj^^  -  6.2  X  10 


pressures  are  capacitance-manometer  corrected. ) 


PeoK  area /arb.  units 


Figure  28.  C3H3'*'  Decay  Curves  for  the  Reactions  with  C3H3I  and  C4H2. 
(Pressures  are  the  same  as  given  for  Figure  27.) 
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TABLE  23.  PERCENTAGES^  OF  REACTIVE  C^H^'*’  OBSERVED  IN  THE  REACTION  WITH  C^H2-*^ 


Ionization 

Precursor 

technique 

propargyl 

iodide 

Propargyl 

bromide 

Propargyl 

chloride 

Electron  impact 

(15  ev) 

75% 

30% 

5% 

Charge  transfer 
by  Xe'*’ 

75% 

65% 

5% 

^  Estimated  error  is 
^  Pn  u  -  4.8  X  10"^ 

±5%. 

torr. 

TABLE  24.  CHANGES  IN,  C3H3'*'  REACTIVITY^ 

D I ACETYLENE. ° 

AT  DIFFERENT 

PRESSURES  OF 

pressure  of 

C4H2/10”^  torr 

%  of 

unreactive  C3H3‘'* 

0.8 

16 

1.6 

17 

4.8 

25 

7.2 

32 

8.0 

35 

9.6 

40 

®  I-C2H3"*'  ions  were  produced  from  propargyl  iodide  by  chemical  ionization 

charge  transfer  with  Xe"*",  (Pc3H3i  ■  i’i  ^  ^  torr;  pj^g  was  adjusted  to 

maintain  a  constant  total  pressure  of  2.6  X  10  ^  torr  as  read  on  the 
ionization  gauge) . 

^  All  pressures  are  capacitance-manometer  corrected. 
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3.  Discussion 


Although  the  results  of  this  work  are  not  in  agreement  with  the 
earlier  report  (43)  of  C3H3'*‘/C2H2  reactivity,  the  discrepancy  is  probably  due 
to  limitations  of  the  older  pulsed  ICR  (43)  instrumentation  for  studying 
ion/molecule  reaction  pathways  in  complicated  systems  when  compared  to  newer 
FTICR  capabilities.  Facile  ejection  of  all  ions,  except  the  one  whose 
ion/molecule  reactions  are  being  investigated,  offers  a  very  clean  monitoring 
opportunity  for  product-parent  relationships  even  in  complicated  consecutive 
and  competitive  reaction  systems.  Various  alternative  pathways  for  the 
production  of  and  which  have  been  described  above  probably 

contributed  significantly  to  the  intensities  of  these  ions  seen  in  the  earlier 
work.  Additional  support  for  the  low  reactivity  of  C2H3'*’  with  C2H2  is  found 
in  a  recent  report  (83)  in  which  the  rate  coefficient  for  this  reaction  is 
reported  to  be  less  than  0.01  X  10~^  cm^/s,  although  the  isomeric  form  of 
C3H3''’  is  not  given.  It  is  also  possible  that  the  highest  pressures  used  in 
this  work  did  not  reach  those  of  the  earlier  study  (43) .  Higher  pressures 
might  lead  to  third-body  stabilization  of  C5H5''’  collision  complexes. 

The  most  likely  mechanism  of  the  observed  isomerization  of  C3H3‘''  ions  by 
collisions  with  acetylene  is  a  "reactive”  rather  than  a  "nonreactive"  one. 
That  is,  it  results  from  an  intimate  encounter  of  the  ion  and  neutral  in  the 
C5H5'''  collision  complex.  This  hypothesis  is  confirmed  by  the  fact  that  deu- 
terated  forms  of  C3H3'*’  were  produced  when  C2D2  was  the  neutral  reactant  (see 
Figure  26).  In  most  cases  the  C5H5''’  collision  complex  dissociates  to  give  the 
cyclic,  unreactive  C3H3'*’  isomer.  Instead  of  the  linear  form.  The  possibility 
of  non-reactive  colllsional  isomerization  of  linear  C3H3'*’  to  the  cyclic  isomer 
has  been  ruled  out  because  experiments  at  elevated  pressures  of  xenon  (to  ca. 
1  X  10“^  torr)  showed  no  interconversion.  Similar  interconversion  of 
ions  from  a  linear  to  cyclic  form  has  also  been  reported  (84)  in  the  reaction 
with  C2H2  and  has  also  been  shown  to  take  place  via  complex  formation  by  using 
isotopically  labeled  C2R2'  To  confirm  the  hypothesis  that  higher  energy,  re¬ 
active,  (linear)  C3H3'^  ions  Interconvert  to  more  stable,  unreactive  ones, 
cyclic  C3H3''’  ions  were  also  reacted  with  C2D2  and  no  isotope  exchange 
reactions  were  observed. 
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plots  of  C3H3''’  ion  intensity  versus  time  for  reaction  with  C^H2  (Figure 
28)  indicate  a  10-12  percent  increase  in  the  intensity  of  the  unreactive 
isomer  relative  to  the  reaction  when  only  the  parent  precursor  is  present. 
Isomerization  of  reactive  C3H3''*  was  also  seen  when  different  precursors  were 
used  (compare  Tables  22  and  23).  A  similar  mechanism  involving  complex 
formation  may  be  responsible  for  this  isomerization  as  well,  although  it  was 
not  investigated  in  any  detail. 

The  extensive  ion/molecule  condensation  reaction  sequences  (Reactions  (59- 
65))  observed  when  03113'*’  reacts  with  C4H2  suggest  that  this  aspect  of  the 
proposed  (l4a,b)  ionic  path  to  soot  formation  is  quite  credible.  The  rate  co¬ 
efficient  determined  in  this  work,  for  C3H3"*’  disappearance  (1.4  ±  0.7  x  10”^ 
cm^/s)  is  in  good  agreement  with  that  determined  earlier  (43)  (1.0  ±  0.5  x 
10“^  cm^/s).  Some  of  the  product  ions  formed  (e.g.  C5H3"'’  and  Cyll^'*’)  have  been 
seen  to  be  abundant  in  both  nonsooting  and  sooting  flames  (10)  .  On  the  other 
hand,  the  observation  of  €3113'*’  isomerization  and  not  condensation  with  acety¬ 
lene  suggests  that  the  proposed  sequential  acetylene  addition  reactions  with 
^3*^3^  in  the  ionic  soot  formation  mechanism  be  reconsidered.  Other  channels 
such  as  direct  reaction  of  neutral  aromatics  with  03113'*’  (9,76)  may  be  as 
Important  in  the  formation  of  small  polycyclic  ions. 

These  experiments  and  those  reported  earlier  (43)  were  carried  out  at 
relatively  low  pressures  (<  5  X  lO'^  torr).  It  is  quite  possible  that  third 
body  collisions  in  atmospheric  pressure  flames  can  stabilize  a  fraction  of  the 
CjHj"*’  reaction  complexes  before  the  "reactive"  deactivation  observed  in  this 
work  can  take  place.  However,  direct  bimolecular  reaction  of  linear  C3H3''’ 
with  acetylene  to  form  C5H3'*'  and  C5H5'*’  does  not  appear  to  occur  readily. 

D.  REACTIONS  OF  C3H5'*'  AND  C5H3'*'  IONS  WITH  ACETYLENE  AND  DIACETYLENE 

1.  Introduction 

As  noted  in  Section  IIIC.,  study  of  the  reactions  of  C3H3'*'  with  acety¬ 
lene  and  dlacetylene  did  not  reveal  facile  formation  of  C5H5^  or  C3H3^  by  re¬ 
action  of  this  ion  with  acetylene,  although  C3H3'''  did  react  readily  with  di¬ 
acetylene,  yielding  C5H3'*'  and  C7H3'’’.  All  three  of  the  postulated  products  of 
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C3H3'*'  reactions  witn  acetylene  and  diacetylene  C3H3'*',  have  been 
identified  by  mass  spectrometry  (9,10)  in  flames,  but  it  remains  to  be  deter¬ 
mined  which,  if  any,  of  them  may  be  important  in  soot  formation  mechanisms. 
Thus,  as  a  part  of  this  project,  investigations  were  extended  to  include  study 
of  the  reactions  of  C3H3''’  and  ions  with  certain  flame  neutrals. 

Two  previous  experimental  studies  on  the  structures  of  C5H3'*’  involved 
determination  of  the  heats  of  formation  of  the  C5H3''’  ion  from  different  pre¬ 
cursors  (85,86).  The  results  of  these  studies  suggested  the  presence  of  two 
different  linear  structures  for  C3H3''’»  each  produced  from  a  different  pre¬ 
cursor  (85)* 

Many  possible  structures  exist  for  the  C3H5'*’  1-on,  and  despite  a  number  of 
theoretical  and  experimental  studies  involving  it,  few  definitive  results 
exist  regarding  the  specific  relative  energies  of  various  isomeric  forms 
(87,88).  Experimental  studies  concentrated  on  determining  the  appearance  po¬ 
tential  and  heats  of  formation  of  ions  from  different  sources  by  mass 
spectrometric  methods  (89-95).  These  resulted  in  heat  of  formation  values 
ranging  from  239  to  309  kcal/mole,  depending  on  the  source  and  technique  of 
C5H5'^  production.  Early  ICR  experiments  in  these  laboratories  were  carried 
out  to  identify  structures  of  C5H5’*’  according  to  their  reactivity  with  dif¬ 
ferent  neutrals  (87,88,96).  Several  precursors  were  used  in  the  formation  of 
C5H5'*’  ions  and  the  results  of  the  reactivity  studies  indicated  the  possibility 
of  four  different  isomers  of  C5H5'‘'  (88).  Acetylene  reacted  quite  slowly  with 
the  ions  (87,88)  while  diacetylene  and  aromatics  with  side  chains 
reacted  at  an  appreciable  rate  (87).  As  a  result,  it  has  been  suggested  that 
soot  nucleation  may  proceed  by  adding  a  few  large  molecules  rather  than 
through  addition  of  many  smaller  ones  (76,87).  No  definitive  assignment  of 
C5H5'*’  isomeric  structure  was  possible  in  the  earlier  studies  (87,88).  Proton- 
transfer  reactions  involving  one  relatively  unreactive  C5H3'*’  isomer  gave  (87) 
a  proton  affinity  of  227.9  *  0.3  kcal/mol  for  the  neutral  which  remained 
after  proton  transfer.  When  combined  with  estimates  of  the  heats  of  formation 
of  possible  species,  the  results  were  consistent  with  (but  did  not  con¬ 
clusively  prove)  a  vinyl  cyclopropenyllum  form  for  the  CgHj"*”  Isomer. 

A  number  of  theoretical  studies  have  examined  C5H5'*’  structures®  (59- 
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62,97).  Schleyer  and  co-workers  located  two  minimum  energy  forms  on  the 
potential  surface,  the  more  stable  one  corresponding  to  planar  cyclo- 
pentadienyl,  the  other  one  to  a  square-based  pyramid  structure  (60,61). 
Similar  results  were  reported  by  Stohrer  and  Hoffman  (59),  although  they 
proposed  the  pyramidal  structure  as  the  more  stable  form  compared  to  the 
planar  cyclopentadienyl.  Calculations  carried  out  under  this  contract  are 
reported  in  Section  IIF. 

In  this  work,  the  kinetics  and  reaction  mechanisms  of  C^H^^  and  ions 

produced  from  different  precursors  were  studied  with  acetylene  and  diacetylene 
in  order  to  identify  different  structures  and  to  obtain  rate  coefficients. 
Also,  the  energetics  of  C5H5‘*‘  formation  from  two  different  precursors  were 
studied  to  investigate  the  different  pathways  of  formation  which  could 

lead  to  different  structures  as  identified  by  different  reactivity. 

2.  Results 

a.  Reactions  of  C5H3'^« 

A  very  low  number  of  C5H3'*’  ions  were  produced  from  2,4-hexadiyne 
at  electron  energies  above  30  eV.  Reactivity  of  these  ions  was  monitored  at 
an  electron  energy  of  50  eV  although  the  ton  intensity  was  still  very  low. 
Even  though  almost  all  of  the  05113'^  ions  produced  reacted  with  the  2,4- 
hexadiyne  precursor,  no  reaction  was  observed  with  either  C2H2  or  C^H2. 
Reactions  of  C5H3‘*'  with  2,4-hexadiyne  were: 


C5H3++  CgHg  ^  C6H6++  [C3H3]  (71) 

C5H3‘^  +  CgHg  +  CgHy'*'  +  (C2H2]  (72) 

C5H3‘*-  +  CgHe  C5H5'^  +  [C5H4]  (73) 

C5H3'*'  +  CgHg  +  CyHy'*' +  [C4H2I  (74) 

CgH3'^  OgHg  >  ^11^9  (75) 

CgHy  +  >  ^15^13  (76) 


*^5^3^  ions  produced  as  products  of  the  reaction  of  C2H3''‘  +  C^H2  were  100 
percent  reactive  with  both  propargyl  iodide  and  dlacetylene  (the  neutrals 
present  in  the  reaction  medium),  reactions  with  C^H2  were: 


LOO 


C5H3'''  + 

C4H2 

-  C7H3'^  +  IC2H2] 

(77) 

C5H3  + 

C4H2 

>  C9H5''’ 

(78) 

CyH3''‘  + 

C4H2 

(79) 

G9H3'''  + 

C4H2 

>  Ci3H7'^ 

(80) 

Some  of  Che  product  ions  were  observed  to  react  further  with  propargyl  iodide 
by  displacing  atomic  iodine: 

CyHs'*’  +  C3H3I  +  +  I 

C9H5‘'’  +  C3H3I  >  +  I 

C3H3''’  reactions  with  propargyl  iodide  were: 

C5H3+  +  C3H3I  >  C8H6+  +  I 
^8^6^  ^11%^  ^ 

The  CgHg'*’  ion  reacted  further  with  C^H2J 

Ion  intensity  versus  time  curves  for  the  C5H3'''/(C4H2  C3H3I)  system  are 
shown  in  Figure  29.  The  decay  of  05113'^  ions  seen  in  Figure  29  involved 
reaction  with  both  C^H2  and  propargyl  iodide. 

The  procedure  for  rate  coefficient  determination  used  in  previous  studies 
of  C3H3'*’  and  C3H3'*’  ion/molecule  Reactions  (42)  required  subtraction  of  the  ob¬ 
served  rate  constant  for  Che  reaction  of  the  ion  with  Che  precursor  neutral 
from  the  observed  rate  constant  for  the  sum  of  the  reactions  with  precursor 
and  reactant  neutrals.  In  this  case,  however,  the  reactant  neutral  (C^H2)  had 
to  be  Included  in  the  reaction  mixture  In  order  to  produce  the  ion  of  interest 
(C3H3''’)»  and  thus  the  subtraction  procedure  could  not  be  used.  An  alternative 
method  for  rate  constant  determination  was  thus  required.  If  all  other  ions 

are  ejected  from  the  FTICR  cell,  the  decay  of  C,.H3'*’  as  a  function  of  time  is 
.  .  -(n„k„+njkj)t 

given  by  [C5H3'^]  -  [C5H3]oe  PT  “  a  ^ 

products  of  the  number  densities  and  ion/molecule  rate  coefficients  for 


(83) 

(84) 


(81) 

(82) 
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Normalized  Ion  Intensity  (orb.  units) 


Time  (s) 


Figure  29.  Reactions  of  C5H3'*'  with  C^H2.  C5H3‘''(  O  )>  C7H3''‘(+),  C9H5"'’(  0  ) 

and  (A).  Disappearance  of  C5H3'*’  and  product  ions  includes 

reactions  with  propargyl  iodide.  ^5113^  ions  were  produced  from 
the  reaction  of  C3H3'*"  with  diacetylene  within  30  ms  reaction 
time.  C3H3'*’  ions  were  produced  from  propargyl  iodide  by  charge 
transfer  reactions  with  Xe.  PC3H3I  ”  X  10  ^  torr,  Pc^H2  ” 

X  10“^  torr,  Pxe  “  5.4  X  10“®  torr. 
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propargyl  iodide  and  diacetylene,  respectively.  The  quantity  +  ’^p'^p 

thus  be  determined  from  the  slope  of  a  semilog  plot  of  decay  as  a  func¬ 

tion  of  time.  At  short  reaction  times,  the  following  expressions  hold  true: 

,  -(n  lo+nHk:j)t  ^  -(n  k_+n.ko)t 

dD/dt  =  ndkd[C5H3+]oe  d  d  ^  ^ 

where  D  and  P  refer  to  product  ions  of  the  reaction  with  diacetylene  (CyH^'*’ 
and  C9H5''’)  and  with  propargyl  iodide  (CgH^'*’),  respectively.  Thus,  the  ratio 
n^kg/upkp  was  calculated  from  the  ratio  of  slopes  of  product  formation  as  a 
function  of  time.  Next  was  obtained  using  the  calculated  sum  and  ratio 

of  the  two  rate  constants.  Finally  the  absolute  rate  constant  was  determined 
following  the  usual  procedure  (42).  It  was  found  that  ions  reacted  with 


C^H2  with  a  rate  constant  of  (5.6  i  1*7)  X  10 


-10 


cm- 


^/s. 


b.  Formation  of  CgHg'*'- 

CgHg"*"  ions  produced  from  different  precursors  exhibited  behavior 
indicative  of  both  reactive  and  unreactive  populations  toward  the  precursors 
and  the  reactant  neutrals.  Figures  30  and  31  show  C^Hg'*’  normalized  ion  in¬ 
tensity  versus  time  curves  for  two  (1,3,5-cycloheptatriene  and  l-penten-3-yne) 
of  the  five  precursors  used  for  reactions  with  precursor  neutral,  diacetylene 
and  acetylene.  Following  an  exponential  decay  indicative  of  pseudo-first 
order  kinetics,  a  substantial  fraction  of  unreactlve  ions  remains  at  long 
reaction  times.  Table  25  shows  the  method  of  preparation  of  Cgllg'*’  ions  from 
different  precursors  and  the  percent  of  unreactive  ions  remaining  at  long 
reaction  times. 

Two  (norbornadlene  and  cycloheptatriene)  of  the  five  precursors  mentioned 
above  produced  parent  ions  which  could  lead  to  two  different  channels  of 

CgHg"*"  ion  formation.  Thus,  CgHg"*"  ion  formation  pathways  have  been  studied 
somewhat  in  more  detail  for  these  two  precursors.  When  different  charge 
transfer  gases  were  used  for  ionization  of  norbornadlene,  different  percen¬ 
tages  of  reactive  CgHg"*"  were  observed,  as  noted  in  Table  25.  Also  the  abun¬ 
dances  of  CyHy"*"!  and  ions  were  measured  following  charge  transfer 

ionization.  Different  behavior  was  observed  for  compounds  with  ionization  po¬ 
tentials  in  the  range  of  14  to  16  eV.  Table  26  shows  this  effect  for  two 
reagent  gases  (Kr"*"  and  N2^)  which  have  ionization  energies  of  14.0  and  15.7  eV 
respectively. 
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NormoJized  Ion  Signal  (orb.  units) 


Time  (s) 


Reactions  with  Cycloheptstrlenel  O  ). 

-^o  CcHu'*’  ^0“  Curves  Cycloheptatrlene  and 

Figure  30.  C5H5  Acetylene(-»-) ,  and  cycx  k  .  ^  ^ 

Cycloheptatrlene  produced  from  cycloheptatrie 

Dlacetylene(O).  C5H5^  ions  were  p  ^  ^  ^  3  ,  ,0'^  torr. 

Charge  transfer  reactions  -Uh  -  g 

1.1  X  10'*  'C4H2 


i’CjHj* 

«  #N*0 


- -  -  - 

^  p  +  -  1.3  X  10"' 

*  , 

6.6  X  10"'  torr,  P^r  * 


10"^  torr. 
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Time  (s) 


Figure  31.  C5H5'*’  Ion  Decay  Curves  for  Reactions  with  l-Penten-3-yne(  □  ),  l- 

Penten-3-yne  and  Acetylene(+),  and  l-Penten-3-yne  and 
Diacetylene(  0  )•  ^5^5^  ions  were  produced  from  l-penten-3-yne  by 

charge  transfer  reactions  with  Xe'*’.  H  ^  ^  ^  10~^  torr. 


1.4  X  10"^  torr,  P, 


•  C5H6  ■ 

8.4  X  lO""^  torr,  Pj^^ 


5.4  X  10"^ 


TABLE  25.  PERCENTAGES^  OF  UNREACTIVE  C5H5'''  FOUND  FROM  VARIOUS  PRECURSORS  BY 
CHARGE  TRANSFER  CHEMICAL  IONIZATION  MONITORED  BY  OBSERVING  REACTION 
WITH  THE  PRECURSOR  NEUTRAL. 


Precursor  neutral 

Ionizing  technique 

Percent  unreactive  C^Hj"*" 

dicyclopentadiene 

N.j'^/Ar"*"  charge  transfer 

20 

cyclopentadiene 

Xe"^  charge  transfer 

17 

l-penten-3-yne 

Xe"*"  charge  transfer 

65-70 

cycloheptatriene 

Ar"*"  charge  transfer 

0-5 

norbornadiene 

N2^/A.r"*‘  charge  transfer 

18-20 

Kr"^  charge  transfer 

5-10 

^  Estimated  error 

is  ±5%. 

TABLE  26.  CHANGES  IN  ION 
NORBORNADIENE. 

ABUNDANCES  AT  TWO 

DIFFERENT  IONIZATION 

ENERGIES  FOR 

Ionizing  energy/eV 

Abundances  as  a  fraction 

of  total 

ion  signal 

(ionizing  technique) 

CvHg'*' 

C7H7+ 

^5H6"‘ 

^5^5'*' 

14.0 

0.07 

0.45 

0.32 

0.16 

(Kr"*"  charge  transfer) 

15.7 

0.09 

0.28 

0.16 

0.47 

(N2''’  charge  transfer) 

To  further  investigate  this  behavior,  the  relative  ion  abundances  versus 
electron  energy  for  cycloheptatriene  and  norbornadiene  were  obtained  (see 
Figures  32-35).  Electron  impact  ionization  was  used  for  this  study,  there¬ 
fore,  the  energy  scale  in  Figures  32-35  should  be  considered  as  approximate 
with  at  least  ±1  eV  uncertainty. 
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Relative  Intensity  (arb.  units) 


Figure  32.  Relative  Intensities  of  C7Hg‘’‘(+)  and  C7H7'*‘(  □  )  Ions  Produced  from 
Cycloheptatrlene  as  a  Function  of  Electron  Impact  Energy. 
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Relative  Intensity  (orb.  units) 


Figure  33. 


ReUtl,e  intensities  of  and  O  )  !<>"»  fro. 

Norbornadiene  as  a  Function  of  Electron  Impact  Energy. 


108 


Relative  Intensity  (orb.  uni 


Figure  34.  Relative  Intensities  of  C5H^^(  0  )»  C5H5^(+),  and  C3H3^(  D  )  Ions 
Produced  from  Cycloheptatriene  as  a  Function  of  Electron  Impact 
Energy. 
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Relative  Intensity  (orb.  units) 


5  15  25  35  4-5 

Energy  (eV) 


Figure  35. 


Relative  Inteneltles  of  0  ).  C5H5*(^)  •  °  ’ 

Rrodoeed  fro.  Norbotnedlene  as  a  Function  of  Electron  Impact 

Energy. 
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c.  reactions  with  diacetylene. 

After  ejection  of  all  ions  except  following  charge  transfer 

chemical  ionization  of  a  mixture  of  diacetylene  and  a  precursor  compound,  the 
ion/molecule  reactions  as  a  function  of  time  were  monitored.  Independent  of 
the  precursor  used,  the  main  reaction  was  the  addition  of  C^H2  to  to 

produce  the  CgHy"*"  ion. 

Other  minor  reactions  observed  were: 


+  C^H2  +  "h  [C2H2]  (87) 

CjHg'*'  +  C^H2  -*■  CyHy'*'  +  [  C2  ]  (88) 

+  C^H2  ^11^"*"  (89) 

CgH-y"^  +  ^  ^13^"^  (90) 


For  C5H5'''  ions  produced  from  cyclopentadiene,  the  C2  addition  Reaction  (87) 
was  not  observed  with  C^H2* 

The  rate  coefficient  for  the  disappearance  of  05115^  was  calculated  as 
described  earlier  (42)  and  the  values  found  for  ions  produced  from 

different  precursors  are  given  in  Table  27. 

d.  OgHg'*’  reactions  with  acetylene. 

ions  formed  from  four  of  the  five  precursors 
(cyclopentadiene,  l-pentea-3-yne,  norbornadiene,  cycloheptatriene)  reacted 
very  slowly  with  C2H2 ,  producing  very  small  amounts  of  C2H7^  and  CgHg"*"- 

C5H5'''  +  C2H2  -*■  CyHy'*'  (91) 

+  ^2^2  ^  CgHg"^  (92) 

Rate  constants  for  the  05115'*’  +  C2H2  reaction  were  about  one  order  of  magnitude 
less  than  those  for  reactions  with  C^H2  (see  Table  28). 
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TABLE  27.  RATE  COEFFICIENTS  FOR  THE  REACTION  OF  C5H5'''  IONS  FROM  DIFFERENT 
PRECURSORS^  WITH  DIACETYLENE. 


precursor 

Absolute  rate  coef f icient/ ( 10~^®  cm^/s) 

cyclopentadiene 

1.0  ±  0.5 

l-penten-3-yne 

2.0  ±  1.0 

norbornadiene 

2.9  ±  1.4 

cycloheptatrlene 

3.3  ±  1.9 

dlcyclopentadiene 

1.6  ±  0.8 

^  C5H5'^  tons  were  produced  by  charge  transfer  chemical  ionization  using 

different  gases  as  shown  in  Table  25. 


TABLE  28.  RATE  COEFFICIENTS  FOR  THE  REACTION  OF  ACETYLENE  WITH  C5H5''‘  IONS 
PRODUCED  FROM  DIFFERENT  PRECURSORS®. 


precursor^  Absolute  rate  coef f icient/(10“^^  cm^/s) 


cyclopentadiene 


(4.8  ±  1.9) 


l-penten-3-yne 


(2.7  ±  1.0) 


norbornadiene  (3.1  i  1.3) 


cycloheptatrlene  (1.8  i  0.7) 

®  ions  were  produced  using  charge  transfer  chemical  ionization  using 

different  gases  as  shown  in  Table  25. 

^  Ion  signal  for  dlcyclopentadiene  was  too  small  to  produce  reproducible 
results  for  the  rate  coefficient. 


112 


3-  Discussion 


Two  different  structures  have  been  identified  earlier 

(85,86)  by  measuring  the  heats  of  formation  of  the  C5H5''’  ions  from  1,3- 

pentadiyne  and  2,4-hexadiyne.  In  the  work  reported  here,  C5H3'*'  ions  were 
produced  from  2,4-hexadiyne  and  from  reaction  of  C2II3'*’  with  diacetylene. 
Although  ions  produced  either  way  were  reactive  toward  their  precursor 

neutrals,  C5H3'*'  ions  produced  from  2,4-hexadiyne  were  unreact ive  with  both 

C^H2  and  C2H2,  while  those  formed  as  ion/ molecule  reaction  products  reacted 
with  C^H2  with  the  relatively  high  rate  constant  of  (7.79  ±  3.26)  X  10”^® 
cm^/s.  Dannacher  et  al.  (85)  suggested  the  structures  CH^-C  =C-C rC"*"  and 
HC=C-C=C-CH2'’’  for  the  C5H3'*'  ions  produced  from  2,4-hexadiyne  and  1,3- 

pentadiyne  respectively.  Experimental  values  of  the  heat  of  formation  of 
these  ions  showed  that  the  former  had  a  higher  heat  of  formation  (1431 
kj/mole)  (86)  than  the  latter  (1317  kj/mole)  (85).  In  this  study,  higher  re¬ 
activity  of  CgHj'*’  ions  produced  by  the  C3H3'*'  +  C^H2  reaction,  compared  to 

those  produced  by  electron  ionization  of  2,4-hexadiyne  suggests  that  a  higher 
energy  isomer  of  is  formed  in  the  former  case.  It  has  been  found  for 

small  hydrocarbon  ions  that  enhanced  reactivity  often  correlates  with  a  higher 
heat  of  formation. 

b.  ion  structures  produced  from  different  precursors. 

CgHg'*’  ions  produced  from  different  precursors  all  exhibited  at 
least  two  populations,  one  reactive,  other  unreactive.  This  behavior  for  nu¬ 
merous  other  ionic  reactants  has  been  used  previously  (79,98-100)  to  argue  for 
the  existence  of  at  least  two  different  structural  isomers;  one  reactive  and 
one  nonreactlve.  These  results  agree  with  the  earlier  reactivity  study  of 
ions  from  different  precursors  (88) .  To  distinguish  the  reactive  C^Hg'*' 
structures,  their  reaction  mechanisms  and  reaction  rates  with  acetylene  and 
dlacetylene  were  studied.  Although  they  reacted  with  their  precursor  neutrals 
with  different  rates,  they  all  reacted  with  acetylene  and  dlacetylene  at  simi¬ 
lar  rates  (within  the  experimental  error  limits) ,  as  shown  in  Tables  27  and 
28.  Thus  no  isomeric  differentiation  based  on  reactivity  can  be  made.  Rate 
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coefficients  for  the  reaction  with  acetylene  were  similar  to  those  determined 
earlier  (87,88),  i.e.  in  the  range  of  10”^^  cm^/s.  Reaction  mechanisms  were 
similar  for  all  ions  with  the  exception  of  cyclopentadiene,  which  didn't 

give  C2  addition  reactions  with  C^H2*  More  definitive  identification  of  C5H5'*’ 
structures  may  be  forthcoming  when  better  theoretical  calculations®  have  been 
reported. 

To  compare  the  behavior  of  ions  produced  from  1 , 3 ,5-cycloheptatriene  and 
norbornadieue  precursors  (both  having  the  formula  Cytig),  relative  abundances 
as  a  function  of  electron  energy  shown  in  Figures  32-35  were  obtained.  As 
shown,  the  curves  for  CyHy"^  (Figures  32-33),  ^5115'*’  and  C3H3'*’  (Figures  34-35) 
are  very  similar  both  in  terms  of  energetics  and  of  general  shape.  Slight 
differences  in  their  relative  intensities  are  probably  due  to  the  normaliza¬ 
tion  procedure  used  in  each  case.  Increase  in  relative  intensity  of 
ions  in  both  cases  in  the  energy  range  of  15-25  eV  and  the  fact  that  the  CyHy^ 
curve  becomes  flat  in  this  energy  range  implies  that  most  of  the  ions 

are  produced  from  CyHy'*’.  Similar  behavior  exists  for  ions,  i.e.,  the 

intensity  of  C3H3'^  increases  near  22  eV  (on  the  energy  scale  shown  on  Figures 
34-35)  where  the  rate  of  increase  in  intensity  of  C5H5'*'  lessens.  This  implies 
that  formed  from  these  two  compounds  behaves  similarly,  i.e.,  follows 

the  same  fragmentation  pathways  in  the  same  energy  range.  Similar  behavior 
for  CyHy'^  ions  from  these  two  compounds  has  also  been  reported  earlier  as  a 
result  of  a  collisional  activation  study  (101).  Formation  of  CjHg'*’,  both  in 
terms  of  energetics  and  relative  intensity,  is  different  for  cycloheptatriene 
and  norbornadiene,  indicating  that  C7H3''’  ions  formed  from  these  compounds  by 
electron  Impact  have  different  structures.  Based  on  the  higher  relative  in¬ 
tensity  of  CgHg'*’  for  norbornadiene,  there  exists  the  possibility  of  formation 
of  C5H5'*’  ions  from  the  C5H3'*’  channel  too.  In  fact,  C5H5'''  ions  formed  from  cy¬ 
cloheptatriene  in  the  energy  range  of  12  to  50  eV  show  similar  behavior  in  re¬ 
activity,  i.e.,  the  percent  of  unreactlve  C5H5'''  remains  constant  at  5  percent 
while  in  the  norbornadiene  case,  percent  of  unreactivity  increases  to  20  per¬ 
cent,  indicating  that  formed  from  the  channel  has  a  high  percent 

of  unreactlve  Isomer.  This  trend  is  also  apparent  from  the  results  shown  in 
Tables  25  and  26  for  norbornadiene,  which  indicates  that  when  Ar'*'/N2^  replaces 
Kr‘*‘  as  the  charge  donor,  the  unreactive  population  of  the  observed  in¬ 
creases  while  there  is  a  corresponding  decrease  in  the  yield  of  the  ion. 


Similar  results  have  also  been  reported  earlier  for  norbornadiene  (80,102)- 

In  general,  a  similar  trend  in  ion/molecule  reactivity  with  acetylene  and 
diacetylene  has  been  found  for  0^115^  ions  as  was  seen  (Reference  42  and 
Section  IIIC.)  for  C3H3'*’  ions:  little  or  no  reactivity  with  acetylene,  but 

extensive  and  rapid  reactivity  with  diacetylene.  The  results  of  bimolecular 
reactivity  studies  on  both  C3H3‘''  and  C5H5''’  suggest  that  the  proposed  ion/mole¬ 
cule  condensation  reactions  with  diacetylene  in  the  ionic  soot  formation 
mechanism  are  quite  credible  while  the  proposed  sequential  acetylene  addition 
reactions  should  be  reconsidered.  These  experiments  and  those  reported 
earlier  (42)  were  carried  out  at  relatively  low  pressures  (<  5  X  10~^  torr) . 
It  is  thus  quite  possible  that  third-body  collisions  in  atmospheric  pressure 
flames  can  stabilize  a  larger  fraction  of  the  C7H7'*’  reaction  complexes  and 
result  in  increased  reactivity  with  acetylene.  In  fact,  such  collisional 
stabilization  of  the  association  complexes  for  the  reactions  of  C3H3''’  and 
*^4^4*^  with  C2H2  has  been  shown  to  occur  in  higher  pressure  selected  ion  flow 
tube  (SIFT)  studies  (103,104). 

E.  REACTIONS  OF  GASEOUS  CyHy'*’  IONS 

1.  Introduction 

The  CyHy'*’  ion  has  been  proposed  as  one  species  Involved  in  the 
ion/molecule  mechanism  of  soot  formation  (14a, b).  According  to  the  proposed 
ion/molecule  mechanism,  forms  in  reactions  with  acetylene. 

However,  work  reported  in  Section  HID.  on  the  reactions  of  C5H5''’  with  acety¬ 
lene  did  not  reveal  significant  formation  of  CyHy'*’.  Much  higher  concentra¬ 
tions  of  than  observed  experimentally  were  predicted  (14a)  by  a  simple 

model  which  assumed  a  benzyl  structure  for  the  ion.  To  account  for  this  dif¬ 
ference  between  experiment  and  theory,  it  was  suggested  (14a)  that  the  actual 
ton  in  the  flame,  presumably  formed  by  ion/ molecule  reactions  of  C3Hx^  and 
C5Hx^  ions,  is  some  other  CyHy"*"  isomer  less  stable  than  benzyl.  To  help 
understand  this  question  we  have  extended  our  earlier  studies  on  hydrocarbon 
ions  which  might  be  important  In  an  ionic  mechanism  of  soot  formation^ 
(42,76,96,105)  to  follow  the  reactions  of  ions  with  certain  flame 

neutrals  and  to  try  to  Identify  different  Isomers. 
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There  have  been  a  significant  number  of  both  theoretical  and  experimental 
studies  relating  to  structural  characterization  of  ions.  Theoretical 

studies  (106,107)  which  reported  heats  of  formation  for  C7H7''’  ions  indicated 
chat  tropylium,  I  ( AHf  =*  207.9,  195.6  kcal/mol)  (106,107)  and  benzyl,  II  (AHf 
=»  217.1,  220.4  kcal/mol)  (106,107)  are  the  most  stable  structures  with  an 

II 

activation  energy  of  32.7  kcal/mol  (107)  for  II  — >  I  isomerization.  Experi¬ 
mental  measurements  of  ionization  potentials  for  tropylium  and  benzyl  radicals 
by  monoenergetic  electron  impact  gave  values  of  aHj  *  209  and  211  kcal/mol 
(108-110)  respectively.  Studies  of  metastable  ions  (111-113),  collislonal 
activation  (101,114,115)  gas  phase  radiolysis  (116-119),  photoionization/pho- 
todlssociation  (95,120-124)  photoelectron  photolon  coincidence  (PEPICO) 
(125,126),  and  ion/molecule  Reactions  (79,99,100,127,128)  of  07117"^  ions  from 
different  precursors  were  used  to  characterize  the  different  structures. 
General  experimental  evidence  including  Che  determination  of  appearance 
potentials  and  heats  of  formation  (122-126)  supported  the  idea  that  these  ions 
exist  in  the  gas  phase  with  at  least  two  structures;  tropylium,  I  and  benzyl, 
II.  (Jnreactive  and  reactive  07117^  observed  in  ion/molecule  reactions  have 
been  attributed  to  structures  I  and  II,  respectively  (79,99,100,127,128). 
Furthermore,  it  was  proposed  that  there  is  an  equilibrium  between  these  two 
structures  (114,129)  similar  to  that  proposed  between  the  toluene,  III  and 
cycloheptatrlene,  IV  (114,80)  ions. 

<0>-CH3 

III 

Several  groups  studied  the  effect  of  internal  energy  above  the  €7117^  formation 
threshold  on  the  predominance  of  one  structure  over  the  other  (80,112,121, 
128).  Their  results  in  general  were  supportive  of  the  equilibria  between 
different  structures  mentioned  above.  Although  McLafferty  et  al.  (101,112) 
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reported  evidence  for  stable  toiyl  and  norbornadienyl  structures  in  earlier 
collisional  activation  studies  no  definitive  indication  for  stable  structures 
other  than  benzyl  and  tropylium  was  provided  in  later  work  (79,80,114,115). 
In  this  work  CyHy'*’  ions  were  produced  from  three  different  precursors  using 
different  charge  transfer  gases.  Reaction  pathways  and  the  rate  coefficients 
for  the  reaction  of  CyHy"^  ions  with  precursor  neutrals,  acetylene  and  diacety¬ 
lene  near  room  temperature  were  investigated.  Results  are  discussed  in  rela¬ 
tion  to  previous  work  involving  CyHy'*’  structures  and  the  proposed  ionic  route 
to  soot  formation. 

2.  Results 

a.  Effect  of  Precursor  Neutrals. 

CyHy"'"  ions  produced  from  different  precursors  were  found  to  have 
both  reactive  and  unreactive  populations.  Figures  36,  37,  and  38  show 
normalized  ion  intensities  of  CyHy"*"  and  product  ions  vs  time  for  reactions 
with  precursor  neutrals;  toluene,  norbornadiene  and  cycloheptatriene, 
respectively.  Following  an  exponential  decay  indicative  of  first-order 
kinetics,  a  substantial  fraction  of  unreactlve  ions  remains  at  long  reaction 
times.  Table  29  shows  the  charge  transfer  agents  used  to  produce  CyHy"*"  ions 
from  different  precursors,  the  maximum  internal  energy  imparted  to  the  ions, 
and  the  percent  of  unreactlve  ions  remaining  at  long  reaction  times. 

The  major  reaction  of  CyHy^  produced  by  charge  transfer  from  toluene  with 
the  toluene  neutral  gives  CgHg'*'  (Reaction  (93)).  When  norbornadiene  is  used 
as  a  precursor  and  neutral  reactant.  Reaction  (94)  is  observed. 

CyHy'^  +  CyHg - >  CgHg'^  +  CgH^  (93) 

+  CyHg - >  CgHg"*'  +  C5Hg  (94) 

CyHy^  from  norbornadiene  also  reacts  to  give  small  amounts  of  higher  molecular 
weight  products.  On  the  other  hand,  CyHy'*’  produced  from  cycloheptatriene 
reacts  with  its  precursor  neutral  to  give  significant  amounts  of  higher 
molecular  weight  product  Ions  along  with  smaller  amounts  of  CgHg^  and  CgHg"** 
product  ions: 
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NORMALIZED  INTENSITY 


•w - 1 - 1 - 1 - 1 - 1  I 

0  2  4-6 


TIME  (s) 


Figure  36.  Reaction  of  Produced  from  Toluene  by  Kr+  Charge  Transfer 

with  its  Precursor  Neutral.  Ptoluene  *  ^  torr,  - 

6.3  X  10"^  torr.  (  ■»• )  and  CgHg'*’  (□) . 
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NORMALIZED  INTENSITY 


Figure 


0  0.2  0.4  0.6  as  1  1.2 
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37.  Reaction  of  CyHy"*"  Produced  from  Norbornadiene  by  Xe"^  Charge 
Transfer  with  its  Precursor  Neutral.  Pnorbornadlene  “  5.4  x  10 
torr,  P^g  ■  4.3  X  10  ^  torr.  (D)  (‘v>)* 
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NORMALIZED  INTENSITY 


Figure  38.  Reaction  of  CyHy'*’  Produced  from  Cycloheptatriene  by  Xe'*'  Charge 

Transfer  with  its  precursor  Neutral.  ^cycloheptatriene  “  3.6  x 
10“^  torr,  p^g  -  4.6  x  10~^  torr.  (Q),  CgHg'*’  (+■) ,  CgHg"*" 

(O) ,  and  higher  molecular  weight  product  ions  (A). 
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TABLE  29.  PERCENTAGES  AND  MAXIMUM  INTERNAL  ENERGIES  OF  UNREACTIVE  C^Hy'*'  IONS 
FORMED  FROM  VARIOUS  PRECURSORS  BY  DIFFERENT  IONIZATION  METHODS 


MONITORED 

BY  OBSERVING  REACTION 

WITH  THE  PRECURSOR 

NEUTRAL . ^ 

Precursor 

Ionizing  Charge 

Maximum  Internal 

Percent 

Neutral 

Transfer  Agents 

Energy  (eV)^ 

Unreact ive 

Toluene 

Xe+ 

1.4 

45-50 

Kr+ 

3.3 

25-30 

Ar"*" 

5.0 

25-30 

Norbomadiene 

C4H2'^ 

0.9 

70-75 

Xe+ 

2.5 

30-35 

Kr'*’ 

4.4 

25-30 

Cycloheptat riene 

1.1 

70-75 

Xe+ 

2.8 

35-40 

Kr'*' 

4.6 

35-40 

*  For  those  ions  produced  by  charge  transfer,  reactions  rfith  both 

precursor  neutral  and  diacetylene  were  monitored. 

^  Based  on  CyHy"*”  appearance  energies  given  by  Traeger  and  McLouglin  (1978). 


C7H7+  +  C^Hg - >  CgHg-^  +  CgHg  (95) 

CyHy"^  +  ^7^8  - - ^  CgHg"^  +  CjHg  (96) 

C7H7+  +  C7Hg - >  (97) 

0787'*'  +  C7Hg - >  +  CgHg  (98) 

C7H7'^  +  Cydg  — — — >  ^12^11"^  ^2^4  (99) 

0787'’'  +  C7Hg - >  CigHii'*'  +  CH4  (100) 

C7H7'''  +  C7Hg - >  +  2H2  (101) 

C7H7'^  +  C7Hg  — — — >  ^14^13"^  ^2  (102) 


b.  Deuteration  Studies. 


In  order  to  further  investigate  the  mechanisms  of  some  of  the 
reactions  shown  above,  deuterated  toluene  was  used  both  to  produce  for 
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reactions  with  neutral  precursors  and  as  a  neutral  reactant  for  reactions  of 
produced  from  different  precursors.  Reactions  of  with  toluene  and 

norbornadiene  (NBD)  were  respectively; 

CyDy"*"  +  CyHg  (toluene) - >  C8H7D2''’  +  CgHDg  (103) 

+  C7H3  (NBD) - >  0907112'*'  +  CgHg  (104) 

When  C7D7'*’  reacted  with  cycloheptatriene  (CHT),  the  major  product  observed  was 
unreactive  C7H7'*’  due  to  Reactions  (105)  and  (106). 

C7D7'''  +  C7H8  (CHT) - >  C7H8'‘'  +  C7D7 

C7Hg''’  +  C7Hg  (CHT) - >  C7H7''‘  (unreactive)  +  C7H9 

Other  reactions  observed  were: 


C7D7‘''  +  C7H3  (CHT) 

- >  C9D3H3'’'  +  C5H5D 

(107) 

C7D7'''  +  C7H3  (CHT) 

- >  CioDeHs"^  +  C4H5D 

(108) 

+  C7Hg  (CHT) 

- >  CuDeHs''  +  C3H5D 

(109) 

C7D7''’  +  C7H8  (CHT) 

— — >  0^^207!!^'*'  +  ^2^4 

(110) 

C7H7'*'  ions  produced  from  toluene  reacted  with  deuterated  toluene  as  shown  by 
Equation  (111): 


(105) 

(106) 


C7H7'*'  +  0703 - >  C3D7H2'‘'  +  C3H5D  (111) 

Very  low  intensities  of  C3D7H2^  were  also  observed  in  the  mass  spectra  of 
C7H7'*’  Ions  produced  from  cycloheptatriene  and  norbornadiene.  The  rate 
coefficients  for  the  disappearance  of  C7H7'*’/C7D7"*’  by  reactions  with  precursor 
neutrals  were  calculated  as  described  earlier^  (42)  and  the  values  are  shown 
in  Table  30. 


c.  O7H7'*'  Reactions  with  Acetylene. 

C7H7'*'  ions  produced  by  Xe"*"  dissociative  charge  transfer  of 
norbornadiene  and  cycloheptatriene  reacted  slowly  with  acetylene  to  produce 
very  small  amounts  of  C9H9'*’* 
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TABLE  30.  RATE  COEFFICIENTS  FOR  DIFFERENT  C7H7‘'‘/C7H8  REACTION  SYSTEMS. 


k/l0“^®(cm^/ s) 

from 

CjHj'*'  from 

C7H7'^  from 

Neutral 

Toulene-dg 

NBD 

CHT 

Toluene-dg 

(3.5  ±  1.0) 

(0.8  i  0.3) 

(1.7  ±  0.5) 

NBD 

(3.2  ±  1.0) 

(5.4  ±  1.7) 

CHT 

(4.6  ±  1.5) 

(8.2  ±  2.6) 

C7H7'^  +  C2H2 - >  CgHg'*'  (112) 

The  following  reaction  was  also  observed  with  the  Ions  produced  from 
cycloheptatrlene ; 


C7H7'*'  +  ^2^2  ^2  (113) 

The  rate  coefficients  for  the  disappearance  of  07117’^  Ions  produced  from 
norbornadlene  and  cycloheptatrlene  were  found  to  be  (1.0  ±  1.0)  x  10~^^  and 
(3.6  ±  1*3)  X  10“^^  cm^/s,  respectively.  However,  when  toluene  was  used  as  a 
precursor  for  production  by  Kr"*"  dissociative  charge  transfer,  no 

reaction  with  acetylene  was  observed. 

d.  C7H7'*’  Reactions  with  Diacetylene. 

C.7H7''’  ions  from  different  precursors  were  reacted  with  both 
diacetylene  and  deuterated  diacetylene  in  order  to  identify  the  reaction 
products  from  those  formed  by  reaction  of  C7H7'*’  with  the  precursor  neutral 
(for  instance  the  reactions  with  dlacetylene  and  cycloheptatrlene  both  give 
CiiHg^).  The  following  reactions  of  C7H7'*’  with  diacetylene  were  observed: 

C7H7‘’'  +  C4H2 - >  CgHj'*'  +  C2H2  (114a) 

C7H7'^  +  C^D2 - >  C9H5D2‘^  +  C2H2  (114b) 
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C7H7+  + 

C4H2 - > 

(115a) 

C4D2 - >  ^11^7°2'*’ 

(115b) 

C9H7'^  + 

C4H2 - > 

(116a) 

'^9^5°2''' 

+  C4D2 - >  C]^3H5D4'*' 

(116b) 

In  the  toluene  case,  the  following  reactions  were  also  observed: 

CgHg  +  C^H2 - y  ^12^9  ^2  (117  a) 

+  C^D2 - >  C^2^7°2'^  ^2*  (117b) 

Ion  Intensity  versus  time  curves  of  (produced  by  Kr"*"  charge  transfer 

from  toluene)  and  product  ions  for  the  reaction  with  diacetylene  are  shown  in 
Figure  39.  The  rate  coefficients  for  the  disappearance  of  by  the  sum  of 

Reactions  (114)  and  (115)  for  different  precursors  are  given  in  Table  31. 


TABLE  31.  REACTION  RATE  COEFFICIENTS  FOR  CyHy'*'  +  C4H2  REACTION. 


Precursor 

k/10~^®(cm^/s) 

Toluene 

(1.50  i  0.50) 

Norbornadlene 

(0.25  ±  0.19) 

Cyclohpetatriene 

(0.75  ±  0.25) 

When  CyHy"^  ions  were  produced  in  a  100  ms  reaction  time  by  Xe'*’  dissocia¬ 
tive  charge  transfer  from  norbornadlene  and  cycloheptatriene  in  the  presence 
of  C4H2,  production  of  unreactive  was  observed  by  the  following 

reactions : 

Xe'*'  +•  C4H2 - >  C4H2'’'  +  Xe  (118) 

^4^2^  +  C7H8 - >  C7H7'''( unreactive)  +  C4H2  +  H.  (119) 

Thus,  to  obtain  accurate  information  on  reaction  kinetics  of  the  C7H7'^  +  ^4^2 
reaction,  the  C4H2^  ion  was  continuously  ejected  when  C7H7''’  ions  were  formed 
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TIME  (9) 


Figure  39.  Reactions  of  C7H7'*’  Produced  from  Toluene  by  Kr"*”  Charge  Transfer 
with  Dlacetylene.  ^toluene  ”  torr,  PC4H2  * 

torr,  -  6.3  X  10"^  torr.  C7H7'*’  (□),  C9H7‘'‘  (+),  (0), 

and  (A). 
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by  dissociative  charge  transfer  using  Xe'*'. 

Also  in  this  study,  the  relative  abundances  of  the  ions  produced  from 
toluene  as  a  function  of  electron  impact  energy  were  determined  for  the 
purpose  of  comparison  with  earlier  results  obtained  for  norbornadiene  and 
cycloheptatriene^ .  Relative  abundances  of  and  as  a  function  of 
energy  for  toluene  are  shown  in  Figure  40.  Results  of  the  earlier  study^  for 
norbornadiene  and  cycloheptatriene  are  also  included  in  the  figure  for 
comparison. 

3.  Discussion 

Since  numerous  earlier  investigations  (79,100,114,121,123)  have 
demonstrated  that  ions  have  at  least  two  stable  structures,  one  reactive 
and  the  other  unreactlve,  the  results  of  this  study  are  interpreted  in  terms 
of  the  possibility  of  other  stable  structure(s) . 

The  nature  of  the  reaction  products  observed  for  the  reactions  of  CyHy'*’ 
ions  with  their  precursor  neutrals,  toluene,  norbornadiene  and  cycio- 
heptatriene,  were  significantly  different,  as  shown  in  Figures  36,  37  and  38, 
respectively.  In  order  to  determine  whether  these  differences  were  due  to  the 
neutral  or  ion  structures,  deuterated  toluene  was  used  to  monitor  the  reaction 
mechanisms  with  different  neutral  precursors.  As  seen  in  Equations  (103-110), 
CyDj'*'  ions  reacted  with  the  precursor  neutrals,  producing  deuterated  analogues 
of  the  same  product  ions  observed  for  each  nondeuterated  ion/precursor 
reaction  system.  This  result  showed  that  the  differences  in  reaction  products 
are  the  result  of  the  different  structures  of  neutral  precursors,  and  not 
necessarily  of  different  ion  structures.  Then  deuterated  toluene  was  used  as 
a  reactant  neutral  for  ions  produced  from  different  precursors.  As  seen 
from  Equation  (111),  CyHy"*"  ions  produced  from  all  three  precursors  reacted 
with  C7D3  to  produce  CgD7H2'’’f  although  very  low  ion  intensities  were  produced 
with  cycloheptatriene  and  norbornadiene  precursors.  The  reaction  mechanism 
which  gives  this  product  ion  is  also  operative  in  Reaction  (103)  and  is 
believed  to  be  indicative  of  the  benzyl,  II,  structure  for  ions 
(79,95,121,128).  This  mechanism  has  also  been  studied  earlier  by  similar 
isotopic  labeling  expex'lments  (79,99).  The  only  important  difference  in  the 
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Figure  40.  05115^  and  C2H3'*’  Ion  Abundances  as  a  Function  of  Electron  Energy. 

CjEj"*"  (X)  and  03113^  (7)  from  Cycloheptatriene,  C5H5'*’  (  □  )  and 
^3^3^  (+)  from  Norbornadlene  and  (O)  and  C3H3'*'  (4)  from 

Toluene. 


behavior  of  reactive  CyHy"*"  ions  produced  from  different  precursors  was  their 
reaction  rates  with  precursor  neutrals.  ions  produced  from 
toluene/deuterated  toluene  reacted  at  similar  rates  with  toluene, 
norbornadiene  and  cycloheptatriene  whereas  CyHy"*"  ions  produced  either  from 
norbornadiene  and  cycloheptatriene  reacted  at  much  slower  rates  with  toluene 
than  with  their  own  precursor  neutrals  (see  Table  30)  . 

Reactions  of  and  product  ions  with  diacetylene  and  deuterodiacety- 
lene  are  represented  by  Equations  (114-116),  independent  of  the  precursor 
used.  An  Interesting  mechanism  is  observed  for  C2  and  addition  reactions 
shown  by  Equations  (114b)  and  (117b)  which  is  indicative  of  an  ion/molecule 
complex  from  which  C2H2  and  H2,  respectively,  rather  than  C2D2  and  D2  are  the 
apparent  neutral  products.  Since  complete  isotopic  scrambling  of  product  ions 
is  not  observed,  preferential  retention  of  C2D2  and  D2  is  probably  due  to  a 
specific  structure  of  the  intermediate  complex  which  does  not  rearrange  fast 
enough  to  give  complete  scrambling  before  falling  apart  to  give  the 
products.  Reaction  rates  with  dlacetylene  followed  a  different  trend  than 
those  with  precursor  neutrals  for  CyHy"**  ions  from  different  precursors.  In 
this  case,  ions  from  toluene  reacted  faster  than  those  from 
norbornadiene  and  cycloheptatriene  (see  Table  31)* 

Different  behavior  was  observed  when  C2H2  used  as  a  reactant 
neutral.  from  norbornadiene  and  cycloheptatriene  reacted  slowly  with 
acetylene  as  shown  by  Equations  (112)  and  (113).  On  the  other  hand,  no 
reaction  products  were  observed  with  CyHy'*'  from  toluene. 

Table  29  shows  the  percent  of  unreactive  CyHy'*’  as  a  function  of  the  in¬ 
ternal  energy  of  the  ion  formed  from  different  precursors.  At  the  lowest  in¬ 
ternal  energies  above  formation  threshold,  CyHy'*’  from  toluene  had  a  much 
smaller  percentage  of  unreactive  isomer  than  CyHy^  from  norbornadiene  and  cy¬ 
cloheptatriene.  Formation  of  unreactive  CyHy"*"  from  all  three  precursors  de¬ 
creased  with  increasing  internal  energy  and  stayed  constant  at  high  enough  en¬ 
ergy.  This  behavior  is  different  from  that  observed  previously  for  toluene 
and  norbornadiene  (121,80)  in  which  the  unreactlve  percent  first  decreased  and 
then  Increased  with  Internal  energy.  The  previous  results  were  Interpreted  in 
terms  of  the  interplay  between  the  two  equilibria  Involving  different  Cyng"*" 
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and  Cytiy''’  structures  (114,118,129,80).  The  different  behavior  observed  in  an 
earlier  ICR  study  (79),  which  is  similar  to  that  observed  in  this  work,  was 

suggested  (114,80)  to  be  the  result  of  the  displacement  of  the  equilibrium 
between  two  CyHy"*"  structures  by  the  removal  of  the  reactive  CyHy'*’  ions  through 
chemical  reaction. 

In  an  earlier  study^,  the  relative  abundances  of  ions  produced  from 

norbornadiene  and  cycloheptatriene  as  a  function  of  ionizing  electron  energy 
were  determined.  Similarities  of  CyHy'*’,  appearance  curves 

were  interpreted  to  imply  that  CyHy"*"  formed  from  these  two  neutral  molecules 
behaves  similarly,  i.e.,  follows  the  same  fragmentation  pathways  in  the  same 
energy  range.  Differences  in  the  appareance  curve  of  for  norbornadiene 

and  cycloheptatriene  were  interpreted  as  suggesting  that  CyHg'*'  formed  from 

these  compounds  by  electron  Impact  may  have  different  structures^.  Extensive 
fragmentation  of  CyHg"^  to  give  was  also  reported  earlier  for 

norbornadiene  (80).  In  this  study,  the  relative  abundances  of  the  ions 
produced  from  toluene  as  a  function  of  energy  were  determined  for  the  purpose 
of  comparison  with  earlier  results  obtained  for  norbornadiene  and  cyclohep¬ 
tatriene.  In  contrast  to  the  observed  similarity  mentioned  above  between 
norbornadiene  and  cycloheptatriene  fragment  ion  appareance  curves,  those  of 
toluene  were  significantly  different  both  in  terms  of  energetics  and  of  gen¬ 
eral  shape.  For  example,  the  onsets  of  formation  for  both  C5H5''’  and 
were  different  than  those  values  observed  for  the  other  two  precursors.  Fur¬ 
thermore,  the  relative  intensities  for  and  03113'^  increased  much  more 

slowly  with  energy  relative  to  those  from  norbornadiene  and  cyclohep¬ 
tatriene.  The  maximum  relative  intensities  for  C5H5'*’  and  C3H3'*’  ions  were 

0.017  and  0.005  for  toluene  and  0.10  and  0.035  for  norbornadiene  and  cyclohep¬ 
tatriene,  respectively  (see  Figure  40).  This  different  behavior  in  the  frag¬ 

mentation  of  CyHy"*"  produced  from  toluene  suggests  the  involvement  of  a  dif¬ 
ferent  structure  (benzyl"*",  n)  for  CyHy"*"  in  the  case  of  toluene. 

Based  on  the  differences  in  behavior  observed  in  this  study  of  CyHy'*’  ions 
produced  from  three  precursors,  the  possibility  of  structures  other  than  ben¬ 
zyl  and  tropylium  must  be  considered.  The  unreactlve  structure  which  was 

observed  in  all  three  cases  has  been  shown  to  be  cyclic  "tropylium"  (I)  by  a 
number  of  previous  experimental  (79,99,100,122-128)  and  theoretical  (106,107) 


L29 


studies.  Similarly,  the  reactive  isomer  produced  from  toluene  was  identified 
as  the  "benzyl”  (II)  structure  numerous  times  previously  (79,99,100,121- 
128).  However,  the  results  of  this  study  have  shown  that  it  behaves  somewhat 
differently  than  the  reactive  07^1^'*’  isomers  produced  from  norbornadiene  and 
cycloheptatriene.  Although  no  definitive  assignment  for  reactive  CyHj'*'  ion 
structure(s)  produced  from  norbornadiene  and  cycloheptatriene  can  be  made  at 
this  point,  the  most  probable  structures  will  be  discussed  considering  the 
energetics  of  C7H7''’  ions. 

a.  Cycloheptatriene. 

AMI  calculations  for  the  C7Hg''‘/C7H7"*’  surface^^  suggest  the  struc¬ 
ture  produced  by  hydrogen  loss  from  the  bicyclo[3.2.0]heptadiene  radical 
cation  (BClT*')  (V),  is  a  possible  stable  CyHy"**  isomer,  accessible  at  about 


2.8  eV  above  the  appearance  potential  of  tropylium  from  cycloheptatriene. 
This  structure  is  in  a  potential  well  which  is  deep  enough  that  ii.  might  rea¬ 
sonably  be  expected  Co  live  for  an  extended  time^^.  It  was  also  suggested^^ 
that  this  structure  is  most  probably  formed  from  cycloheptatriene  by  isomeri¬ 
zation  on  the  radical  cation  surface  followed  by  hydrogen  loss  from  v.  Isom¬ 
erization  of  (BCH-H)'*',  VI  (AHf  »  274.4  kcal/mol)^^  to  tropylium  is  considered 
unlikely  due  to  the  requirements  of  orbital  symmetry,  which  imposes  no  con¬ 
straints  for  the  odd  electron  radical  parent  ion.  Although  the  barrier  to 
form  II  is  lower  than  that  for  V,  it  was  suggested^2  that  the  simplicity  of 
the  rearrangement  to  the  latter,  compared  to  the  convoluted  surface  leading  to 
benzyl  or  toluene  radical  ion,  may  make  formation  of  BCH'*',  V,  a  faster  rear¬ 
rangement  at  high  enough  energy.  Consistent  with  these  calculations,  our  ex¬ 
perimental  results  show  that  the  percent  of  unreactive  C7H7''"  from  cyclohep¬ 
tatriene  decreases  very  rapidly  with  increasing  internal  energy.  This  can  be 
Interpreted  as  a  result  of  Increased  isomerization  of  the  cycloheptatriene 
cation  to  V.  At  high  energy,  when  the  equilibrium  between  these  two  parent 
ions  is  established,  the  percent  of  unreactive  isomer  stays  constant  due  to 
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the  absence  of  isomerization  between  two  CyH^'*’  structures  with  an  energy  bar¬ 
rier  of  14.1  kcal/mol^2.  Other  experimental  evidence  for  the  involvement  of 
different  structures  for  CyHg'*'  ions  formed  from  toluene  and  cycloheptatriene 
comes  from  the  strikingly  different  photodissociation  curves  of  these  ions 
(130). 

b.  Norbornadiene. 

In  spite  of  the  general  similarity  of  the  behavior  of  ions 
from  cycloheptatriene  and  norbornadiene,  differences  in  reactivity  (see  Tables 
29  and  30)  suggest  the  formation  of  a  CyHy'*'  structure  different  from  the 
structure  formed  from  cycloheptatriene.  A  simple  bond  breakage  in 
norbornadiene,  (VII)  can  directly  form  an  ipso-protonated  benzyl  radical. 


VIII 


0^115(11^)012  ion,  in  which  a  1,2-hydrogen  shift  leads  to  the  intermediate,  5- 
methylene-1 ,3-cyclohexadiene  ion,  VIII,  with  a  theoretical  -  219  kcal/mol 
(131).  The  latter  ion  was  proposed  to  account  for  the  ring  expansion  of  the 
toluene  ion  in  the  III  — >  IV  isomerization  process  (131).  The  formation  of 
VIII  from  other  precursor  cations  at  least  with  a  transitory  existence  has 
been  indicated  by  previous  experimental  evidence  (132).  After  the  availabili¬ 
ty  (133)  of  5-methylene-l ,3-cyclohexadiene,  experimental  determination  of  its 
ionization  potential  indicated  a  heat  of  formation  18  kcal/mol  higher  than 
that  of  the  toluene  ion  (134). 

Once  this  intermediate  Ion  is  formed,  it  can  rearrange  to  another  struc¬ 
ture  or  decompose  directly  to  the  CyHy'*’  ion.  In  fact,  it  was  suggested  that 
direct  decomposition  of  this  ion  is  in  competition  with  rearrangement  to  other 
structures,  cycloheptatriene  cation  being  the  most  probable  rearrangement 
product  (135).  The  competition  between  these  two  processes  is  believed  to 
depend  on  the  Internal  energy  with  rearrangement  being  favored  at  lower  ener¬ 
gies  (136).  Thus  at  low  energies,  fragmentation  to  CyHy^  after  rearrangement 
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Co  the  cycloheptatriene  ion  is  expected  due  to  Che  lower  activation  barrier 
for  fragmentation  of  the  cycloheptatriene  ion  (121).  This  is  consistent  with 
our  results,  which  indicate  a  high  percent  of  unreactive  CyHy'*’  ( tropylium"^,  I) 
formed  from  norbornadiene  at  low  energy  (see  Table  28).  Studies  of  Che  photo- 
dissocxction  of  ions  suggested  about  half  of  the  norbornadiene  cations 

convert  to  cycloheptatriene  cations  (130).  On  the  other  hand,  it  has  been 

reported  that  there  is  no  rearrangement  of  this  ion  to  the  toluene  ion  on  a 

time  scale  of  10”^  s  (137). 

photodissociation  spectra  of  n-butylbenzene  and  2-phenylethanol  ions  which 
are  believed  to  form  VIII  were  studied  to  identify  any  possible  rearrangement 
products  (136).  The  results  were  consistent  with  the  retention  of  the  VIII 
structure.  Although  there  was  not  any  evidence  for  ready  interconversion 
among  CyHg^  structures  for  nondecomposing  methyl enecyclohexadiene  ions,  it  was 
suggested  that  the  more  energetic,  decomposing  VIII  ions  may  interconvert  to 
other  structures.  In  this  study,  the  decrease  In  the  percent  of  tropylium 

structure  with  an  increase  in  energy  is  interpreted  as  a  result  of  the  fact 
that  decomposition  of  VIII  becomes  competitive  with  rearrangement  to  IV. 
Thus,  we  suggest  the  CyHy'*’  structure  formed  by  hydrogen  loss  from  VIII  as  the 
reactive  Isomer  formed  from  norbornadiene. 

ICR  studies  (134)  involving  the  formation  of  VIII  from  its  precursor 

neutral  at  30  eV  electron  energy  indicated  ion/molecule  chemistry  different 
from  that  of  toluene.  produced  from  VIII  was  reported  to  have  a  dif¬ 

ferent  rate  coefficient  and  reaction  efficiency  for  the  formation  of  Cgllg"^. 
It  is  Interesting  to  note  that  the  rate  coefficient  reported  (0.6  x  10“^*^ 
cm^/s)  (134)  for  this  reaction  agrees  very  closely  with  the  value  (0.3  ±  0.3) 
X  10”^®  cm^/s  found  in  this  work  for  the  reaction  of  0^117'*’  formed  from 
norbornadiene  with  toluene  (see  Table  30).  Although  a  definitive  structure 
can  not  be  assigned  to  CyHy"*"  formed  from  VIII,  one  possible  structure  is  the 
norcaradienyl  ion,  IX,  which  was  calculated  as  a  stable  structure  with  AHf  = 
235.9  kcal/mol  (107). 
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Supporting  experimental  evidence  for  the  existence  of  stable  C7H7'*’ 
tures  other  than  benzyl  and  tropylium  comes  from  a  UV  absorption  study  of 
ions  in  solid  argon  (129)  •  A  trend  in  reactivity  was  observed  for 
reactions  of  CyHy'*'  ions  from  norbornadiene  and  cycloheptatriene  with  precursor 
neutrals,  acetylene,  and  diacetylene.  In  all  three  cases,  the  reactive 
structure  formed  from  norbornadiene  reacted  more  slowly  than  that  formed  from 
cycloheptatriene.  This  is  consistent  with  the  general  observation  that  less 
reactive  structures  of  small  hydrocarbon  ions  have  lower  heats  of  formation 
since  structure  formed  from  VIII  has  a  higher  heat  of  formation  than 

that  formed  from  V. 

Finally,  in  relation  to  the  ionic  soot  formation  mechanism,  the  existence 
of  a  different  Ion  less  stable  than  benzyl  in  flames  would  account  for 

the  low  CyHy"*"  concentration  found  In  flames  (14a)  .  As  with  and 

(42),  also  reacts  rapidly  with  diacetylene,  but  does  not  react  with 

acetylene.  Furthermore,  in  this  study,  product  ions  with  increased  C:H  ratios 
having  as  many  as  13  carbon  atoms  have  been  observed  in  the  reaction  of  CyHy'*’ 
with  diacetylene.  There  is,  however,  still  the  possibility  of  stabilization 
of  a  larger  fraction  of  the  Cgllg'*'  reaction  complexes  by  the  greater  number  of 
third-body  collisions  in  atmospheric  pressure  flames.  In  fact,  such 
collisional  stabilization  for  C5H5'*’  and  reaction  complexes  has  been 

observed  in  higher  pressure  SIFT  studies  (103,104). 
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SECTION  IV 


KINETIC  MODELING  OF  THE  REACTIONS  OF  C3H3'''  WITH  ACETYLENE, 
DEUTEROACETYLENE,  AND  DIACETYLENE 

A.  INTRODUCTION 

The  energetics  of  small  hydrocarbon  ions  such  as  C3H3'*’,  C^H3'*’,  C^H5'*’>  and 
CyHy'*’  and  their  reactions  with  a  variety  of  neutral  hydrocarbons,  particularly 
acetylene  and  diacetylene,  have  been  reported  in  Sections  II  and  III. 

As  reported  in  Section  IIIC,  bimolecular  reactions  of  the  propargylium 
form  of  C2H3'^  with  acetylene  most  often  result  in  an  isomerization  to  the 
cyclopropenylium  isomer.  To  help  understand  this  isomerization  process,  C3H3'^ 
reactions  with  deuterated  acetylene  were  investigated  (42)  •  These  studies 
showed  that  the  isomerization  proceeds  via  the  C5H3'*’  ion/molecule  reaction 
complex,  which  is  sufficiently  long-lived  under  the  experimental  conditions 
employed  that  deuterium  exchange,  as  well  as  isomerization,  takes  place.  Thus 
with  time  the  reactive,  propargylium  C3H3'^  isomer  is  converted  to  both 
reactive  and  unreactive  species  containing  one,  two,  and  three  deuterium 
atoms. 

To  better  understand  the  isomerization  which  converts  the  reactive  to  the 
unreactive  form  of  C3H3'*’,  kinetic  modeling  studies  of  the  ion  Intensity  versus 
time  curves  reported  earlier  (42)  were  carried  out.  It  was  also  hoped  that 
fitting  procedures  would  produce  improved  ion/molecule  reaction  rate 
coefficients.  Quantum  mechanical  calculations  on  C2H3'*'  and  C^Hj"*"  structures 
and  reactivity,  reported  in  Sections  IID-IIF,  were  used  to  guide  the  modeling 
effort . 

B .  METHODOLOGY 

Nonlinear  least-squares  fitting  routines  employing  Marquardt's  algorithm 
(138),  implemented  on  two  different  computers^^,  were  used  for  kinetic 
modeling.  Complete  analytical  solutions  were  obtained  from  the  chemical 
models  developed  below  for  the  systems  C3H3'*’  +  C2H2  and  C3H3'*’  +  C^H2*  A 
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complete  analytical  solution  was  not  possible  when  an  isotope  effect  was 
included  in  the  chemical  model  for  the  C3H3''’  +  C2D2  system.  Numerical 
integrations  used  the  finite  difference  method  (138).  The  independent 
variable  step  size  for  the  numeric  portion  was  generally  set  to  0.001  s. 
Increasing  the  step  size  to  0.002  s  had  negligible  effect  on  the  results. 

C .  THEORY 


In  studying  reactions  it  is  desirable  to  have  knowledge  of  the  internal 
energy  distribution,  P(E),  and  the  microcanonical  rate  constants,  k(E).  With 
this  information  the  time-dependent  expression  for  the  jth  rate  constant, 
kj(t),  can  be  written  as: 

k  (t)  =  P(E)k  (E)exp{[-zk^(E)]t}dE,  (119) 

in  which  Eq  is  the  threshold  energy  for  the  jth  reaction  and  the  summation  in 
the  exponential  term  is  over  all  competing  reactions,  including  the  jth 

reaction  (139,140).  In  fast  unimolecular  ionic  decomposition  reactions  P(E) 
might  not  be  characterized  by  a  Boltzmann  distribution. 

The  reactions  studied  here  are  bimolecular  addition  reactions  followed  by 
unimolecular  decomposition.  Given  the  ionization  potential  of  Xe  (12.1  eV) 
and  the  reported  (44)  appearance  potential  of  the  propargylium  ion  from  C^H^I 
(10. 5  eV),  1.6  eV  of  excess  energy  is  available.  Much  of  this  excess  energy 
will  be  converted  into  translational  motion  of  the  heavy  Xe  and  I  neutrals 
resulting  from  the  charge  transfer  in  the  collision  process.  Under  typical 
experimental  conditions  (ca.  2  -  3  X  10“^  torr  total  pressure) ,  about  125  ms 
was  allowed  for  the  charge  transfer  process  and  for  ejecting  intermediates. 

For  these  conditions  the  ions  collided  a  number  (10-15)  of  times  with 

the  excess  Xe  present  in  the  FTICR  cell,  leading  to  near  thermalization  of 

internal  energy  before  the  ion/molecule  reactions  were  monitored.  Since  the 
reaction  time  scale  was  on  the  order  of  seconds,  we  believe  that  any  slight 
initial  deviation  from  Boltzman  behavior  presents  no  serious  error. 

The  long  reaction  time,  due  to  the  bimolecular  nature  of  the  process, 
allows  for  energy  randomization  very  early  within  the  time  scale  of  the 
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reactions.  Therefore  simple  rate  constants  and  rate  expressions  can  be  used 
in  Lieu  of  detailed  applications  of  Equation  (119). 


D.  RESULTS 


1.  Models  of  C2H3'^  +  ^2^2 


As  reported  earlier  (42),  collision  of  the  propargylium  cation,  i- 
C3H3'*’,  with  acetylene  forms  the  cyclopropylium  cation,  C-C3H3'*’,  which  is 
unreactive  on  the  time  side  of  the  experiments,  given  the  pressures  attainable 
in  the  FTICR  cell.  Experiments  with  C2D2  showed  that  an  encounter  complex 
which  allows  for  isotopic  scrambling  is  formed.  Hence,  whatever  the  isomeric 
form  of  this  complex,  a  structure  having  the  chemical  formula  C5H3'''  can  be 
postulated.  Since  no  species  of  m/z  65  is  observed  in  the  mass  spectrum,  the 
species  must  be  in  steady-state  and  of  low  concentration.  The  simplest 
scheme  which  takes  into  account  this  information  is: 


2,— C3H3"^  +  ^2.^2  — k~^ 

C3H3I  — sink 


->  i-C3H3'^  +  C2H2 

c— C3H2  4"  ^2^2 


(120) 


Applying  steady-state  kinetics  to  scheme  (120)  yields 


I(t)  -  1°  -  kpl°  (1  -  e"®'^)/9. 


(121) 


in  which  1°  and  I(t)  are  ion  intensities  initially  and  as  a  function  of  time. 


respectively,  refers  only  to  the  linear  form,  k'  *  h  I* 

®  ‘‘cU  *  '‘p'- 


(122) 


Equation  (121)  was  fitted  to  several  kinetic  runs  reported  earlier  (42). 

Table  32  shows  results  of  these  fits.  A  plot  of  9  versus  P_  „  should  be 

^2^2 

linear,  as  Implied  by  Equation  (122),  and  this  is  demonstrated  In  Figure  41. 
Results  yield  kfkg/(k^  +  k^,)  -  (2.3  ±  0-2)  x  10^  torr~l  s'^  and  k'  -  (1.4  ± 
0.2)  s"’-. 
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TABLE  32.  RESULTS  OF  FITS  OF  EQUATION  (121)  TO  KINETIC  DATA  FOR  THE  C3H3'''  + 
C2H2  REACTION.^ 


P  /torr 

^2^2 

k'  I°/Arb  units  •  s  ^ 

^  1 

9/s~^ 

0 

1.53( .06) 

1.58( .08) 

4.4  X 

1.65(  .20) 

2.40(  .31) 

1.0  X  10“^ 

1.74( .12) 

3.60( .24) 

1.6  X  10’^ 

1.71(  .08) 

4.82(  .23) 

2.0  X  10"^ 

L.71(.ll) 

6.33(.38) 

^  The  standard  error  of  estimate  computed  by  the  fitting  program  is  shown  in 
parentheses. 


2.  Models  of  C3H3‘*‘  +  C2D2. 

Reaction  of  C3H3'*’  ^2^2  complicated  by  the  observation  that 

isotopic  scrambling  occurs  and  isotope  effects  are  possible.  Several  models 
have  been  tried,  including  those  which  allow  for  complete  isotopic  scrambling 
and  those  which  allow  only  partial  scrambling.  The  simplest  model  allows  for 
complete  scrambling  and  no  isotope  effects^^. 

Complete  scrambling  occurs  when  fragmentation  of  the  complex,  C5L5^  (L  * 
H,  D) ,  yields  precursors  having  a  statistical  distribution  of  hydrogen  and 
deuterium  atoms.  C5H3D2^  will  then  yield  the  following  ratios  of 
precursors;  C3H3'^:  031120'^  :C3HD2''’  *  l.'6:3.  Using  these  statistics  for 

obtaining  the  isotopic  distribution  in  the  C3L3'*’  precursors,  we  find: 

kf 

i-C3H3'*'  +  C2D2  - C5H3D2‘'' 

k 

Z-C3H3'*'  +  C3H3I  — sink 
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+  10  ‘^2°2  T  * 

T  S'®  +  To  '-S'^z"  ^  W  SS 


SSS^  ^  TO  ^  TO  SS  T  '-'=3S'>'‘  ^ 


I  S'®  *  TO  ‘=-S‘®2^  TO  SS 


f 

l-C^^2^'^  ■*■  ^2°2  - <^5fi2°3^ 


)1-C3H2D'''  +  C3H3I  — ^->  sink 


sss*  — >  TO  ’■-^iK  *  TO  SS  ♦  T  '-'=3“s^  ^ 


f  S'®  *  TO  ‘-SS”'’  TO  SS 


SSS*  — >  TO  '-"sS*  *  TO  SS  T  '-'3‘®2* 


T  S'®  *  TO  "-SS”"  TO  S”2 


CV£ 

a-C3HD2'’'  +  C2D2  - >  G5HD4''' 

+ 

Z-C3HD2  +  C3H3I  — —>  sink 


C.HD,'^  — —>  i  i-C.HD-'*’  +  i  C-D-  +  i  i-C-D,'*’  +  C-HD 
54  5  32  5225  33  52 

k 

— 1->  C-C2HD2'’'  +  3-  C2D2  -••  y  ^-<^303^  +  y  C2HD 


2-0303*  3-  C3O3  -i->  C3D5* 


l-C  0*^  +  0  H»I  — —>  sink 
3  3  3  3 
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•V 


S“5  — ^  '-^“3'’  *  S“2 

+  ^  + 

S°5  — ^  ■"  ‘=2‘>2 


(123) 


With  the  assumption  again  that  all  four  C5L5'^  complexes  are  in  steady-state  a 
full  analytical  solution  is  possible  for  the  set  of  kinetic  differential  rate 
equations.  The  solutions  to  f  e  set  of  equations  are  given  in  Appendix  B. 
Figure  42  shows  the  best  fit  curves  to  a  typical  data  set.  Table  33  shows  the 
fitted  parameters,  errors,  and  residual  sum  of  squares  from  fits  of  a  typical 
data  set:  these  results  are  labelled  F  »  1. 


An  examination  of  Figure  42  shows  that  the  model  given  in  scheme  (123) 
does  not  adequately  explain  the  production  of  C2HD2^*  attempt  to 

examine  this,  incomplete,  or  partial  Isotopic  scrambling  was  next  assumed.  To 
do  this  correctly  requires  a  detailed  knowledge  of  the  chemistry  of  the 

system,  which  is  not  available.  A  somewhat  crude  application  of  isotope 
effects  applies  mtiltipllcatlve  factors  to  the  individual  rate  constants  and 
this  procedure  requires  but  a  nominal  knowledge  of  the  structure  of  the 

species  lnvolved^5,  apply  this  simple  scheme  as  described  below. 


Quantum  mechanical  calculations^'^  Indicate  that  ^2^2  does  not  react  with 
the  cyclopropenyliun  cation,  but  does  react  with  the  propargylium  cation  (£- 
^3^3^^  without  barrier  with  formation  of  the  following  four  possible  products: 
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Figure  42.  Model  Fit  (scheme  (123))  to  Typical  Data  Set  for  C3H3'*’  +  ^2^2 
Reactions. 
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TABLE  33.  RESULTS  OF  MODEL  FITS  (SCHEMES  123  AND  124)  FOR  THE  SYSTEM  €3113'^  +  C2U2  UNDER  VARIOUS  EXPERIMENTAL 
CONDITIONS . 
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Theraodynamically ,  only  I  is  stable  with  respect  to  decomposition  to  C-C3H3''’  + 
02^1-7  •  However,  since  the  Ji-C3H3'*’  +  C2H2  -*■  I  reaction  is  at  least  60  kcal/mol 
exothermic,^  in  the  absence  of  stabilizing  collisions  this  energy  then  permits 
many  different  isomeric  forms  to  be  sampled  before  decomposition  back  to  C3H3'*' 
+  C2H2*  If  fbe  sampling  of  all  isomers  is  fast,  complete  isotopic  scrambling 
is  expected.  Given  the  uncertainty  in  structures,  total  equivalence  of 

carbon  atoms  in  the  complex  is  assumed.  For  ease  of  understanding,  a  cyclic 
C5H5'*’  complex,  in  which  all  C  atoms  are  sp^  hybridized,  might  be  visualized. 
With  these  assumptions,  it  is  a  simple  matter  to  determine  which  carbon  atoms 
undergo  hybridization  changes  during  the  reaction.  For  instance,  in  the 

attack  of  C2H2  on  il-C3H3^*  three  carbon  atoms  change  from  sp  to  sp^ 

hybridization.  Upon  fragmentation,  some  carbon  atoms  remain  sp^  hybridized; 
others  become  sp  hybridized. 

In  general,  for  a-secondary  Isotope  effects,  an  increase  in  hybridization 
in  going  from  the  reactant  state  to  the  transition  state  yields  an  inverse 
isotope  effect  (k^  <  k^) ,  whereas  a  decrease  in  hybridization  yields  a  normal 
isotope  effect  (kfj  >  kp)  (142).  la  our  system  the  bimolecular  addition 
reaction  will  have  an  inverse  isotope  effect,  and  the  fragmentation  reaction 
should  have  a  normal  isotope  effect.  Figure  43  shows  the  energetics  for  the 
reaction  of  £-C3H3^  ^2^2  isotope  effects  for  the  elementary 

reactions  can  arise.  For  the  forward  addition  reaction,  Epjj  >  Ep^  and  for  the 
fragmentation  reaction,  Ejyj  <  Ej^q,  in  which  the  subscripts  ”F"  and  "R"  refer 
to  forward  and  reverse,  respectively.  Since  the  bimolecular  addition  is  very 
exothermic  for  the  formation  of  most  C3L5^  isomers®’^  it  might  be  expected 
that  the  "average"  transition  state  structure  might  closely  resemble  the 

reactants  and  not  any  of  the  C3L5'*’  isomers.  This  in  turn  implies  that  Ej^^  - 

^  ^FH  “  ^FD  (zero-point  effects).  However,  the  excess  energy  in  the 

reaction  will  allow  longer  sampling  times  for  the  more  energetic 
isomers.  The  less  energetic  isomers  which  are  sampled  will  undoubtedly  be  in 
higher  rotational  and  vibrational  levels.  The  overall  effect  is  to  lessen  the 
normal  isotope  effects  expected  in  the  fragmentation  of  the  C5L5'*'  Isomers, 
unless  there  are  sufficient  collisions  to  stabilize  the  isomers  prior  to 
fragmentation. 
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+  3J 


Figure  43.  A  Schematic  Diagram  of  the  Reaction  of  £-€3113'’’  with  02^2'  The 
values  given  in  the  Figure  are  kcal/mol.  Since  experiments 
indicate  that  C-C3H3'*’  can  be  found,  cannot  be  much  greater 

than  28  kcal/mol.  AG2*  Is  unknown  as  the  studied  reaction  leads 
to  no  stable  C5H5'*’  species  In  the  absence  of  colislons. 
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In  the  model  development  outlined  below,  only  o-secondary  isotope  effects 
will  be  considered  important.  B-secondary  isotope  effects,  arising 
predominantly  from  hyperconjugation,  can  sometimes  be  important  (143),  but 
will  be  assumed  here  to  be  minor  compared  to  the  cr-effects.  It  is  also 
possible  that  hydride/deuteride  shifts  may  be  occurring  in  the  transition- 
state.  However,  hydride  transfers  often  exhibit  small  isotope  effects  (144), 
and  hydride  shifts,  if  they  occur  at  all,  will  be  assumed  here  to  give 
negligible  contributions  to  the  isotope  effects. 

Even  though  the  different  reactions  represented  by  Equations  (113)  will 
have  different  isotope  effects,  the  introduction  of  an  independent  fitting 
parameter  for  each  reaction  is  not  justified.  Though  crude,  only  one 
additional  parameter  is  Introduced  into  the  model  to  account  for  all  potential 
a-secondary  isotope  effects.  The  method  for  introducing  this  parameter  is 
outlined  as  follows. 

If  a  carbon  atom  bearing  a  deuterium  atom  undergoes  a  hybridization  change 
from  sp  to  sp2  (force  field  becoming  stronger),  the  "isotope  effect  factor", 
F,  is  introduced  as  a  multiplicative  factor  in  the  rate  constant.  For  two 
deuterium  atoms,  F^  is  the  multiplicative  factor,  etc.  If  the  deuterium  atom 
is  attached  to  an  atom  changing  from  sp^  to  sp  hybridized,  the  rate  constant 
is  divided  by  F,  for  two  deuteriums,  F^,  etc.  Introducing  the  same  factor  for 
both  addition  and  fragmentation  reactions  implies  a  constraint  which  is  at 
best  only  qualitatively  correct.  Applied  in  the  numerator  the  factor  corrects 
for  a  single  deuterium  atom  (F^  for  two,  etc.)  attached  to  a  center  undergoing 
hybridization  change  from  sp  to  sp2  in  the  transition-state  complex.  Applied 
in  the  denominator,  it  corrects  for  a  change  from  sp^  to  sp  hybridization. 
Since  the  transition-state  complex  has  a  stronger  force  field  than  the 
reactant  state  (^-€213"^  +  C2L2)  at  these  centers,  the  ct-secondary  isotope 
effect  should  be  "inverse"  and  F  should  therefore  be  greater  than  unity. 

Using  the  structural  notation, 

-  [  [J ^  C-C-C-D  , 

I-HDC^d'^  =  [  p^C-C-C-D]'*’,  and  °  ^  C«C-C-H  , 
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the  full  kinetic  scheme  is: 


i> 


I 

+  G^H^I  — L.>  sink, 
k 

c  H  d/  — i->  +  4,  i-nwy  +  j-  + 

^  ^  ^  lOF^  3  3  5F  3  gp2  2  3 

5F  ‘-“^3“^  ^  TO  '-“2'=3“* 


S«3“2''  — >  ^  -=-<=3“3''  ^  If  ‘=-«3“2''''  ^  TO  '=^3““2* 

F^k 


^■^2^3°^  +  C2D2  - ^>  C5H2D3'^ 


i-HDC^H'*’  +  C3H3I  — 2->  sink 


t-H^CjD  +  C3H3I  — L.>  sink 

^5^2°3  W  ■*■  i-HDCjD'*’  +  i-D^C^n'*’  + 

5F 

a-HDC  h'*'  +  e-H  c  d'*’ 

5F^  ->  lOF^  ^  ^ 

S‘'2“3^  — >  TO  ‘=-<=3‘>3''  *  If  '-'3'®2''  *  7^  "^3“2'>* 

^  lOF 

F-^k 

I-HDC.d'^  +  C-D-  - ->  C.Hd/ 

3  2  2  5  4 


£-D2C3H  +  - i-> 


Z-HDC^d’'"  +  C3H3I  - sink 


Z-D^C^h'^  +  C3H3I  — sink 


CcHD  — —>  ^-C,D  £-HDC  H-D-C.r'*' 

5F^  ->  3  J  3p2  2  3 
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^c  ,  2  . 

HD,  - 

— >5f"'S°3 

f\ 

^  ^f  + 

+  >  C  n. 

3  3 

2  2  5  5 

■C  D 

3  3 

Fk 

J_  P  U  T  P  N 

3  3 

5F^ 


k.  /F^ 

CjDj-"  ^_>  C-O3D3 


(124) 


The  neutral  fragmentation  products,  C2H2.  C2HD,  and  C2D2,  have  not  been 
included  for  brevity.  It  is  assumed  that  i-C3L3‘'’  undergoes  the  same  kind  of 
hybridization  changes  when  reacting  with  C3H3r  as  it  does  with  C2L2*  Earlier 
work  (42,76)  has  shown  that  the  products  of  the  +  C3H3I  reaction  are 
of  higher  mass  and  do  not  enter  further  into  the  kinetic  schemes  modeled  here. 


Steady-state  conditions  are  applied  for  all  C5L5''’  isotopic  species  as 
before  (without  regard  for  isomeric  differentiation)  and  the  differential  rate 
Equations  which  represent  Equations  (124)  are  given  in  Appendix  B.  Table  33 
shows  fitting  results  for  several  data  sets  using  Equations  (124).  Figure  44 
shows  plots  of  the  best  fit  of  this  model  to  the  same  data  set  as  fitted  in 
Figure  42. 


Using  FTICR  double  resonance  techniques  (69)  some  experiments  were 
conducted  where  certain  ions  were  ejected  from  the  analyzer  cell  as  they 
formed.  Among  the  ions  ejected  were  C3H2D'*'  and  €31102'^.  Without  further 
fitting,  this  model  was  used  to  predict  the  behavior  of  the  kinetic  system  if 
there  were  ejected.  Figure  45  shows  predicted  results  and  data  points. 


3*  Models  for  €3113'^  +  C^H2. 

The  reaction  of  C3H3'*'  with  C4H2  (diacetylene)  is  kinetically  more 
complicated  than  that  of  C3H3'^  with  C2H2.  There  are  more  isomeric 
possibilities,  and  ion/ molecule  reaction  products  of  higher  m/z  are 
detected.  Several  models  were  tried  in  attempts  to  fit  the  experimental  data, 
and  the  best  of  these  made  the  assumptions  that  ^2115^  and  the  excited  forms  of 
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Figure  45.  Data  from  Ejection  Studies  (€3113^  +  C2D2  reaction)  and 
Prediction  (C3HD2'*'  ejected  -  top;  C3H2D'*'  ejected  -  bottom). 


Model 
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C9H5'''  steady  state,  while  CyH^"*"  ^re  not.  In  order 

to  account  for  Che  build-up  of  CgH^"*"  and  ions,  it  was  necessary  to 

include  collision-induced  stabilization  steps  for  the  excited  forms  of  these 
ions.  A  further  assumption  was  made  that  the  stabilized  forms  of  and 

are  not  reactive  within  Che  time  frame  (2s)  of  the  experiments. 


The  resulting  scheme  is: 

k, 

+  C.H.  — —> 

3  3  4  2  7  5 


+  C3H3I 


+  '^-12  + 

C7H5  — -h 

+  ^-13  + 

C_H.  — c-C-H-  +  C.H., 
7  5  3  3  4  2 


+  f2  + 

+  C.H.  -^>  C^U. 

5  3  4  2  9  5 


+  C-H,I  -HL.>  C-H,!'*' 
5  3  3  3  00 


+  -21  + 


+  -22  + 

C.H.  — C-H-  +  C,H_ 
9  5  5  3  4  2 


(stabilized) 


C^„3+  +  -«-> 


CjHj*  *  W  =10"6^^ 


— >  W  *  Vz 
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(stabilized) 

k 

W*  — >  ^ 

k 

CgHg'*'  - ^>  sink  (125) 


Reactions  involving  CgHgl"*"  are  included  because  CgHg'*'  is  detected  in  the  ICR 
experiments.  Applying  steady-state  assumptions  to  the 

and  ions  yields  a  model  with  a  full  analytical  solution  (Appendix  B)  . 

Figure  46  shows  results  of  the  model  fit  to  a  typical  data  set.  The  kinetic 
parameters  have  been  grouped  as  exponential  and  pre-exponential  terms  in  the 
equations,  yielding  9  fitting  parameters.  Table  34  shows  results  of  model 


fits  of  two  data  seta.  The  term  kp2Pc  h  I  *^-21'l*2  *^l't*2^» 

hence  can  be  calculated  from  the  parame^e^s.  The  results  of  the  first  data 


fit  yield  in  Table  34  yield  h  I  ” 

second  data  fit,  k  P  „  -  (1.3  ±^.^)  s“^ . 

p2  C^H^I 


(1.2  ±  0.3)  3~  and  those  from  the 


E.  DISCUSSION 


1 .  ^3^3*^  ^2^2^ C2D2 

Several  insights  into  the  mechanism  of  the  CgH-j'*'  +  C2H2/C2D2  reactions 
are  obtained  from  these  modeling  studies.  First,  a  value  of  (1.4  ±  0.2)  s~^ 
Is  obtained  for  k^  from  the  intercept  of  Figure  41  (values  for  k^  are  a  little 
higher  from  model  fits  shown  in  Table  33).  Assuming  a  cell  temperature  of  363 
K  (3)  and  given  that  P_  »  3.3  X  10“®  torr  (after  system  factor 
correction),  kp  »  (1.6  ±  0*2)^X  10”^  cm^  *8"^ ,  in  good  agreement  with  previous 
work  in  these  laboratories.  Second,  the  model  gives  branching  ratios  for  the 
complex  C5L5'*'.  The  ratio,  k^^/k^,  is  ca.  0.34  (.03)  (Table  33),  implying  that 
CgL5''’  fragments  to  l-CgHj"'’  at  a  rate  three  times  that  of  fragmentation  to  c- 
^3^3^*  c-C'jHg''’  is  thermodynamically  more  stable  than  25  kcal 
mol“^  (16);  however,  only  I  can  decompose  to  c-C-jH-j^.  Apparently,  there  is 
enough  "sampling  time"  for  I  to  allow  c-C^H^'*’  to  be  formed  in  significant 
amounts.  Any  factor  affecting  the  internal  energy  of  C5L5'^»  such  a 
collisional  stabilization,  will  undoubtedly  affect  the  ratio  k(,/kj^.  Hence, 
k^/kjj^  is  probably  very  sensitive  to  experimental  conditions  (method  of 


TABLE  34.  ElESULTS  OF  MODEL  FITS  (SCHEME  125)  FOR  THE  SYSTEM  CoH-i  +  C.Ho  UNDER  VARIOUS  EXPERIMENTAL  CONDITIONS. 
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ionization,  total  pressure,  and  temperature).  Finally  the  "direction"  and 
magnitude  of  F  (1.67  (.10))  suggest  not  only  that  the  model  is  reasonable,  but 
that  the  isotope  effect  plays  an  important  part  in  the  overall  kinetics. 

It  should  be  remembered  that  F  is  not  an  overall  isotope  effect;  isotope 
effects  are  ratios  of  rate  constants.  The  magnitude  of  F  does,  however, 
indicate  the  kinetic  contribution  of  replacing  one  hydrogen  atom  with  a 
deuterium  atom  on  a  hydrogenic  site  undergoing  a  hybridization  change  in  the 
transition-state  structure.  F~^  is  the  ct-secondary  kinetic  isotope  effect  for 
the  elementary  bimolecular  addition  reactions  per  hydrogenic  site  (only  those 
sites  giving  rise  to  a-secondary  effects  are  included).  F  is  the  isotope 
effect  for  the  elementary  fragmentation  reactions  (to  and  C-C2H3''’)  ?er 
hydrogenic  site  (again  only  those  sites  giving  rise  to  ct-secondary  effects  are 
Included).  Because  of  the  complexity  of  scheme  (124),  overall  isotope  effects 
cannot  be  gleaned  by  determining  product  ratios,  as  is  often  done  in  isotopic 
studies  of  more  elementary  ion/molecule  reactions.  However,  F  can  be  compared 
to  isotope  effects  found  in  elementary  ion/molecule  reactions. 

oc-secondary  kinetic  isotope  effects  can  be  quite  large  (145,146).  The 
larger  effects  found  in  ion/molec\ile  reactions  are  due  to  a  narrow  (non- 
Boltzmann)  energy  distribution  centered  close  to  the  threshold  energy  (147). 
Even  though  the  secondary  isotope  effects  for  the  elementary  reactions  are 
fairly  large  in  our  system,  they  are  not  nearly  as  large  as  in  some  of  the 
other  systems  studied  at  low  pressure  in  ICR  experiments.  The  magnitude  of 
the  isotope  effect  factor,  F,  indicates  a  fair  amount  of  bond  rearrangement  in 
the  transition  state  structures.  Due  to  the  number  of  assumptions  and  the 
Indirect  method  for  obtaining  F,  we  have  not  attempted  to  calculate  Isotope 
effects  from  molecular  models  and  calculated  vibrational  frequencies. 

2 .  C3H3'*'  +  C^H2 

Kinetic  data  and  modeling  of  this  system  reveal  some  insight  into  the 
mechanism  for  the  reaction  of  C3H3'*’  and  C^H2»  Although  scheme  (125)  does  not 
account  for  specific  Isomeric  reactivities,  the  scheme  does  account  reasonably 
well  for  the  data.  The  signal-to-noise  ratios  of  the  available  data  do  not 
warrant  further  refinements  of  the  model.  Several  crude  models  were  at- 
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tempted,  and  the  model  reported  here  may  not  be  the  best  or  only  model  which 
can  adequately  explain  the  results.  However,  of  the  variations  of  scheme 
(125)  which  were  tried,  the  one  reported  here  gave  the  best  fitting  results. 
In  particular,  a  scheme  was  tried  in  which  stable  forms  of  CgH5''’  and 
were  formed  from  the  corresponding  bimolecular  additions  reactions  rather  than 
requiring  stabilization  steps  of  "hot"  and  ions.  This  scheme  did 

not  fit  the  data.  Apparently,  only  in  the  ions  larger  than  are  there 

enough  internal  vibrational  modes  to  give , rise  to  an  ion/molecule  collision 
complex  of  sufficient  lifetime  to  allow  collisional  stabilization  at  these 
pressures . 

Scheme  (125)  does  not  show  and  CyH3''’  fragmenting  to  any  products. 

Undoubtedly  any  isomers  of  and  CyH^"*"  which  form  are  all  resonance- 

stabilized.  Any  fragmentation  product  would  be  highly  energetic  and  unlikely 
to  form,  even  at  the  energy  levels  available  in  these  reacting  systems. 
Hence,  once  C3H3^  and  ions  are  formed,  they  will  react  only  in  bimolecu¬ 

lar  addition  reactions,  which  are  orders  of  magnitude  slower  than  the  uni- 
molecular  fragmentation  reactions  of  CgH5'*’>  and  hence,  there 

Is  a  non- steady- state  build  up  of  C5H3‘'’  and  ions  at  the  first  part  of 

the  reaction. 


SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  experimental  and  theoretical  work.  reported  here  contributes 
significantly  to  an  understanding  of  the  reactivity,  energetics,  and 
structures  of  small  hydrocarbon  ions.  However,  it  does  not  prove  or  disprove 
an  ionic  mechanism  of  soot  formation.  It  does  lead  to  several  clarifications 
of  the  ionic  mechanism,  which  will  be  discussed  below. 

The  low  bimolecular  reactivity  of  acetylene  with  all  of  the  small 
hydrocarbon  ions  studied  indicates  that  it  may  not  be  as  important  as 
postulated  in  soot  nucleation.  On  the  other  hand,  as  reported  in  Sections 
IIIB,  IIIC,  and  HID,  dlacetylene  reacts  at  rates  approaching  the  Langevin 
limit  with  all  of  the  small  hydrocarbon  ions  studied,  producing  larger  ions 
with  high  carbon- to-hydrogen  ratios.  The  general  reaction 

CjjHy'*'  +  C4H2  Cx+2Hy‘^  +  C2H2  (126) 

appears  to  occur  in  a  facile  manner  and  could  be  quite  important  in  soot 
nucleation. 

The  experimental  studies  carried  out  here  were  at  significantly  lower 
temperatures  and  pressures  than  those  found  in  most  combustion  environments. 
Studies  of  selected  rate  constants  over  a  limited  range  found  no  significant 
temperature  dependence.  At  high  pressures,  significant  colllsional 

stabilization  of  ion/molecule  collision  complexes  can  occur.  Thus,  the 
complexes  formed  by  reaction  of  acetylene  with  small  hydrocarbon  ions  might  be 
stabilized,  leading  to  significant  termolecular  reactivity,  even  though  no 
bimolecular  reactivity  was  observed  in  the  work  reported  here. 

B.  RECOMMENDATIONS 

Continuation  of  this  work  in  several  directions  could  have  a  major  impact 
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on  assessing  the  importance  of  an  ionic  mechanism  of  soot  formation.  First, 
experiments  should  be  carried  out  at  higher  pressures,  using  a  high  pressure 
source  for  the  FTICR  mass  spectrometer,  to  assess  the  importance  of  third-body 
stabilization  of  ion/molecule  collision  complexes.  Actual  sampling  of  ions 
from  flames  into  the  FTICR  mass  spectrometer,  in  order  to  probe  their 
structures  and  reactivity,  would  extend  the  isolated  studies  done  here  to  an 
actual  combustion  environment. 

A  crucial  argument  in  favor  of  an  ion/ molecule  soot  nucleation  route  is 
the  ease  with  which  cyclic  ions,  which  could  lead  to  polycyclic  aromatic 
hydrocarbons,  are  formed.  Yet  the  work  reported  here  did  not  identify 
definitively  when  the  first  cyclic  ions  were  formed.  The  predicted  spectra 
from  theoretical  calculations  on  ions  were  obtained  after  the  end  of  the 

contract  period,  and  similar  data  which  can  be  calculated  using  the  same 
methodology  for  ^5^3^  ions  should  be  combined  with  laser 

spectroscopic  studies  on  these  ions  to  determine  which,  if  any,  have  cyclic 
structures.  Also,  theoretical  studies  which  probe  the  geometries  of  the 
ion/ molecule  collision  complexes,  such  .s  those  carried  out  here  for  C3Hj'^  + 
acetylene,  should  be  extended  to  other  important  reactions,  e.g.  with 
diacetylene,  and  the  nature  of  the  proposed  chemi-ionization  reaction  which 
forms  the  C|jH3^  Ion  (see  the  Introduction  Section)  should  be  investigated. 
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APPEMDIX  A 


If  Internal  coordinates  y  are  used,  this  choice  of  =  I  for  cartesian 
coordinates  x  corresponds  to  Gy^  =  where  y  =  Bx  and  = 

(B'*'mb)~^-B~*'in  where  ^  is  an  arbitrary  (3N  -  6)  x  (3N  -  6)  matrix;  see,  for 
example,  E.  B.  Wilson,  J.  C.  Decius  and  P.  C.  Cross,  in:  Molecular 
Vibrations  (McGraw-Hill,  New  York,  1955). 

-  Quantum  Chemistry  Interactive  Program  Utility  (QUIPU),  Quantum  Theory 
Project,  Departments  of  Chemistry  and  Physics,  362  Williamson  Hall, 
University  of  Florida,  Gainesville,  FL  32611. 

n 

J.  S.  Binkley,  M.  J.  Frisch,  D.  J.  Defrees,  K.  Raghavachari,  H.  B. 
Whiteside,  H.  B.  Schlegel,  E.  M.  Fulder  and  J.  A.  Pople,  Carnegie  Mellon 
Univ.  Pg,  PA  15213. 

^  F.  Wiseman,  F.  Ozturk,  J.  R.  Eyler  and  M.  C.  Zerner,  see  Section  IV. 

^  F.  Ozturk,  M.  Moini,  F.  Brill,  J.  R.  Eyler,  T.  J.  Buckley,  S.  G.  Lias  and 
P.  J.  Ausloos,  submitted. 

^  A.  Cameron,  J.  Leszczynski,  M.  C.  Zerner  and  B.  Weiner,  submitted. 

^  F.  L.  Wiseman,  F.  Ozturk,  M.  C.  Zerner  and  B.  Weiner,  submitted. 

®  J.  Feng,  J.  Leszczynski  and  M.  C.  Zerner,  submitted, 

q 

J.  Leszczynski,  M.  C.  Zerner  and  F.  Wiseman,  submitted. 

The  uncertainty  given  is  the  95  percent  confidence  limit  of  the  mean  of  5 
determinations.  Reported  uncertainties  for  the  correction  factor  and  rate 
coefficients  throughout  this  Section  are  also  95  percent  confidence 
limits. 

The  large  uncertainties  reported  for  the  correction  factor  and  rate  co¬ 
efficients  calculated  using  it  in  this  work  are  primarily  due  to  the  wide 
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range  of  reported  values  for  the  C2H4'*'  +  *"2^4  r'eaction  used  to  determine 
the  correction  factor.  The  even  more  widely  studied  "standard"  calibra¬ 
tion  reaction  of  with  CH^  could  not  be  employed  because  the  lowest 

mass-to-charge  ratio  accessible  to  the  FTMS-1000  Instrument  is  m/z  18. 

1 

Bartmess,  J.  E.  private  communication. 

Tektronix  Model  4054  and  Hewlett  Packard  Model  150. 

Lampe  and  Field  (F.  W.  Lampe,  and  F.  H.  Field,  Tetrahedron,  7,  189  (1959)) 
studied  the  reaction  of  CD^'*'  +  C2H4,  and  observed  the  following  yields: 
G3HD4''‘  :  1/10,  C3H3D2'^  :  "■'/5,  C3H,D3''‘  :  2/5,  and  C3HD4'‘’  :  1/10  (+L.  + 
L2).  Statistical  yields  wxth  no  Isotope  effect  would  have  been:  C3H4D'*’: 
1/14,  C3H5D2^  •  3/7,  C3H2D3'''  :  3/7,  and  C3HD4'''  :  1/14,  very  close  to  the 
observed  values.  Hence,  since  the  C3Lg'*’  complex  exhibited  almost  complete 
scrambling,  it  is  reasonable  to  expect  the  complex  modeled  here  to 

undergo  complete,  or  nearly  complete,  scrambling. 

Although  somewhat  crude,  applying  a  multiplicative  factor  to  the  rate 
constant  for  each  deuterated  site  is  in  keeping  with  experimental 
evidence.  For  instance,  in  the  acetolysis  of  some  identical  tosylates, 
each  deuterium  atom  substitution  changed  the  rate  constant  by  ca.  0.84 
(Reference  141). 
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APPENDIX  B 


ANALYTICAL  EXPRESSIONS  FOR  KINETIC  MODELING  IN  SECTION  IV. 


1.  Scheme  (123)  in  Section  IV. 


I  ,  =  1°  +  8<t)/(lO0  )]  [exp(9.  t)  -  l] 

^  +  [exp(9  t)  -  exp  (9,t)]  -  ^  64,1°  {[exp(9  t) 

C2H2D  i-C^H^  Z-Cyl^ 


-l]/9j^  -  3[exp(92t)  -  l]/92} 


I  +  “  +[exp(9.  t)  -  2exp(6»t)  +  eKp(9-t)]  + 

C3HD2  .-C3H3 


LO 


64.1 


Z-C3H3 


_^{[exp(9^t)  -  l]/9^  -  6[exp(92t)  -  l]/92 


+  6[exp( 93 1)  -  1  ]/93  } 


S°3 


^”^3^3 


^[3  exp(92t)  +  exp(9^t)  -  exp(9^t)  -  3  exp(93t) 


-  3<Cl°  {[exp(9.  t)  -  1]/(109  )  -  9[exp(9»t)  - 

Z-CjHj 

1]/(1092)  '''  ~  1]/(593)  ~  [exp(9^t)  -  l]/9^}, 


in  which  6  “  ”  ^£^f^C2D2^^^ll  ^  *^c^‘  “  ^^^^2^2  ~  '*P^C3H3l, 

92  *  3k^4,/10  -  kfPc2D2  "  S’^C3H3l‘  ®3  *  3k.j^4,/5  -  -  kpP^^fj^j,  and  64  - 

kj^4,  -  kfPc  D  -  kpPc  H  I-  The  fitting  parameters  for  this  model  are 
o  2  2  3  3  Q 

I  ,,  B,  kfPp  n  ,  and  k_Pp  u  t*  ^  j.  fixed  at  its  experimental 

t-CjHj*  ‘  '=2“2  ’’  W  c^h/ 

value  if  the  initial  point  is  at  t  ■  0;  otherwise  it  is  also  a  fitting 
parameter. 


2.  Differential  Rate  Equations  for  Scheme  (124)  in  Section  IV. 


dl  “ 

91 

i-C3H3 

dl  "* 

31 

C-C3H3 

*^“^3^3 

dl  /dt 

-  -  (F^k' 

^-HDC3H 

.L 

FI 

i-H^C 

dl  /dt 

-  -FCF^k' 

Z-H2C3D 

FI 

il-H2C 

dl  ./dt 

-  381^2^1 

C-C3H2D 

/(lOF") 


+ 

Z-HDC^H  Z-HDC^H 


)/(5F^)  +  2<j,  I  /(5F) 

D  ^-C3H3 


^.  +  ()>2  + 
il-H2C3D  i-HDC3H 


)/(10F^)  +  4,  I  ./(5f2) 

y  £-C3H3 


A-HDC3H  1-^2^ 


FI  ^)/(2F)  + 


ll  +]/(5F) 


il-C3H3 


dl 


£-HDC3D 


./dt  -  -F(F%  +  k  )I  +  2^  (FI 

T  t  P  „T  J 


Z-HDC^D 


Z-HDC3D 


I  +)/(5F^)  +  24,2(1 


Z-D2C3H 


Z-HDC3H 


FI  ^.)/(5F^)  +  <j.^I  ^/(5F) 

t-H2C3D  i-C3H3 
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dl 


, /dt  =  -(F^k.  -t-  k  )1 

-■T  ID 


f  ■  "'p'"  +  + 

i-D^C^H  ^  i-HDC^D 


JI-D2C2H 


^)/(5F^)  + 


Jl-HDC^H 


FI  .  )/(5F)  +  <j,,I  /lO 

A  ft  r\  *  ^  ’ 


i-H2C3D 


dl  /dt  -  3S[<|).(FI  +  I 


C-C3HD2  ^  ^  Z-HDC^H^ 


i-H,C,D 


3  3 


dl 


i-CjDj 


/dt  -  [d./r^  -  f(fV  +  k')]i  +  2»,CFI 

^  ^  Pa^TN*  »%.  ' 


Z-HDC^D 


+  I 


^)/(5F^)  +  ^2(1  +  +  +)/(10F) 


dl 


C-C3D3 


^/dt  -  6[.^4l  ^/F"  +  2 ^3 (FI 


£-€303 


£-HDC3H  £-H2C3D' 


£-HDC3D 


I  .)/(5F)  +  ^  (I  +  FI  )/10, 

A  TT  '  ^  A  Tt  ’  A  1Y  TV* 


£-D2C3H 


2-HDC3H  £-H2C3D 


In  which  9  -  (j)j^/(10F  )  -  F  -  k^,  k^  -  ^f^C2D2’  S  “  S^C3H3l’  h  * 

(X^  +  8X2),  4>2  *  +  8X^),  *3  -  F^k^/(X^  +  8X5),  -  F^k'/Cl  +  F6), 

and  the  X's  are  given  as: 


Xj^  -  (1  +  6/F  +  3/F^)/10 
X2  -  (3  +  6/F  +  1/F^)/10 
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=  X2/F 

-  (3  +  2/F)/5F^) 
\5  =  (2  +  3/F)/(5F). 


The  fitting  parameters  for  this  model  are  I 
1°  is  fixed  if  the  initial  point  is  at 

C3R3 

fitting  parameter. 


0; 


3,  kp,  and  F. 
otherwise  it  is  also  a 


3. 


I 


I 


C3H3 


C5H3 


C7H3 


Analytical  Expressions  for  Scheme  (125)  in  Section  IV. 
»  1°  .  +  1°  .  [1  +  ][exp(9,  t)  -  1] 

C3H3  £-€3^3 

*  *^-11^1^°  +[exp(9j^t)  -  exp(  92t)  ]/9^2 

•  ^-11^-21  +[exp(9j^t)/9j^3  -  exp(92t)/923  + 


(1/923  -  l/9|^3)exp(93t)]/0j^2 


I  ^  *  ^-11*^31 ‘<*1<<’2^°  +{[exp(9j^t)  -  lj/9^  -  [exp(92C)  - 

(s)  Z-C3H3 


I  ^  *  ^-11*^-21*^82*^1  *^2 *^3^  [®5q>(  9j^t)  -  l]/(9j^9j^3) 

(3)  ^"S^3 

[exp(92t)  -  1]/(92923)  -  (^/923  “  l-/9i3)  [exp(  93t)  - 
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^  4.  =*  [exp(9,t)/9,  -  exp(9-t)/9»  +  (1/9, 


-l/9^)exp(-k2t)  ]/9^2> 


in  which  ((jj^  =  ^-12  ’^-13^’  ^2  '''  ^-22 

®3  “  ^f3^C,H.' ^  -31  '*'  ^-32  ''’  ^s2^’  -0j^,  -021  and  -0^  are  parameters  (1), 
(2),  and  (3),  in  t^at  order,  in  Equations  (125),  9^2  “  9l  “  92*  ®13  “  ~  ®3> 


023  =  02  “  931  9{  *  9^  +  k^,  and  02  “  ®2  "**  '^d*  fitting  parameters  are  9]^, 

®2’  ®3*  ^  +'‘^  ^  *  ^—11*^1^  +^^12’  2l'^2’  ^sl*^2’  ^s2'^3’ 


kj.  It  can  be  seen  that  k^2  “  “9^^  -  k_2]^<))2  ~  ^31*^2’  which  kp2  = 
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